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ABSTRACT 


The Naval Postgraduate School possesses software designed to perform thermal 
analysis of electronic components. At the core of this package is a model builder whose 
purpose is to generate a thermal model for use in steady state and transient thermal 
analyzers. The current version of the model builder requires excessive amounts of time 
for data input and model verification. This thesis describes the development of a model 
builder specifically designed to reduce the time required to model a printed circuit board 


containing up to four copper layers. 


TABLE OF CONTENTS 


LiceI NTRODUCTION 2.0... 09%. oy 2s 2 enn cece I 
I]. A JUSTIFICATION FOR THERMAL ANALYSIS 22...) See -~ 
A. RELIABILITY™ . oP re 5 
B. MATERIAL .......... 5. 5 see ss 0 eee y 


C. THERMAL FATIGUE AND CATASTROPHIC THERMAL FAILURE ..7 


Il. HEAT TRANSFER oe... cei rerereereenr reenter 9 
A. CONDUCTION .........26.5..2.3 2 22s = eee one 2 

I. General Equation of Heat Conduction >. 75a 10 

2. Single Plane Slab Modeling 2.2... 2 3 Sys 1] 

B. CONVECTION (00000 ee ce oe eo eee 12 

l. Electrothermal Analog ....:........ 255 13 

C. RADIATION 2.2... 0.000 See ae cs ee I3 

1. Transformation of the General Radiation Equation” ©. 2). 14 

IV. FINITE DIFFERENCE ANALYSIS .. 2. 18 
A. FUNDAMENTALS OF FINITE DIFFERENCE ANALYSIS ee 18 

I. Furst and Second Derivative Approximation 3... .. 23) I8 

V. THE MODEL BUILDER ..... 2): 2 26 
A. FEATURES .... 2.20.5. 50 mee © seine ents 26 

B. THE THERMAL ANALYZER INPUT DATA FILE ee 28 

C. PCB MODEL BUILDER SAMPLE PROBEEM ee 3] 

1. Printed Circuit Board Specifications  . . (eyes. 31 

2. Terminal Session ....06.-.-..02+ 000) 00g 32 

3. The output data file ... 2.2.90) eee 42 

4. Thermal analyzertoutput ..... 2.5...) ss 42 

D. PCB MODEL BUILDER VALIDATION (ee 42 

E. MODEL BUILDER LIMITATIONS AND POSSIBLE IMPROVEMENTS 45 

1. PCB Model Builder Limitations ....... . . 9) > eee 45 


Pee Oss ibicmampuomwements to the Thermal Model Builder. ............. 45 


ONCLUSIONS 2... ee ee ee eee 46 
APPENDIX PROGRAM LISTING FOR PCB MODEL BUILDER ......... 47 
Meer REFERENCES «2 .2.......5. EROS Cee eae ne eee a 145 
IDRIS SU eTOM GS) 1) 0S) tte re 146 


LIST OF TABLES 





Table 1. MATERIALS USED IN PCB PRODUCTION 


LIST OF FIGURES 


meee 1, 1ypical printed circuit board configuration. ....sm..........05.. 08 3 
Figure 2. Failure rate as a function of time for typical component. ............ 6 
Figure 3. Radiation, convection, and conduction electrothermal equivalent. .... 16 
ievies.4, First and’second derivative approximation.» 21... ..cc cee eve ee aes 19 
Pelice olrdphicalrepresentation of asingle node, .....s..cs.e0s.......008- 22 
Meeure 6. Node arrangement on typical printed circuit board. ............... 23 
iene). oainplic PCR partial output data file. .me.......... WW. ee eee 29 
meunes &. Printed circuit board modeled in the terminal session. ............. 33 
mieure 9. Partial output data file of PCB terminal session. ................. 43 
meu wOutputdata file of FASS thermal analyzer. ................00065 44 





I. INTRODUCTION 


The Naval Postgraduate School possesses software designed to perform thermal 
analysis of electronic components. The software package contains two routines de- 
signed to generate a thermal model (in the form of an ASCII input data file) to be read 
by the thermal analyzer program. The first routine is a generalized model builder used 
in developmental stages, as well as an editor used to modify existing models. The second 
routine was designed to work with models that have a specific geometric configuration. 

In order to accurately model electrical components, the structure must be subdivided 
into small, equal sized subvolumes. The centroid of each subvolume is referred to as a 
node, and due to the assumed isothermal nature of each subvolume, a node is considered 
as representative of the total subvolume. 

Producing an accurate thermal model, requires the design engineer to deal with 
enormous amounts of equations and temperatures, each describing an indiviual node. 
Additionally, each node is connected to adjacent nodes by thermal conductances. As the 
desired accuracy increases, the number of nodes and equations also increases. Modeling 
the electronic component without the aid of a computerized model builder 1s a task re- 
quiring inordinate amounts of time. 

Current versions of the thermal analyzer contain provisions for the generation of the 
node equations and node interrelationships; however, data input is still a manual, time- 
consuming process. At present there exists a need for a program that will generate a 
data-input file for the thermal analyzer that is both generated through a menu-driven 
routine, and allows the design engineer enough flexibility to model the electronic com- 
ponent to suit his or her needs. 

This thesis describes the development of a model builder designed specifically to re- 
duce the time required to model the copper and epoxy layers of a printed circuit board. 
A typical printed circuit board configuration is shown in Figure 1. The printed circuit 
board may contain up to four copper layers. Additional features of the model builder 


covered are: 
1. The capability of working in SI or English units. 
2. The choice of a total of sixteen aspect ratios. 
3. The provision for up to 740 nodes. 
4 


. The ability to specify the percent of copper coverage for each layer. 


2 - Saas 


The ability to specify dimensional parameters for each individual copper layer. 
The ability to input heat dissipation using several methods. 
The provision for six ambient temperatures. 


The automatic calculation of conductance values based on user input. 





EPOXY LAYER 
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Typical printed circuit board configuration. 


Figure 1. 


Source: Reference 2. 


Il. A JUSTIFICATION FOR THERMAL ANALYSIS 


A printed circuit board is a conglomeration of organic and inorganic materials with 
external and internal wiring, which allows electronic components to be physically sup- 
ported, and electrically connected. [Ref. 1] 

Over the past several decades, printed circuit board technology has developed sub- 
stantially. Early printed circuits were fabricated by printing a pattern of polymer resist 
on a copper plane and then chemically etching. Holes drilled in the laminate held the 
component leads that were soldered to the copper-printed patterns. The technology has 
progressed in developing the sophistication and uses of the printed circuit board inter- 
connections. 

Today, the basic functions of the printed circuit board are the same; the intercon- 
necting copper signal lines join two I/O leads from two different components. The 
components may be resistors, inductors, capacitors, or semiconductor chips. When 
applying multichip technologies, there may be hundreds of components attached to the 
printed circuit board. The ever-increasing level of complexity of printed circuit boards 
has forced primitive boards using a few yards of printed wiring in the 1960s, to evolve 
into sophisticated, multi-layered structures requiring kilometers of printed wiring. This 
increase in the circuit board level of sophistication can be attibuted to the integration 
of semiconductors and an increased need for I/O capabilities. [Ref. 2] 

The most complex printed circuit boards contain kilometers of copper intercon- 
nection, roughly 50 to 100 microns wide, and half as thick. These boards distribute 
KW|m? of power internally in very densely packed layers of copper. The drive for higher 
performance means that there is a greater requirement for power-handling and cooling 
capabilities. [Ref. 2] 

It is the responsibility of the designers to ensure that cooling on the printed circuit 
board is adequate under all possible load conditions in order to allow proper perform- 
ance of the individual components and the board as a whole. Therefore, it is imperative 
that the designers understand and are able to predict the temperature distribution on 
multilayered structures prior to prototype production. The overriding reasons for per- 
forming a precise thermal analysis are to increase component reliability, ensure proper 
material selection, reduce the possibility of catastrophic thermal failure, and guarantee 


electrical performance. 


A. RELIABILITY 

Reliability 1s defined as the probability that a component is functioning as designed, 
while failure is defined as the probability that a component is not functioning as de- 
signed. There is a predictable relationship between the operating temperature of elec- 
tronic components and reliability. The materials used in the fabrication of these 
components have thermal limitations, and should these thermal limitations be exceeded, 
the physical and chemical properties of the material are affected, and the device fails. 
[Ref. 1] 

For a large number of components, a typical plot of failure rate as a function of time 
is shown in Figure 2. Failures at short times are called early fails, or infant mortality, 
while failures at long times are called wearout fails because they result from usage. At 
all times, failures can occur from intrinsic mechanisms, or from random overstress. [Ref. 
3] 

Provided the device has been adequately designed, early failures can arise as a result 
of manufacturing defects. Defects that occur early on, or the “burn-in” period, are 
considered to be the result of poor or inadequate quality control mechanisms in the 
manufacturing process. [Ref. 2] 

Of greatest concern are the failures that occur during the useful life of the device 
because the probability of failure during this period should be nearly zero. Should a 
device fail during its useful life, the. probable cause of failure would be due to a variety 
of external factors, and are unpredictable. [Ref. 2] 

As time and usage progress, the terminal period of wearout is encountered, usually 
well past the the end of system life (EOL), in which the failure rate increases. Con- 
ductor electromigration is a typical example of a wearout mechanism, in which the 
electron flow itself causes irreversible mass flow, which causes the formation of voids and 
consequent conductor failure. Printed circuit boards seldom have significant failure rates 
when they are produced with sufficient characterization and control. When failures do 
occur, they can normally be attributed to either manufacturing defects or lack of integ- 
rity to the design specifications. The most notable concerns are the resin and laminate 
effects that can give rise to insulation integrity. As printed circuit boards become more 
densely packaged and continue to spread into more unconventional environments such 
as homes, automobiles or marine engines, there is a need to establish improved materials 
and process controls for even better reliability. Thermal analysis of the printed circuit 
board becomes one of the fundamental facets of the design process in order to better 


characterize the printed circuit board. [Ref. 2] 


The infant mortality,’or early fail 
portion of the curve is the result 

of failures during the "burn-in" 
period. The high failure rate shown 
on the right side of the curve is due 
to usage, or ''wear-out." 
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Figure 2. Failure rate as a function of time for typical component. 
Saprce: References: 


B. MATERIAL 

The fabrication of printed circuit boards results in the joining of different materials. 
The materials selected can have a significant impact in thermal properties of the printed 
circuit board. In the largest circuit boards hundreds of amps may be switched at once. 
As packaging densities increase thermal, mechanical, electrical, and chemical coupling 
becomes very strong. In view of the complexity of electrical structures today, computer 
modeling of total thermal response of a printed circuit board is a requirement conducive 


to understanding one of many interactions. [Ref. 2] 


C. THERMAL FATIGUE AND CATASTROPHIC THERMAL FAILURE 

Printed circuit boards are comprised of dissimilar materials that expand at different 
rates of heating. Table 1 shows the thermal coefficients of expansion of materials com- 
monly used in the fabrication of printed circuit boards. The differential expansion of 
mismatch must be accommodated by the various elements on the board. The increasing 
levels of packaging densities and board complexity dictate the need for designing a 
thermal environment that can accommodate the diverse components that are in close 


proximity to each other. [Ref. 4] 


Table 1. MATERIALS USED IN PCB PRODUCTION 


Coefficient of Thermal Expansion 10-7/°C 
C—O 
170 (ay) 600 (2) 


Source: Reference 2. 
















As previously mentioned it 1s very important that the board designers have an 
understanding of the operating environment in which the board will be operating in or- 
der to incorporate into the design the tolerances that will allow the product to operate 
reliably. 

Catastrophic thermal failure is defined as an immediate, thermally induced, total loss 
of electronic function in a specified component. This type of failure comes as a result 
of excessive temperature, or a thermal fracture. Catastrophic failure comes about as a 
result of many factors including the operating environment, equipment history, me- 
chanical loading, and operational modes of the component. Although it is difficult to 
predict the temperature at which thermal failure may occur, it is possible to establish 
with the aid of thermal analysis the boundaries at which the board can be expected to 
operate reliably and within its useful operating life cycle. [Ref. 1,2] 


Il. HEAT TRANSFER 


The degradation of the heat flow capabilities of a printed circuit board can lead to 
reliability problems due to excessive operating temperatures. It is imperative that de- 
signers incorporate into the board the capability to maintain temperatures within upper 
Operational limits while operating in all possible environments in which the board will 
be exposed. [Ref. 2] 

Heat transfer is defined as all energy flows that arise as a result of temperature dif- 
ferences. Because the components mounted on printed circuit boards and, indeed, the 
printed circuit boards themselves are not one hundred percent efficient, heat is gener- 
ated. The primary modes of heat transfer are conduction and convection. Conductive 
modes include mechanical thermal contact and solid thermal interfaces between materi- 
als, such as copper, solder, or epoxy layers. Convective modes include natural and air- 
forced cooling, and forced liquid cooling. Radiation is also a factor; however, it is not 
as significant as conduction and convection at the temperatures in which printed circuit 


boards operate. [Ref. 5] 


A. CONDUCTION 
Conduction is the transfer by molecular motion of heat between one part of a body 
to another part of the same body or one body and another in physical contact. [Ref. 1] 


For the case of conduction, the heat flow equation is the basis for understanding this 


behavior 

mar (1) 
where 

QO = heat flow along the thermal gradient, af 


A = the area through which the heat 1s flowing, m? 


k = the thermal conductivity of the material, — 


AL = change in temperature per unit length, & 





Rearranging Equation (1) leads to 
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which is the thermal resistance 
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Equation (3) demonstrates that thermal resistance will increase with an increase in 


path length for heat flow, AX, with a decrease in area of heat flow, A, or change in 


conductivity & to a lower value. 
1. General Equation of Heat Conduction 


The general equation of heat conduction is 
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where 
k 
p= density, — 
m 


! y 
C= specific heat, ke*C 
7 = temperature, °C 


x,y, and z= cartesian coordinates, m 


t=time, sec 


- W 
k = thermal conductivity, ——> 
oe mora 

. : 

q = internal heat generation, a 


(4) 


Assuming k, C and p are independent of temperature, direction, and time, the resulting 


equation 1s 
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where 
an k mr 
a = thermal diffusivity, TG’ SEC 


Several variations of the general conduction equation exist. Fourier’s equation, 


which contains no heat sources 1S 
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Another variation, known as Poisson’s equation, solves a system in which tem- 


perature 1S not time dependent 
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The last variant of the general equation of conduction 1s intended for « system 
Operating in a Steady-state condition, and does not contain any heat sources. This 


equation, known as La Place’s equation Is 
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2. Single Plane Slab Modeling 
For illustrative purposes, a single plane slab, with 7, and 7, being the face 


temperatures, 1s considered. If Equation (8) 1s limited to only one coordinate then 


aT 
a ( (9) 
dx’ 


If Equation (9) is integrated twice, and boundaries are established, then the temperature 


distribution across the slab can be expressed as follows 


1] 


T= T,-+(T, -T,) (10) 


If Equation (10) is substituted into Equation (1), then a solution for the heat flow across 
the single plane slab can be obtained 


seo Jy) 
¢— A 


kA 
pe =) (11) 
It is known that Ohm’s law relates the voltage drop across a resistor with the 
current flow through the resistor, V=/R. Ohm’s law is analogous to Equation (1), 
where the current flowing through a resistor is equivalent to heat flow, and the electrical 
resistance is equivalent to the thermal resistance expressed in Equation (2). 


Ohm’s law can be expressed as 
i= (12) 


where the analogy between Ohm’s law and Equation (1) is 


Current /<>Heat flow q 
Potential V<>Temperature difference AT 


Resistance Re>Thermal Resistance R 


It can be observed that for the heat flow across a simple plane slab described 


by Equation (11), the thermal resistance is described by Equations (2) and (3). 


B. CONVECTION 

Convection is defined as the process by which thermal energy is transferred to or 
from a solid by a fluid flowing past it. If the movement of fluid is a result of a temper- 
ature differential, then the process is called free or natural convection. When natural 
convection is present, the movement of fluid can be accelerated by increasing the tem- 
perature differential. When a pressure differential is introduced to force the movement 
of the fluid by using a pump or a fan, the process is called forced convection. [Ref. 
eS), 

Newton's law of cooling states that the heat flow through a body ts proportional to 
the normal area and the temperature difference between the body and the surrounding 


fluid [Ref. 2]. To make it an equality, a proportionality constant, A, is introduced. The 


proportionality constant is the surface heat transfer coefficient. Newton's law of cool- 


ing can be expressed as 
q = hA(Ty — Ty) (13) 
where 


his a proportionality factor known as the surface heat transfer coefficient 


Comparing Newton's law of cooling with Fourier’s law leads to an expression re- 
lating the surface heat transfer coefficient to thermal conductivity, the surface fluid 


temperature differences, and the wall temperature gradient of the fluid 


eg 
= Coram 

. ae (14) 

AAT AT 

Consequently, correlating heat transfer coefficients must be based on the depend- 
ence of A on the thermal conductivity of the fluid and on the ratio of the wall temper- 
ature gradient to the temperature difference. [Ref: 1] 

1. Electrothermal Analog 


For the case of convective heat transfer the thermal resistance is represented by 


i's 
R=7- (15) 


Consequently, the total thermal resistance is now defined for a single slab and convective 


heat transfer on its two faces as follows 
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The heat transfer equation can now be represented as 


(T, — Ty) 
i. =a Lj - 
antes, 


C. RADIATION 
The third mechanism for heat transfer is radiation in the form of electromagnetic 


waves. The rate at which a body radiates thermal energy 1s proportional to the area 
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of the body and to the fourth power of the absolute temperature [Ref. 5]. This result, 
found empirically by Josef Stefan in 1879, is written as 


q=ecAT (18) 
where 


g = power radiated, HV 
A = area, m? 


eé = emissivity of the body, a value between 0 and 1 


o = Stefan-Boltzmann constant, 5.6703 x lon: 
meK 

When radiation falls on an opaque body, part of the radiation is reflected and part 
is absorbed. Light-colored bodies reflect most of the radiation, whereas dark bodies 
absorb most of it. [Ref. 5] 

Materials employed in the manufacture of electronic components are classified as 
gray. Gray bodies are diffusely-reflecting opaque surfaces. These surfaces reflect equal 
amounts of energy over the thermal radiation spectrum in all directions. [Ref. 2] 

1. Transformation of the General Radiation Equation 

The use of the thermal radiation equation in analytical studies 1s a difficult task 
due to the fourth power relationship with temperatures. The complexity of the calcu- 
lations involved with radiation dictate the need for the aid of a computer [Ref. 2]. The 


general equation for radiation interchange 1s 
q = oF ,F,A(T; — T,) (19) 
where, as before, o 1s the Stefan-Boltzmann constant and 


F,= shape factor accounting for source and receiver arrangement. 
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F,= emissivity factor accounting for properties of the source and receiver. 


7, = temperature of the source, °K 


T, = temperature of the receiver, °K 


It must be observed that for the radiation case absolute temperature is the norm. 
Radiation based heat transfer is represented by transforming the general radi- 

ation equation, Equation (19), into a form compatible with Fourier’s law. Linearization 

of the general radiation equation is the method used to produce the desired result. The 


difference between the two fourth powers can be reduced to: 

= 7, )=(1; + TET, ) 

=(T; + T;)(T, + T,(T,- T,) (20) 
Inserting this into Equation (19) results in 

q = oF ,F_A(T, + T-)(T, + T,\(T; — T,) (21) 


A radiative heat transfer coefficient may be defined as 


h, = oF ,F{T; + T-)(T, + T,) (22a) 
or 
oF ,F(T. + TT. + T.T + T°) (226) 


Substitution of h, into Equation (21) shows that radiative heat transfer can now be 
treated similarly to convection at the boundary [Ref. 6]. Thermal resistance in the case 


of radiation heat transfer can now be denoted as 


Re am (23) 
Now that all three methods of heat flow have been discussed, a graph depicting 
the electrothermal equivalent is shown in Figure 3. Here the plane slab has both radi- 
ation and convection on its two faces, and the nonlinearity of the radiation coefficient, 
h, , leads to a more detailed analysis procedure. 
When considering the calculation of heat transfer by radiation, it is usually 
necessary to approximate real material behavior with gray-body idealization. There is 


also a great deal of data required in order to model radiation for a real body not only 








T 


Figure 3. Radiation, convection, and conduction electrothermal equivalent. 


Source: Reference 8. 


because the surfaces are non-black, but also because the non-uniformity of the body’s 
surfaces creates greater complexities in the collection and processing of data. [Ref. 7] 
Basic problems, such as the one shown in Figure 3, defy a simple solution be- 
cause of h, Therefore, a computer-aided approach becomes necessary in solving even 
elementary problems. Models exist for aiding the user in solving radiative heat transfer 


problems; however, their development is not within the boundaries of this thesis. 
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IV. FINITE DIFFERENCE ANALYSIS 


Finite difference methods represent the most appropriate approach in determining 
the temperature distribution within an electronic coinponent. By using a finite differ- 
ence approach, complex problems involving difficult geometries, non linearities, or 
complicated boundary conditions can be successfully approached. The purpose of this 
section is to explain the fundamental concepts behind finite difference analysis and its 


use in solving differential equations. [Ref. 6] 


A. FUNDAMENTALS OF FINITE DIFFERENCE ANALYSIS 
The general equation for heat transfer by conduction, Equation (4), can be ex- 
pressed in a three-dimensional variant as 
é Be 0 OT 
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where all the equation components have been previously defined. In order to formulate 
a problem using Equation (24), it is first necessary to approximate the first and second 
derivatives. 
I. First and Second Derivative Approximation 

The derivative of a function at a point can be formulated by a finite difference 
approximation. It is first necessary to determine the slope of a line tangent to the point 
of interest on the plot of temperature as a function of x, y, or z. (Figure 4). In choosing 
x as the coordinate system, it is observed that the slope of the line (7m) 1s defined as 
= I(xo + Ax) — flxo) 


li 


Ax—0 Ax oe 


The term of Equation (25), Lx, + Ax) —/(%)J/Ax is referred to as the difference quo- 
tient, and is the ratio of the change in the value of the function at x, and x, + Ax to the 
change in x. The limit of the difference quotient is called the derivative of the function 
dle INel). 6) 


f(x) = Jim, < (26) 





T=S(x) we 


Figure 4. First and second derivative approximation. 


Source: Reference 8. 
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By letting f(x) = 7(x), the derivative for this particular curve at the point x, can be de- 


rived 


F (4) = im Tet Ba) = Te) an 


As Ax > 0 Equation (27) becomes 


_ aT _— Txt Ax) — Te) 


a8 (Xq) dx Ax (28) 
Thus, for Ax finite, but very small 

; ol 
I°(%0) = A (29) 


The second derivative, using the definition of a derivative, will be the ratio of the 
change in value of the first derivative to the change in Ax [Ref. 8]. Again, looking at 


Figure 4, the second derivative is 


Jf" (Xo + Ax) — f'(%) 
ey ce ee 


if (Ce), = sien ne (30) 
P 2T (xy + 2Ax) — 2T (xy + Ax) + T(x) 
& () = (31) 


For both Equations (28) and (31), the approximations are valid as long as Ax is small. 
Taylor’s theorem states that a function can be approximated by a polynomial 


of the form 


SG 
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fix) = fla) + f'(a)(x — a) +9 (e- a) ++ 
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where the polynomial is for the function expanded about x =a [Ref. 8]. The Taylor se- 
ries for the function 7(x) at T(x) + 2Ax) is 


dT(x, + Ax a’ T(x, + Ax 
T(x9 + 2Ax) = Theiy + x) ¢ OEE) gg LEEDS at (33) 
Xx ax 
a: Tien ad" Tx, + Ax 
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ane for the function /(x) at 7(x,): 


T(Xo) = T(xp + Ax) - 


dT (x, + Ax d? T(x, + Ax 
(xp Lact (xp Re. (34) 


ax ake 
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If Ax is kept small, then the terms above second order become negligible. When 
Equations (33) and (34) are added, an expression for the second derivative can be ob- 
tained 


d Tixy + Ax) 5 
T(x + 2Ax) + T(x) = 2T(x) + Ax) + — a Ax 
x 


(35) 
d°T(xy + Ax) — T(xp + 2Ax) — 27 (xp + Ax) + Tx) | 


* —3 f"(20) (36) 


As previously stated, the configuration to be analyzed can be subdivided into 
small, finite subvolumes considered to be isothermal. The centroid of each subvolume 
is called a node, and the node is representative of the entire subvolume. Nodes may are 
connected to their adjacent nodes through thermal resistances and the nodal analysis for 
the node temperatures can best be solved by means of a computer aided model [Ref. 8]. 

The first law of thermodynamics states that energy can neither be created nor 
destroyed, but can be transformed from one form to another. An energy balance may, 
therefore, be formed on the typical node shown in Figure 5. For the sake of simplicity, 
the environment is composed by a single node, node 5 (Figure 6). Close observation of 
node 5 shows that it is connected to nodes 2, 4, 6, 8, 14, and 100. If any energy 1s di- 
rectly applied to or removed from this node, it would become part of the equation. The 


node equation for node 5, with rate of heat input g, becomes 
92 + 9a + 96+ 98 + 914 + F100 - % = 9 (37) 


where each one of the g values with a numerical subscript represents the rate of heat 


flow from node 5 to the node indicated by the subscript. [Ref. 1,6,8] 
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Figure 5. Graphical representation of a single node. 


Source: Reference 6. 





Equation (37) does not display the node temperatures; however, g = A7T/R and 
Equation (37) can be combined to produce 




















a= Af aT =“ (15-1) (38a) 

g=AAar=— (7,- 1) (386) 

=A r= (1-1) (38c) 

os = 4a ay a Cis 3) (384) 

nen *Aar=— (7, - To (38e) 

9qio0 = Aa AT= St (Ts — T\00) (38/) 
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Where 7}, 18 on the top face of node 5 and, Ax, Ay, Az = the distance between nodes in 
the x,y, and z directions respectively. Assuming Ax = Ay=Az=1, then the node 


equation becomes 
A(T; — Ty) + k(T; — Ts) + A(T; — Ts) + (39) 
k(Ts — Tg) + K(Ts — Ty4) + 2K(Ts — Tyo) = 4 


Or 
— Ty — Ty +175 — Ts — Te — Tr = + 2T hop (40) 


where the terms on the right hand side of Equations (39) and (40) are known values. 
{Ref. 1,6,8] 

A node equation is written for each node in the configuration, and in Figure 
8, Which contains 18 nodes, 18 node equations would be required to obtain the 18 un- 
known node temperatures. A larger problem having several hundred nodes would re- 


quire prohibitive amounts of time to solve or even enter data for computer-aided thermal 
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analysis. That 1s why a computer-aided model builder capable of constructing the 
needed input data for thermal analysis, 1s so desirable. The thrust of this thesis is not 
to solve the equations, but to allow the user to model a printed circuit board to his or 
her specifications with minimal effort, and then produce a data file compatible for use 


in the thermal analyzer. 
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V. THE MODEL BUILDER 


As previously indicated, the model builder currently in use in the thermal analysis 
software package requires time-consuming manual data entry. This thesis proposes that 
a model builder, PCB, can be incorporated into the software package. PCB is a menu- 
driven, user-friendly program which will assist the user in developing a thermal model 
of a printed circuit board with pre-defined geometries. PCB will generate a properly 
formatted ASCII output data file for use with the thermal analyzer. The current version 
of the thermal analyzer is designed to handle a maximum of 300 nodes [Ref. 9]. An 
updated version capable of handling the amount of nodes generated by the model builder 
is under development; however, abridged versions of the model builder have been suc- 
cessfully tested with the current version of the model builder. This chapter addresses 
many aspects of the model builder, including a terminal session demonstrating the design 
process used in modeling a circuit board, and the interface between the output data file 
and the thermal analyzer. Also addressed is the methodology used in validating the PCB 
model builder output data file. Finally, limitations and possible improvements to the 


model builder are discussed. 


A. FEATURES 

PCB is a menu driven interactive program providing the user with a variety of sys- 
tematic choices progressing in a logical sequence that eventually will lead to the devel- 
opment of the model. Upon entering the program, the user will be provided with the 
option of viewing an optional overview. First-time users should consider reading the 
overview. The program is designed to operate with upper-case letters, therefore, it 1s 
suggested that the caps-lock key is pressed before starting the program. 

Following the overview option the user will be asked to choose either SI or English 
units. For this choice as well all others, the program provides the user with the op- 
portunity to correct any erroneous entries. 

The third section queries the user for the number of copper layers desired in the 
design. The allowable number ranges from one to four layers. For every choice of 
number of copper layers, there are four possible board configurations the user can select. 
If one copper layer is selected, the program will display the following choices of aspect 


ratios: 


1. 10 BY 36 PROVIDES A 1:3.6 RATIO WITH 720 NODES. 
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2. 1OBY 35 PROVIDES A 1:3.5 RATIO WITH 700 NODES. 
3. 12 BY 25 PROVIDES A 1.2:2.5 RATIO WITH 600 NODES. 
4. 15 BY 20 PROVIDES A 1.5:2.0 RATIO WITH 600 NODES. 


When two copper layers are selected, the user will be provided with the following 


choices: 
1. 10 BY 18 PROVIDES A 1:1.8 RATIO WITH 720 NODES. 
2. 9 BY 20 PROVIDES A 1:2.11 RATIO WITH 720 NODES. 
3. 8 BY 20 PROVIDES A 2:5 RATIO WITH 640 NODES. 
4. 10 BY 15 PROVIDES A 1:1.5 RATIO WITH 600 NODES. 


If three copper layers are selected, the user will have the following choices: 
1. 5 BY 24 PROVIDES A 1:4.8 RATIO WITH 720 NODES. 
2. 10 BY 10 PROVIDES A 1:1 RATIO WITH 600 NODES. 
3. 12 BY 10 PROVIDES A 1.2:1 RATIO WITH 720 NODES. 
4. 8 BY 15 PROVIDES A 1:1.875 RATIO WITH 720 NODES. 


By selecting four copper layers, the program will provide the user the following 
options: 
1. 5 BY 18 PROVIDES A 1:3.6 RATIO WITH 720 NODES. 
2. 9 BY 10 PROVIDES A 1:1.11 RATIO WITH 720 NODES. 
3. 8 BY 10 PROVIDES A 1:1.25 RATIO WITH 640 NODES. 
4. 13 BY 6 PROVIDES A 1:1 RATIO WITH 624 NODES. 


After selecting the desired aspect ratio, the user must specify the epoxy character- 
istics. The program will prompt the user to enter the epoxy layers length, width, 
thickness and thermal conductivity. An option available to the user will be the choice 
of a default epoxy thickness value, 0.0625 in., (0.15875 cm.). 

Once the epoxy characteristics have been entered, the program will then display the 
corresponding length and thickness for the copper layers selected. The length and width 
for the copper layers will be the same as those of the epoxy layers. The user will then 
be provided with the option of entering the thickness of each selected copper layer either 
by length or by weight. If input by length is selected, the units should be consistent with 
the system of measurement selected; however, if entry of thickness by weight 1s selected, 
then the entries should be in ounces. A copper layer weighing one oz. is equivalent to 
a copper layer 0.0014 in. thick with a surface area of 144 in? . The program adjusts 
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thelaver thickness to the surface area provided automatically. Entry of thickness by 
weight is a standard practice used in industry. [Ref. 2] 

After entering the copper layer thickness the program will query the user for the 
thermal conductivity of copper. After this entry, the program will redisplay all infor- 
mation pertaining to the copper lavers and will require the user to confirm all entries in 
Orden (Gene ceced: 

The following section requests the input of initial and ambient temperatures. PCB 
has six ambient temperatures. 

After the temperatures have been entered the program requires the user to input the 
external heat sources. PCB provides for heat input into the upper copper layer. There 
are four methods of external heat input from which to choose. The first alternative al- 
lows for a total rate of heat applied to the upper surface. An entry for this choice would 
be divided by the number of nodes and distributed uniformly. The second choice pro- 
vides for the entry of average heat per unit area. The third alternative gives the user the 
ability to enter heat in specifically designated nodes. The last option provides for no 
heat input and was developed to test compatibility with the thermal analyzer. 

The following selection is the percent copper coverage for each copper layer. The 
program will prompt the user to enter the percentage of copper coverage for each layer 
in the printed circuit board. 

Finally, the program requires a name and a title for the file that will be created. This 


file will be in the proper format for use with the thermal analyzer. 


B. THE THERMAL ANALYZER INPUT DATA FILE 

The model builder generates an ASCII data file from the physical characteristics of 
the printed circuit board provided by the user. In order for the input data file to be ac- 
ceptable to the thermal analyzer, it must be in a specific format which is compatible to 
the thermal analyzer [Ref. 9]. Because each value and position of the output data file 
has a meaning to the thermal analyzer, and is not readily identifiable to the user, it 1s 
beneficial to describe each line and data set of the output data file and their relationship 
to PCB. Figure 7 shows a partial output data file. 

Line one is the title line. It may be left blank or may contain up to 79 alphanumeric 
characters. The user-selected title appears at the top of the data file. [Ref. 6] 

Line two is the problem data line. Jt has nine entries of which two are under user 
control, the number of nodes under consideration and the unit type. One entry, the 


number of constant temperatures, is preset at six for this specific model. The remaining 
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78.000 78.000 
? 73381 2) 
1.860 2626 
? 11 31 
0626 2626 
? 21 41 
2624 0626 
? 31 $1 
624 2626 
? 4) 6} 
26246 26246 
? $1 71 
2 626 626 
? 6) 81 
626 2626 
7 71 7341 
2626 1.560 
? 8) 101 
624 2624 
? 91 331 
2626 2 626 
? 10) 121 
6246 2624 
? 331 131 
2626 0626 
? 321 14) 
0626 2626 
? 131 1$1 
2626 2626 
? 16) 16) 
0624 2626 
? 1$1 171 
2 6246 2626 
? 16) 181 
0626 ° 6246 
? 171 191 
2624 626 
? 181 20) 
0626 0 626 
? 191 21) 
2626 2 624 
7 20) 221 
2 626 2 624 
? 233 233 
0626 2626 
? 221 24) 
2626 2626 
? 231 251 
2626 0626 
Figure 7. 
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69.000 


7$21 

0276 
7$21 

0276 
7$2) 

2276 
7321 

0276 
7$21 

e276 
7$21 

0276 
7$21 

0276 
7$21 

0276 
11 

0337 
21 

2137 
31 

0337 
6) 

0337 
$} 

0137 
61 

0337 
71 

0337 
91 

0137 
91 

0137 
101 

0137 
113 

0137 
121 

2337 
131 

2137 
341 

013? 
151 

013? 
16) 

0137 


6 6 | 
é ) ) 
8000000 €9.00000 
78.000 89.000 09.000 
91 7$1) @11 9991 
0187 2.787 2.787 0363 
30) $1 $21 9991 
0137 2.787 2.787 oS6S 
111 7$11 831 999) 
0337 2.787 2.787 863 
121 7811 06) 999) 
0137 2.787 2.787 363 
131 7$11 @$1 999) 
0137 2.78? 2.787 363 
161 7311 61 9991 
0137 2.787 2.787 363 
15) 7$11 @7) 9991 
0337 2.787 2.787 863 
161 7311 63) 999) 
0337 2.787 2:78? 363 
171 7813 @91 9991 
0337 2.787 2.787 -363 
181 7311 901 9991 
0337 2.787 2.787 o36S 
19) 7$11 911 9991 
2137 2.787 2.787 -863 
201 7$11 921 9991 
0137 2.787 2.787 S63 
211 7311 93) 9991 
o137 2.787 2.787 363 
22) 7$11 961 9991 
0137 2.78? 2.78? 863 
23) 7831 951 9991 
0137 2.787 2.78? 863 
261 7511 961 999) 
0137 2.787 2.787 $63 
25) 7$11 971 9991 
0137 2.787 2.787 -S63 
26) 7$11 981 9991 
0137 2.787 2.787 -363 
271 7$11 991 9991 
0137 2.787 2.787 S63 
28) 7311 100) 999) 
2337 2.78? 2.78? 863 
291 73811 1033 9991 
0137 2.787 2.787 363 
301 7311 1021 9991 
0337 2.787 2.787 -S63 
311 7311 103) 999) 
0337 2.787 2.78? 863 
32) 7313 1863 9991 
2337 2.787 2.787 
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Sample PCB partial output data file. 


entries have applications to models associated with heaters, unique exponents, secondary 
heat input, temperature coefficients and curves, and nodes controlling fast heat. These 
entries are not applicable to this model and are preset to Zero. [Ref. 6] 

Line three places a zero at three points and is beyond the user’s control. Therefore, 
no further discussion is required. [Ref. 6] 

Line four is the problem capability line. This line defines the maximum values for 
the entries in line two. The first entry 1s 750, which is the number of nodes for which 
the analysis is dimensioned. The number 750 is significant because the first constant 
temperature will be assigned to node number 751. The second entry is 50 which re- 
presents the largest possible number of constant temperatures in accordance with the 
analyzer dimension statement. The third entry is set to 6 and does not change. This 
entry is related to heaters and 1s not applicable to the model. The balance of the entries 
in line four represent a listing of data sets that are required for the particular analysis 
at hand. PCB uses three data sets that will be discussed in what follows. [Ref. 6] 

The fifth line contains five values that relate to the accuracy level that the thermal 
analyzer will achieve. These entries are preset. The first value provides the level of ac- 
curacy between iterations. The accuracy level number is critical because too small a 
tolerance will cause the computer to run for excessive amounts of time, and too large a 
number will provide inaccurate results. The second value is the damping factor used 
between iterations in order to prevent temperature oscillations between iterations. The 
third number provides the maximum number of iterations. If erroneous data 1s entered, 
the computer will not run for excessive amounts of trme. The fourth value is the con- 
vergence factor which adjusts the damping in order to close the critical value. The fifth 
entry is the initial temperature at which the iterative process begins. This value is sup- 
plied by the user. [Ref. 6] 

Line six contains the temperature dependent coefficients, and is not used in this 
model. Line seven contains up to 50 constant temperature inputs. This model has six 
ambient temperatures entered by the user. [Ref. 6] 

The following lines contain all pertinent information concerning the n-node 
equations. Each node requires two lines of data. Even-numbered lines are used for 
specifying the nodes that interact with the node in question and the modes by which this 


interaction takes place. For example, a line of the form: 


6 7551 21 752° 7 Sri 
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is for node number one of aspect ratio selection one. The first entry indicates the num- 
ber of connections to that node. The second entry is indicating that node 755, an am- 
bient temperature, is connected to the node in question and the | indicates that the 
connection is conductive. The same procedure applies to the rest of the values on the 
line. If the entry is 9991 then the heat input 1s external. (Ref. 6] 

Other data sets that appear relate to unique exponents, secondary heat, and 
temperature-dependent heat input curves. These data sets are not used in this model. 
[Ref. 6] 


C. PCB MODEL BUILDER SAMPLE PROBLEM 
The following terminal session is a typical example of the printed circuit board 
model builder, PCB, and its capabilities. Figure 8 depicts the printed circuit board 
modeled in the terminal session. 
1. Printed Circuit Board Specifications 
. Unit system: British. 
. Copper layers in printed circuit board: three. 
. Aspect ratio desired: 12 by 10, with 720 nodes. 
meepoxy laver length: 7.2 i. 
. Epoxy layer width: 4.2 in, 
. Epoxy layer thickness to default value: Yes (.0625 in.) 
. Epoxy layer thermal conductivity: 0.087 Btu/hr/°F. 
Specify copper layer thickness by weight or length: length. 


Oo On HD OO HR WN = 


Thickness for the three copper layers, respectively: .20, .10, .11 in. 


pe 
© 


. Copper thermal conductivity: 243.000 Btu/hr/°F. 


— 
—" 


. Initial board temperature: 87 °F. 


— 
ho 


. Upper surface ambient temperature: 86 °F. 


—" 
WJ 


. Lower surface ambient temperature: 85 °F. 


pa 
SS 


. Right surface ambient temperature: 84 °F. 


15. Left surface ambient temperature: 86 °F. 


—" 
N 


. Front surface ambient temperature: 86 °F. 

17. Rear surface ambient temperature: 86 °F. 

18. Heat input to the upper copper layer: 300.23 Btu/hr/°F. 

19. Percent copper coverage for copper layers, respectively: 87.3, 76.2, 65.7 %. 
20. Name of output data file: LAYERS. 
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21. Header of the output data file: THIS IS THE OUTPUT OF THE THERMAL 
ANALYZER. 
2. Terminal Session 
Using the board specifications and Figure 8, the terminal session for the PCB 
program ts as follows: 


THIS PROGRAM WAS WRITTEN TO INTEGRATE WITH EXISTING 
THERMAL ANALYSIS SOFTWARE AND TO REDUCE THE AMOUNT 
OF TIME REQUIRED FOR DATA ENTRY. 


WOULD YOU LIKE AND OVERVIEW OF THE PROGRAM PRIOR TO 
BEGINNING? ENTER Y FOR YES AND N FOR NO: N 


KREKEEKEKEEREKREEEKEEKEKEEEEEEEKKEKEKEKKKKKKRKKKRKKKKREKSE 
KEEKEKEKEEKEKEKEEEEEEEKEKEEEEEKEEEEKEEKEEKKKKEKKKKKKKKESE 


zk & xx 
** PLEASE SELECT UPPER CASE LETTERS “x 
belied PRIOR TO BEGINNING ax 
kk xx 


KEKKEKEKEEKEEEKEEEEEEEEEEKEEEEEEEEKKKKEKKKEKKKKKRKEKEE 
KRAEKKEKKKEKEKEEKEKEEEEKEKEEKEEKEEKKKKKKKKEKKKKKKKKKKKKSE 


PRIOR TO ENTERING DATA INTO THIS PROGRAM ENSURE THAT 
YOU HAVE A DRAWING OF YOUR DESIGN AND ALL PERTINENT 
DATA. 


PRESS <ENTER> TO CONTINUE 
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epoxy k=0.087 Btu/br/*F 





copper K=2493 Btu/hr/*F 


epoxy thick= 
©.@62e5 in. 


copper jayers vy VEROXxY ayers 


Figure 8. Printed circuit board modeled in the terminal session. 
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THIS PROGRAM IS CAPABLE OF OPERATIONS IN EITHER SI OR 
ENGLISH UNITS. AFTER THE SELECTION OF THE UNITS, ALL 
ENTRIES MUST BE COMPATIBLE. PLEASE MAKE YOUR SELECTION. 


S FOR SI NOTATION 
FE FOR ENGLISH NOTATION: E 
YOU HAVE SELECTED ENGLISH NOTATION. 


IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FOR NO: 


THE PCB CAN HAVE UP TO FOUR COPPER LAYERS 


HOW MANY COPPER LAYERS DO YOU DESIRE? 


PICK A NUMBER 1 TO 4: 3 


YOU SELECTED 3 COPPER LAYER(S) FOR THE PCB 


IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FOR NO: Y 


4 


YOU SELECTED 3 COPPER LAYERS, GIVING YOU THE FOLLOWING ALTERNATIVES: 


1. 5 BY 24 PROVIDES A 1:4.8 RATIO WITH 720 NODES. 
2. 10 BY 10 PROVIDES A 1:1 RATIO WITH 600 NODES. 
3. 12 BY 10 PROVIDES A 1.2:1 RATIO WITH 720 NODES. 
4. 8 BY 15 PROVIDES A 1:1.875 RATIO WITH 720 NODES. 


PLEASE SELECT A NUMBER 1 THROUGH 4: 3 


YOU SELECTED NUMBER 3 OF THE FOLLOWING 4 ALTERNATIVES: 


1. 5 BY 24 PROVIDES A 1:4.8 RATIO WITH 720 NODES. 
2. 10 BY 10 PROVIDES A 1:1 RATIO WITH 600 NODES. 

3. 12 BY 10 PROVIDES A 1.2:1 RATIO WITH 720 NODES. 
4. 8 BY 15 PROVIDES A 1:1.875 RATIO WITH 720 NODES. 


IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FOR NO: re 
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KH KK KK KH KKK KK IKK KKK KKK KKK KKK KKK KEKE KEKE KEK KKK EK KEK KEKEKEKEKK KKK KKK 
KKKKKKKEKKKKEKKKKEEEEKEKEEKEKEPOXY CHARACTERISTICS & kK EK EKKKKKEKKKRKKEKKEKEK 
CETTE CTICTCCTOCCTOCCTCCCTCCTOCTCCTCCCTCCLCCCCLCCCCTCCCCCCC CCC TTT TTT TTT TTT TTT tt 


ALL ENTRIES ARE IN ENGLISH NOTATION 


ENTER EPOXY LAYER LENGTH (in): 7.2 

ENTER EPOXY LAYER WIDTH (in): 4.2 

THE EPOXY LAYER THICKNESS WILL DEFAULT TO 0.0625 in (0.15875 cm) 
DO YOU WANT TO CHANGE THE EPOXY THICKNESS? (Y OR N): N 


ENTER EPOXY LAYER THERMAL CONDUCTIVITY (Btu/hr/F ): 0.087 


YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE EPOXY LAYER: 


1.) LENGTH: 7.2000 in 
2.) WIDTH: 4.2000 in 
3.) . THICKNESS. -0625 in 
4.) k: -0870 Btu/hr/F 


DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR YES AND N FOR NO: N 
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KKEKKKKKKEKKEKEEKEKEKEKEKKKKEKKKEEKEKEKEKEKKKEKEEKEKKEEEKKEEKEKEKKKKKKKEKKKKEKEEKKEE 


kkkkkkkkkkkkkkkkKCOPPER LAYER CHARACTERISTICS EEKEKKKKKKKKKKKKKKKKKKKKKKE 
KKK KKK KKKKEKKKEKKEKKEKKEKKEEKEEEKEEKEEEKEEKKKKEEEEEKKKKKKKK KKK EKKKEEK 


ALL ENTRIES ARE IN ENGLISH NOTATION 


[THE COPPER LAYER LENGTH IS THE SAME AS THE EPOXY LAYER: 7.2000 in 
[THE COPPER LAYER WIDTH IS THE SAME AS THE EPOXY LAYER: 4.2000 in 


{OU SELECTED 3 COPPER LAYER(S) FOR THE PCB. 
OU NOW WILL BE ASKED TO ENTER THE COPPER LAYER THICKNESS 


90 YOU WANT TO SPECIFY THICKNESS BY LENGTH OR WEIGHT? 
SNTER L IF YOU WANT TO ENTER LENGTH, W IF YOU WANT TO ENTER WEIGHT: L 


SNTER THE THICKNESS FOR LAYER 1 (in): .2 


<NTER THE THICKNESS FOR LAYER 2 (in): .1 


SNTER THE THICKNESS FOR LAYER 3 (in): .11 





*NTER COPPER LAYER THERMAL CONDUCTIVITY (Btu/hr/F ): 243 


HAVE MADE THE FOLLOWING ENTRIES FOR THE COPPER LAYER(S). 


LENGTH: 7.2000 in 
WIDTH: 4.2000 in 
THICKNESS LAYER 1: ~2000 in 
THICKNESS LAYER 2: ~1000 in 
THICKNESS LAYER 3: ~-1100 in 
kK: 243.0000 Btu/hr/F 
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KEKKKKKEKEKKEKKKEKEKKKEKKKKEKKEEKEKKEKEKEKKKEKKKKKEKKEKKEKKKEKEKKEKKEKKKKKKKKKKKKKKKKKKK: 


kkkkkkkkkkkkkkkkkKKAMBIENT TEMPERATURE INPUTS&&&kKkkKKKKKKKKKKKKKKKKKKKKKS 
KKKKKKKKKKKKKKKKKKREKEKKKKEK EEE KEEKEEK KEKE EKEKKKKEKEKEKKEEREKEEKEEEKEKKKKKKKKKEK 


ENTER THE INITIAL BOARD TEMPERATURE (F): 87 


ENTER THE UPPER SURFACE AMBIENT TEMPERATURE (F): 86 


ENTER THE LOWER SURFACE AMBIENT TEMPERATURE (F): 85 


ENTER THE RIGHT SURFACE AMBIENT TEMPERATURE (F): 84 


ENTER THE LEFT SURFACE AMBIENT TEMPERATURE (F): 86 


ENTER THE FRONT SURFACE AMBIENT TEMPERATURE (F): 86 


ENTER THE REAR SURFACE AMBIENT TEMPERATURE (F): 86 


YOU HAVE MADE THE FOLLOWING AMBIENT TEMPERATURE ENTRIES: 


1. INITIAL BOARD TEMPERATURE 
2. UPPER AMBIENT TEMPERATURE: 
3. LOWER AMBIENT TEMPERATURE 
4. RIGHT AMBIENT TEMPERATURE 
5. LEFT AMBIENT TEMPERATURE: 
6. FRONT AMBIENT TEMPERATURE: 
7. REAR AMBIENT TEMPERATURE: 


DO YOU WISH TO MAKE ANY CHANGES? 


DO YOU WISH TO MAKE ANY CHANGES? 


87.000 
86.000 
85.000 
84.000 
86.000 
86.000 
86.000 


tes ba] haf bof bay Py Ory 


SELECT Y FOR YES AND N FOR NO: N 


SELECT Y FOR YES AND N FOR NO: N 
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a i i a a ee ee o_o anaana nme eee RMR RRR RRR RRR RAR RAR RR RR RRR RRRARZKZTETKTIEZKIK—IEKEIXKTr—XKEKK RX 
RAEKKKKAEKKKKKKKKKKKKKEKKEKKEKKXKKHEAT INPUTS &XKKKKKKKKKKEKKKKKKKKKKKKKEKKKKKS 
HKKKKKKKKKEKKEKKKKKEKEKKEKKKEKEKEEEKKEKKEKKEKKEKEKKEKKKKKKKKKKKKKKEKKKKKKKKKK 


HEAT INPUT TO THE PCB OCCURS ONLY ON THE UPPER COPPER 
LAYER. HEAT INPUT IS ACCOMPLISHED BY ONE OF THE 
FOLLOWING METHODS: 

1. ENTER AS TOTAL HEAT APPLIED TO THE PCB 

2. ENTER AS AVERAGE HEAT PER UNIT AREA 

3. ENTER HEAT NODE BY NODE 

4. NO HEAT INPUT 

PLEASE SELECT NUMBER 1 THROUGH 4: i 


YOU HAVE SELECTED NUMER 1 OF FOUR ALTERNATIVES. 
IS THIS THE DESIRED SELECTION? (Y OR N): Y 


YOU HAVE SELECTED TO INPUT HEAT AS A TOTAL HEAT 
APPLIED TO THE SURFACE. 


ENTER TOTAL HEAT APPLIED TO THE SURFACE (Btu/hr): 300.23 
IS THIS THE CORRECT ENTRY? (Y OR N): Y 


TOTAL HEAT PER NODE IS: 2.5019 Btu/hr 


<PRESS ENTER TO CONTINUE> 


ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR LAYER 1: 87.3 
ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR LAYER 2: 76.2 


ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR LAYER 3: 65.7 
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YOU HAVE SELECTED: 87.30% COVERAGE, LAYER 1 
76.20% COVERAGE, LAYER 2 
65.70% COVERAGE, LAYER 3. 


IS THIS YOUR DESIRED ENTRY? (Y OR N): y 


YOU SELECTED LAYERS FOR YOUR DATA FILE NAME 


IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FOR NO: 


THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY INTO THE 
EXISTING THERMAL ANALYZER, FURTHERMORE, THIS PROGRAM DOES 
NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE 
THE USER WILL NAME THE DATA FILE FOR EACH RUN OF THIS 
PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS, AND 
SHOULD NOT HAVE ANY SPACES. 


PLEASE ENTER THE DESIRED DATA FILE NAME: LAYERS 


YOU SELECTED LAYERS FOR YOUR DATA FILE NAME 


IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FOR NO: 


ENTER THE DESIRED TITLE TO BE PLACED ON LINE 
NUMBER ONE OF THE OUTPUT DATA FILE: 


THIS IS THE OUTPUT OF THE THERMAL ANALYZER 


DO YOU WISH TO CHANGE THE TITLE OF YOUR OUTPUT DATA FILE? 
ENTER Y FOR YES AND N FOR NO: N 


THE OUTPUT DATA HAS BEEN PLACED IN A FILE NAMED LAYERS 


<PRESS ENTER TO CONTINUE> 
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3. The output data file 
Figure 9 shows a partial listing of the model builder output computed in the 
terminal session. A detailed description of the output data file is provided in section B 
of this chapter. 
4. Thermal analyzer output 
The file obtained in the model builder produces an output for 720 nodes. In 
order to use this data file with the thermal analyzer the data file must first be shortened. 
In this particular case, the output file, LAYERS, shortened to eight nodes was entered 
to the thermal analysis program, TASS. Three files are produced by the thermal analysis 
program, TASS. The first two output files consist of detailed and summarized thermal 
analysis results, respectively. The third file is an error message summary in the event 
there is an error in the input data file. Figure 10 shows the detailed output file of the 


eight node thermal analysis. 


D. PCB MODEL BUILDER VALIDATION 

The output of the model builder was validated by following a three-step process. 
The first step in the validation process concerned format. In order to proceed with any 
other validation steps it was first necessary to ensure that all data fields were in accord- 
ance with the specifications set forth in the user’s manual of the thermal analyzer [Ref. 
9]. The thermal analyzer input data file consists of five lines and as many as seven input 
data sets. The thermal analyzer’s manual model builder program, THANSS, was in- 
strumental in verifying the accuracy of the model builder output data file format. 

The second step entailed content. The output of the model builder was checked for 
accuracy and completeness. Thermal resistances in the output data files were checked 
and compared with manually calculated benchmark models. A benchmark model was 
developed for each copper layer configuration. Other information such as node numbers 
and node relationships were also compared with the benchmark model and verified for 
accuracy. To aid in the verification process, the manual model builder program, 
THANSS, part of the thermal analyzer software was used to verify node relationships 
(Ref. 9]. Other items checked for correctness were specification codes relating to units 
(British or SI), temperature inputs, and thermal analyzer default values. 

The third and final step in the validation process required running the thermal ana- 
lyzer program with the PCB output data file as the program input. Abridged versions 
of the model builder were successfully executed with the current version of the thermal 


analyzer (maximum node capability is 300 nodes). The results of the thermal analysis 
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THIS IS THE OUTPUT OF THE 


THERMAL ANALYZER 


720 6 0 0 0 0 1 
0 0 0 
750 50 € 2 4 6 t's) 
0 
-0500000 .6666670 12 .8000000 87.00000 
86.000 686.000 86.000 84.000 8G.000 85.000 
7 7551 21 7521 1321 7511 1212 
242.445 87.280 41.249 20.625 -7 Ol -7021 2.50 
7 oe 3 7521 141 75121 1221 
87.280 87.280 41.249 20.625 -702 -702 2.50 
7 21 al 7521 152 7511 1231 
87.280 ' 87.280 41.249 20.625 -702 o702 2.50 
7 31 51 7521 161 7511 1241 
87.280 87.280 41.249 20.625 -702 -702 2.50 
a 41 61 75221 7 7511 1252 
87.280 87.280 41.249 20.625 oToy -701 2.50 
7 51 71 7521 181 7512 1261 
87.280 87.280 41.249 20.625 -702 -702 2.50 
7 61 81 7521 191 7511 1272 
87.280 87.280 41.249 20.625 -7021 2701 2.50 
a 71 91 7521 ; 201 7511 1281 
87.280 87.280 41.249 20.625 -702 701 2.50 
7 81 101 7521 211 7511 1291 
87.280 87.280 41.249 20.625 -7021 e700 2.50 
7 91 111 7521 221 7511 1301 
87.280 87.280 41.249 20.625 701 S702 2.50 
a 101 121 7521 231 7511 1311 
87.280 87.280 41.249 20.625 -702 701 2.50 
7 211 7541 75221 241 7511 1321 
87.280 242.445 41.249 20.625 .701 -702 2.50 
7 121 141 11 251 7511 1331 
87.280 87.280 20.625 20.625 Ey Hop! .702 2.50 
7 131 151 22 261 7511 1341 
87.280 87.280 20.625 20.625 -701 -701 2.50 
7 142 162 31 271 7511 1351 
87.280 87.280 20.625 20.625 -7021 -701 2.50 
7 152 171 ai 281 7511 1361 
87.280 87.280 20.625 20.625 -702 701 2.50 
| 161 181 51 291 7511 1371 
87.280 87.280 20.625 20.625 -701 -701 2.50 
7 171 191 61 301 7511 1381 
87.280 87.280 20.625 20.625 -702 -701 2.50 
7 181 201 uo 311 7511 1391 
87.280 87.280 20.625 20.625 701 -701 2.50 
7 192 211 81 322 7511 1402 
87.280 87.280 20.625 20.625 -701 #7 Ol 2.50 
v 201 222 91 331 7511 1411 
87.280 87.280 20.625 20.625 -701 -701 2.50 
7 211 : 231 1021 341 7511 1421 
87.280 87.280 20.625 20.625 -701 e701 2.50 
Z 221 241 111 351 7511 1431 
87.280 87.280 20.625 20.625 270) .7021 2.50 
7 231 251 121 361 7511 1441 
87.280 87.280 20.625 20.625 -701 e701 2.50 
7 241 262 131 371 7511 1451 
87.280 87.280 20.625 20.625 -701 -701 2.50 
7 251 271 142 381 7511 1461 
87.280 87.280 20.625 20.625 .702 -701 2.50 
7 261 281 1521 391 7511 1471 
Figure 9. Partial output data file of PCB terminal session. 


were compared with the results obtained using the analyzer’s manual mode] builder, 
THANSS, and checked for discrepancies. 
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Pepe wc. | 
THIS 1$ THE OUTPUT OF YHE THERMAL ANALYZER 
Tempereturee. Ges 


Mode 1 2 | 4 ) é 
let Tap @s.o¢ e$.o¢ @s.eoc¢ @%.00 es .00 es .6¢ 
Mode 7 8 

det Tao 65.06 @s.0¢ 

Wode 1 2 3 4 5 6 
Mew Tae es .86 es.ec es.e¢ BS .80 @s.8¢ es.8¢ 
Mew-Oled 27998 7998 7998 7998 22997 77997 
mode 7 8 

ew Tap es .6¢ es .6¢ 

wew-Ole 27997 -7998 

Node | 2 3 4 3 é 
Nee Tap OS .%% 65.9% Os .96 Os. 96 es . 96 Os .96 
mewCle elec? 1606 siiieicic ShOSE eleoc 31606 
mode 7? 8 

Mew Tap OS .96 OS .%6 

mew Cia ~l6¢¢ SNK EE 

mode | 2 3 é 5 © 
Mew Tar @6.ec¢ 6.00 66 .o¢ @c.ec @6.0¢ 86.6¢ 


Mew-Cle¢ eo OOSGE-03 o6C36E-01 o40S8E-0) o6061E-01 -4063E-C) »6066F-6) 


Mode ? 8 
MewOle6 e604SE-60) »4046F-0) 


THIS 1S THO CUTIPLUT OF THE TRERMAL ANALYZER 
TEMPER ATUSES. CEGF 


_ 


6 3 lcs c0 G 1CE 96 5 106.06 6 1Co.03 





Figure 10. Output data file of TASS thermal! analyzer. 


E. MODEL BUILDER LIMITATIONS AND POSSIBLE IMPROVEMENTS 
Although the model builder provides the user with a significant advantage in terms 

of time and effort necessary to model a printed circuit board, there are at present certain 

limitations to the program. There is also significant potential for improvements that 


could enhance the capabilities of the current version of the model builder. 
1. PCB Model Builder Limitations 


1. The model builder is restricted to a maximum of four copper layers due to con- 
ventional personal computer operating system limitations. Current versions of 
MS-DOS ® do not allow the user to access extended memory. Access to extended 
memory would allow for a larger copper layer capability for the model. 


2. The model builder will not run when using MS-DOS ®, version 4.01. When run- 
ning PCB on a computer operating with MS-DOS ®, version 4.01, the computer 
will display an “insufficient neemory” message. Writing more compact code for the 
model builder could overcome this problem. The program will run on any other 
version of MS-DOS ®. 


3. The model builder does not account for convective and radiative heat transfer. 


4. The model builder only accepts heat inputs from the top surface of the printed 
circuit board. 


2. Possible Improvements to the Thermal Model Builder. 


1. The present version of the model builder requires extensive amounts of FORTRAN 
code. An updated model builder could be designed using a more compact and ef- 
ficient higher level language, such as C++. Smaller, more efficient source code 
would allow for incorporating more copper layers into the model builder. 


2. The menu driven queries used in the thermal analyzer could be improved by tn- 
corporating assembly language routines enabling the use of a mouse. Attempts 
were made to introduce the use of mouse driven menus; however, memorv handling 
lunitations in the Microsoft ® linker, version 3.55, prevented the use of a mouse. 


3. Allowing for the injection of heat not only from the top surface of the printed cir- 
cuit board, but also from other external surfaces (to include the heat being trans- 
ferred through the wedges supporting the oe circuit board), 1s another 
potential unprovement. 


4. Adding graphics capabilities to the model builder would enable the user to model 
the printed circuit board without the need for preliminary drawings or sketches; 
however, adding this feature would cause the model builder executable code to be- 
Pome excessively large. 


VI. CONCLUSIONS 


The purpose of this thesis is to develop a printed circuit board thermal analvsis 
model builder that will effectively interact with thermal analysis software used by the 
Naval Postgraduate School. [Ref. 9] 

The current version of tlhe model builder enables the user to choose from sixteen 
possible printed circuit board configurations, ranging from one to four copper lavers. 
Other features of the inodel builder include the ability to set the physical characteristics 
of the board including dimensions, thermal conductivitics, percent copper coverage, 
board temperatures, and heat input modes. 

The main goal of the model builder is to enable the user to model a printed circuit 
board with minimum effort, and produce a thermal analysis input data file. Using the 
model builder, PCB, relieves the user of the tedious, time-consuming manual data entry 
required of the thermal analyzer’s current model builder, THANSS. [Ref. 9] 

The expansion possibilities for the model builder are significant. Adding more cop- 
per layers and nodes would enable the user to model more sophisticated printed circuit 
boards. Other potential features include the addition of graphics capabilities and the 
introduction of user-friendly peripherals, such as mouse drivers or digitazing pads. As 
the level of complexity of future model builders increases, it will become necessary to 
Import the source code from FORTRAN compilers linuted to 640 Kilobytes of memory, 
to 32 bit FORTRAN compilers allowing use of all RAM. The source code for PCB 
contains over 5000 lines of code, effectively reaching the memory limits of both the 
compiler and personal computer. 

The possible applications for the model builder are substantial. As the level of 
complexity of printed circuit boards increases, there is a real need for a tool that allows 
the designer to efficiently perform thermal analysis of printed circuit boards during the 
design process. The model builder employed must allow the designer the necessary 
flexibility needed to model a board that meets the desired design criteria. The PCB 
model builder is the first step in creating a tool that enables the designer to effectively 


build a model that meets specified requirements. 
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APPENDIX PROGRAM LISTING FOR PCB MODEL BUILDER 


$ LARGE 
C 

C 

oe TITLE: MODEL BUILDER -- MAIN PROGRAM 
C AUTHOR: LT STEVE GLASER 

C ‘DATE: 09 JUL 1991 

C COMPILER: MICROSOFT VERSION 4.01 

C LINKER: MICROSOFT VERSION 3.55 

C 

C DEFINE REAL VARIABLES 

& 


Poiteaw  WiGinOKE obsow ,l1G1;1201,1T202,13L1,13L2,1T3L3,T4L1,T4L2, 
+T4L3 ,T4L4 ,TiW1,T2W1,T2W2 ,T3W1,T3W2,T3W3 ,T4W1,T4W2 ,T4W3 ,T4W4 , AREA,U 
+N ,VOL1,VOL2 ,VOL3 , VOL4 ,UKC , IBT,UPRT, LWRT,RT,LT,FT,BT,THEAT , THPN , AHE 
+AT , NHEAT , DELX , DELY 


COMMON IBT,UPRT,LWRT,RT,LT,FT,BT,UT,ANSWER 


INTEGER NWIDE ,NDEEP ,NPL,NUM,NUMA,CH,H 


C 
@ 
C DEFINE ALL ONE CHARACTER VARIABLES 
@ 
e 
C 

CHARACTER*1 SELECT,ANS ,ANSL,OVR, ANSWER, ANSN, ANSA, LAY1, LAY2, LAY3,LA 

+Y4,ANSE,ANSEL, ANSEW, ANSET, ANSEK , ANSQ ,WORL, ANSE1, ANST1,ANSK1,ANSQ1, 
+SPEVAL, ANSTEM , ANSTEB , ANSTEU , ANSTEL, ANSTER , ANSLFT , ANSFRT , ANSBCK 

e 
@ 
C CHARACTER VARIABLES OF MORE THAN ONE POSITION 
C 
C 

CHARACTER UK*11,UT*1,UH*13 , UAH*6 ,UL*2 , UW*2 
C 
@ 
C DEFINE MATRICES 
C 

REAL HEAT( 100,100) ,COEF( 1000, 12) 
C 
C 
@ 

INTEGER IH( 1000) ,JH( 1000) ,NCON( 1000, 12) 

C 
C 
C PROVIDE THE USER WITH AN INTRODUCTORY STATEMENT. 
C 
C 
88888 CALL CLS 
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op) (op) py (qe emp ol ap) Cap elem et ap )) (ais Ja gp | (cap J ca al eg a gp Bs ge J oe I > JE a I a> > J > J a> J DG I ep I > I ep ep Bp I a 1 a op) | ap ge) (>) sel (sp) 


001 


002 


WRITE(*, 7001) 


FORMAT(///////,' "Ws 
+! THIS PROGRAM WAS WRITTEN TO INTEGRATE WITH EXISTING ',/, 

+! THERMAL ANALYSIS SOFTWARE AND TO REDUCE THE AMOUNT 7 

+" OF TIME REQUIRED FOR DATA ENTRY.',//, 

+! WOULD YOU LIKE AND OVERVIEW OF THE PROGRAM PRIOR TO ',/, 

+! BEGINNING? ENTER Y FOR YES AND N FOR NO: 2X 
READ(*,7002) OVR 

FORMAT(A1) 


VARIABLE, CONSTANT, AND STRING DEFINITION 
PHYSICAL CHARACTERISTICS 
CL,EL - EPOXY AND COPPER LENGTHS 
CW,EW - EPOXY AND COPPER WIDTHS 
UKE,UKC - EPOXY AND COPPER THERMAL CONDUCTIVITIES 
DELE, T1L1 THROUGH T4L4 - THICKNESS OF EPOXY AND COPPER LAYERS 
DELX - CL/NDEEP 
DELY - CW/NWIDE 
NPL - NUMBER OF NODES PER LAYER 
NWIDE - NUMBER OF NODES WIDE 
NDEEP - NUMBER OF NODES DEEP 
UL - UNITS OF LENGTH (SI OR ENGLISH) 
UK - UNITS OF THERMAL CONDUCTIVITY (SI OR ENGLISH) 


INITIAL AND AMBIENT TEMPERATURES 


IT - INITIAL BOARD TEMPERATURE 

LT - LEFT SIDE AMBIENT TEMPERATURE 

RT - RIGHT SIDE AMBIENT TEMPERATURE 

FT - FRONT AMBIENT TEMPERATURE 

BT - BACK AMBIENT TEMPERATURE 

UPRT - UPPER AMBIENT TEMPERATURE 

LWRT - LOWER AMBIENT TEMPERATURE 

UT - UNITS OF TEMPERATURE (CENTIGRADE OR FAHRENHEIT) 


HEAT INPUT 


THEAT - TOTAL INJECTED HEAT 
THPN - TOTAL HEAT PER NODE 
AHEAT - AVERAGE HEAT OVER A GIVEN SURFACE 
NHEAT - NODE PER HEAT INJECTED NODE PER NODE 
UH, UAH - UNITS OF HEAT (SI OR ENGLISH) 
NUM, NUMA, CH, H - DUMMY VARIABLES 
TOTNOD, NN, NC, N - VARIABLES USED TO ALLOW NODAL HEAT INPUT 
IH, JH - VECTORS USED TO CORRELATE NODE NUMBER WITH 
MATRIX POSITION 
HEAT - MATRIX USED TO HOLD HEAT INPUTS 


COEFFICIENT DEFINITIONS Y IMPLIES WIDTH, AND X IMPLIES DEFT 
Z IMPLIES HEGGuY 
C IMPLIES COPPER, E IMPEEESS Er Gis 
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801 


EYLR - PROVIDES COEFFICIENT FOR LEFT OR RIGHT EDGE NODES 
TO THE EXTERNAL NODE, EYLR IMPLIES EPOXY. 
CYLR - SAME AS ABOVE, COPPER LAYER 


EXFB - PROVIDES COEFFICIENT FOR FRONT AND BACK EDGE NODES 
TO THE EXTERNAL NODE. 
CXFB - SAME AS ABOVE FOR COPPER LAYER 


Es SUN ERNAL COEFFICIENT IN THE Y DIRECTION, EPOXY 
CYY - INTERNAL COEFFICIENT IN THE Y DIRECTION, COPPER 


EXX - INTERNAL COEFFICIENT IN THE X DIRECTION, EPOXY 
CXX - INTERNAL COEFFICIENT IN THE X DIRECTION, COPPER 


SZteme COM? LOLENT FOR EPOXY TO UPPER EXTERNAL@USE 
EZG - COEPPICIENT FOR COPPER TO" EPOAY 
PCC - PERCENT COPPER COVERAGE 


ALL ONE CHARACTER STRINGS ARE SIMPLE 'YES' OR 'NO' TYPE 
VARIABLES. 


CALL INTRODUCTION 


IF (OVR. EQ. 'Y') THEN 
CALL INTRO 

ELSEIF (OVR. EQ. 'N') THEN 
GOTO 899 

ELSE 

GOTO 88888 

ENDIF 


SUPPLY A LIST OF ACCEPTABLE LAYER LEVELS 


BEGIN PROGRAM 


CALL CLS 


WRITE(* , 800) | 
FORMAT(///,' HAKAAAIAAII KIKI RAITT ARTE 


pte! | / 

! Pete Fede de Fee Pe PII RT TORII ICT RICH AIR TIT IARI REAR EEERERS | ; / . 
fer es 
: Ke PLEASE SELECT UPPER CASE LETTERS en 


+ 
+ 

+ 

+! fee PRIOR TO BEGINNING say 
it oe wee! 

a 

+ 


: KAKA AAAK AKAIKE RIAA KERIKERI ERIE ERERERERERER |S | 
: HCH AKIA IK AKIRA KIKI IK RIKER ER EERE RERERERE |) 


WRITE(*, 801) 

FORMAT(//,' ae 

+" PRIOR TO ENTERING DATA INTO THIS PROGRAM ENSURE THAT ',/, 
+" YOU HAVE A DRAWING OF YOUR DESIGN AND ALL PERTINENT ',/, 
i DATA. ee 
+! PRESS <ENTER> TO CONTINUE v 2X) 
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READ(* , 802) ANSWER 
802  FORMAT(A1) 


C 

C 

G 

CALLGie 

357. WRITE(*, 367) 

367 FORMAT(///////,' THIS PROGRAM IS CAPABLE OF OPERATIONS IN EI 
+THER SI OR ',/, 
+! ENGLISH UNITS. AFTER THE SELECTION OF THE UNITS, ALL ',/, 
+! ENTRIES MUST BE COMPATIBLE. PLEASE MAKE YOUR SELECTION. '/ 
+/,' S FOR SI NOTATION’ ,/, 
+! E FOR ENGLISH NOTATION: ',2X, ) 


READ(*,368) ANSN 
368 | FORMAT(A1) 
€ CHECK FOR CORRECT UNIT SELECTION 
CG 
3333 IF(ANSN.EQ. 'S') THEN 
WRITE(*, 369) 


369  FORMAT(////,' YOU HAVE SELECTED SI NOTATION. ') 
ELSEIF (ANSN. EQ. 'E') THEN 
WRITE(*, 370) 
370  FORMAT(////,' YOU HAVE SELECTED ENGLISH NOTATION. ' ) 
ELSE 
CALL ACis 
GOTO 357 
ENDIF 
6 
3711 WRITE(*,3710) 
3710 FORMAT(//,' IS THIS THE DESIRED SELECTION? ENTER Y FOR YES A 
+ND N FOR NO: ‘', ) 


READ(*,372) ANSA 
372 FORMAT(A1) 
C 
IF(ANSA. EQ. 'Y') THEN 
GOTO 378 
ELSEIF(ANSA. EQ. 'N') THEN 
GALL CLs 
GOTO 357 
ELSE 
CALLM@ES 
GOTO 3333 
78 ENDIF 


CALE GES 

SPEVAL='A' 

WRITE(*, 301) 

301 FORMAT(//////,' THE PCB CAN HAVE UP TO FOUR COPPER LAYERS’ ) 
WRITE(*, 302) 

302 FORMAT(//,' HOW MANY COPPER LAYERS DO YOU DESIRE?’ ) 
WRITE(*, 303) 

303 FORMAT(/,' PICK A NUMBER 1 TO 4: ',2X, ) 
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READ(* , 304) SELECT 
304  FORMAT(A1) 


C 
C 
CALL CLS 
IF (SELECT. NE. '1'. AND. SELECT. NE. '2'. AND. SELECT. NE. '3'. AND. 
+SELECT. NE. '4') THEN 
GOTO 7 
ELSE 
10 WRITE(* ,305 )SELECT 
305 FORMAT(//////,' YOU SELECTED ',Ai,' COPPER LAYER(S) FOR THE PCB 
sae! / ) 
ENDIF 
C 
WRITE(* , 306) 
306 FORMAT ( IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FO 
+R NO: ox, ) 


READ(*,307)ANS 
307. FORMAT(A1) 
C 
C 
IF(ANS. EQ. 'N') THEN 
GOTO 7 
ELSE 
CONTINUE 
ENDIF 
IF(ANS. EQ. 'Y') THEN 
GOTO 9 
ELSE 
CALL CLS 
GOTO 10 
ENDIF 


CALL CLS 


aac 2 ©) 


IF(SELECT. EQ.'1') THEN 
GOTO 123 

ELSE 
CONTINUE 

ENDIF 

IF( SELECT. EQ. '2')THEN 
GOTO 124 


ENDIF 

IF( SELECT. EQ. '3')THEN 
GOTO 125 

ELSE 


5] 


C 

C 
ab 
Ma 
308 


309 


310 


3334 


S| 


i 
S12 


S13 
C 
C 


IF( SELECT. EQ. ‘4° )THEN 
GOTO 126 

ELSE 
CONTINUE 

ENDIF 
GOTO 7 


CAL CLS 

WRITE(* , 308) 

FORMAT(//////,' YOU SELECTED 1 COPPER LAYER, GIVING YOU THE FOLLO 
+WING ALTERNATIVES: ' ) 

WRITE(*, 309) 

FORMAT(//,' 1. 10 BY 36 PROVIDES A 1:3.6 RATIO WITH 720 NODES. 
+',/,' 2. 10 BY 35 PROVIDES A 1:3.5 RATIO WITH 700 NODES.',/, 
+' 3. 12 BY 25 PROVIDES A 1.2:2.5 RATIO WITH 600 NODES.',/, 
+' 4 15 BY 20 PROVIDES A 1.5:2.0 RATIO WITH 600 NODES.',//, 
+' PLEASE SELECT A NUMBER 1 THROUGH 4: 2s DK) 

READ(* ,310)LAY1 

FORMAT(A1) 

IF (LAY1.NE. '1'. AND. LAY1. NE. '2'. AND. LAY1. NE. 
+'3'. AND. LAY1. NE. '4') THEN 

GOTO 11 

ELSE 

CALERGLS 

WRITE(*,311)LAY1 

FORMAT(//////,' YOU SELECTED NUMBER 'Al' OF THE FOLLOWING 4 ALTER 
+NATIVES: | ,///, 
+' 1. 10 BY 36 PROVIDES A 
+' 2. 10 BY 35 PROVIDES A 
+' 3, 12 BY 25 PROVIDES A 
+' 4. 15 BY 20 PROVIDES A 
ENDIF 

WRITE(* , 312) 

FORMAT(' IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FO 
+R NO: YO 

READ(*,313)ANSL 

FORMAT(A1) 


3.6 RATIO WITH 720 NODES.',/, 

3.5 RATIO WITH 700 NODES. ',/, 
2:2.5 RATIO WITH 600 NODES.',/, 
5:2.0 RATIO WITH 600 NODES. ',///) 


IF(ANSL. EQ. 'N') THEN 
cone i 


ENDIF 
IF(ANSL. EQ. 'Y') THEN 
GOTO 12 
ELSE 
GOTO 3334 
ENDIF 
GOTO 1234 
CALL CLS 
WRITE(* , 314) 
FORMAT(//////,' YOU SELECTED 2 COPPER LAYERS, GIVING YOU THE FOLL 
+OWING ALTERNATIVES: ' ) 


a2 


B15 


516 


5335 


o1/ 


Ws 
318 


519 
C 
C 


16 


7 
25 
320 


o21 


Sapa 


WRITE(*, 315) 
FORMAT(//,' 1. 10 BY 18 PROVIDES A 1:1.8 RATIO WITH 720 NODES. 
+',/,' 2. 9 BY 20 PROVIDES A 1:2.11 RATIO WITH 720 NODES.',/, 
+' 3. 8 BY 20 PROVIDES A 2:5 RATIO WITH 640 NODES. ',/, 
+' 4. 10 BY 15 PROVIDES A 1:1.5 RATIO WITH 600 NODES.',//, 
+' PLEASE SELECT A NUMBER 1 THROUGH 4: Ok) 

READ(* ,316)LAY2 

FORMAT( Al) 

IF (LAY2. NE. '1'. AND. LAY2. NE. '2'. AND. LAY2. NE. 
+'3'. AND. LAY2. NE. '4') THEN 

GOTO 14 

ELSE 

CAE CLS 

WRITE(* ,317)LAY2 

FORMAT(//////,' YOU SELECTED NUMBER 'Al' OF THE FOLLOWING 4 ALTER 
+NATIVES: ',///, 
+' 1, 10 BY 18 PROVIDES A 1:1.8 RATIO WITH 720 NODES.',/, 
+' 2. 9 BY 20 PROVIDES A 1:2.11 RATIO WITH 720 NODES.',/, 
+' 3. 8 BY 20 PROVIDES A 2:5 RATIO WITH 640 NODES.',/, 
+' 4. 10 BY 15 PROVIDES A 1:1.5 RATIO WITH 600 NODES.',///) 
ENDIF 

WRITE(* ,318) 

FORMAT(' IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FO 
+R NO: wa2x; +) 

READ(* , 319) ANSL 

FORMAT(A1) 


IF(ANSL. EQ. 'N') THEN 
GOTO 14 

ELSE 
CONTINUE 

ENDIF 

IF(ANSL. EQ. 'Y') THEN 
GOTO 16 

ELSE 
BOTOR 235 

ENDIF 

EOTO 1235 


GALL CLS 
WRITE(*, 320) 

FORMAT(//////,' YOU SELECTED 3 COPPER LAYERS, GIVING YOU THE FOLL 
+OWING ALTERNATIVES: ' ) 

WRITE(* , 321) 

FORMAT(//,' 1. 5 BY 24 PROVIDES A 1:4.8 RATIO WITH 720 NODES. 
+',/,' 2. 10 BY 10 PROVIDES A 1:1 RATIO WITH 600 NODES.',/, 
+' 3, 12 BY 10 PROVIDES A 1.2:1 RATIO WITH 720 NODES.',/, 
+' 4, 8 BY 15 PROVIDES A 1:1.875 RATIO WITH 720 NODES.',//, 
+' PLEASE SELECT A NUMBER 1 THROUGH 4: Ox. ) 
READ(* , 322) LAY3 

FORMAT(A1) 

IF (LAY3. NE. '1'. AND. LAY3. NE. '2'. AND. LAY3. NE. 
+'3'. AND. LAY3. NE. '4') THEN 

GOTO 17 
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3336 
3)/,.8 


18 


20 
126 
326 


27 


328 


2557 


329 


7a \\ 
a0 


ELSE 

CALL CES 

WRITE(* , 323) LAY3 

FORMAT(//////,' YOU SELECTED NUMBER 'Al1' OF THE FOLLOWING 4 ALTER 
+NATIVES: ',///, 
+' 1. 5 BY 24 PROVIDES A 1:4.8 RATIO WITH 720 NODES.',/, 
+' 2, 10 BY 10 PROVIDES A 1:1 RATIO WITH 600 NODES.',/, 
+' 3. 12 BY 10 PROVIDES A 1. 2:1 RATIO WITH 720 NODES. ',/, 
+' 4. 8 BY 15 PROVIDES A 1:1.875 RATIO WITH 720 NODES. ,///) 

ENDIF 

WRITE(* , 324) 

FORMAT ( IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FO 
+R NO: 2X a) 

READ(* , 325) ANSL 

FORMAT(A1) 


IF(ANSL. EQ. 'N') THEN 
GOTO 17 
ELSE 
CONTINUE 
ENDIF 
IF(ANSL. EQ. 'Y') THEN 
GOTO 18 
ELSE 
GOTO 3336 

ENDIF 

GOTO 1236 

CALLeCES 

WRITE(* ,326) 

FORMAT(//////,' YOU SELECTED 4 COPPER LAYERS, GIVING YOU THE FOLL 
+OWING ALTERNATIVES: ' ) 

WRITE(* , 327) 

FORMAT(//,' 1. 5 BY 18 PROVIDES A 1:3.6 RATIO WITH 720 NODES. 
+',/,' 2. 9 BY 10 PROVIDES A 1:1.11 RATIO WITH 720 NODES.',/, 
+' 3. 8 BY 10 PROVIDES A 1:1.25 RATIO WITH 640 NODES.',/, 
+' 4. 13 BY 6 PROVIDES A 1:1 RATIO WITH 624 NODES. ',//, 
+' PLEASE SELECT A NUMBER 1 THROUGH 4: wre) 

READ(* ,328)LAY4 

FORMAT(A1) 

IF (LAY4. NE. '1'. AND. LAY4. NE. '2'. AND. LAY4. NE. 
+'3', AND. LAY4. NE. '4') THEN 

GOTO 20 

ELSE 

GALL CLS 

WRITE(*, 329) LAY4 

FORMAT(//////,' YOU SELECTED NUMBER 'Al' OF THE FOLLOWING 4 ALTER 
+NATIVES:',///, 


+' 1. 5 BY 18 PROVIDES A 1:3.6 RATIO WITH 720 NODES.',/, 
+' 2. 9 BY 10 PROVIDES A 1:1.11 RATIO WITH 720 NODES. ,/, 
+' 3, 8 BY 10 PROVIDES A 1:1.25 RATIO WITH 640 NODES.',/, 
+' 4. 13 BY 6 PROVIDES A 1:1 RATIO WITH 624 NODES. ',///) 
ENDIF 

WRITE(*, 330) 


FORMAT(' IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FO 
+R NO: ne ON) 


ob 


2. 


1234 


1235 


1236 


READ(* ,331)ANSL 
FORMAT(A1) 


IF(ANSL. EQ. 'N') THEN 
GOTO 20 

ELSE 
CONTINUE 

ENDIF 

IF(ANSL. EQ. 'Y') THEN 
GOTO 22 

ELSE 
GOEO 3337 

ENDIF 

SOTO 1237 


CONTINUE 
IF( LAY1. EQ. '1')THEN 
NWIDE=10 
NDEEP=36 
NPL=360 
ELSEIF(LAY1. EQ. '2')THEN 
NWIDE=10 
NDEEP=35 
NPL=350 
ELSEIF( LAY1. EQ. '3' )THEN 
NWIDE=12 
NDEEP=25 
NPL=300 
ELSEIF( LAY1. EQ. '4' )THEN 
NWIDE=15 
NDEEP=20 
NPL=300 
ENDIF 
Comewi241 
CONTINUE 
IF( LAY2. EQ. '1')THEN 
NWIDE=10 
NDEEP=18 
NPL=180 
ELSEIF( LAY2. EQ. '2' )THEN 
NWIDE=9 
NDEEP=20 
NPL=180 
ELSEIF( LAY2. EQ. '3' )THEN 
NWIDE=8 
NDEEP=20 
NPL=160 
ELSEIF( LAY2. EQ. '4')THEN 
NWIDE=10 
NDEEP=15 
NPL=150 
ENDIF 
GOTO 1241 
CONTINUE 
IF( LAY3. EQ. '1')THEN 
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1237 


qaqa 


5400 


NW IDE=5 

NDEEP=24 

NPL=120 
ELSEIF( LAY3. EQ. '2' ) THEN 

NWIDE=10 

NDEEP=10 

NPL=100 
ELSEIF( LAY3. EQ. '3' ) THEN 

NWIDE=12 

NDEEP=10 

NPL=120 
ELSEIF( LAY3. EQ. '4' ) THEN 

NWIDE=8 

NDEEP=15 

NPL=120 
ENDIF 

GOTO 1241 

CONTINUE 
IFC LAY4. EQ. '1')THEN 

NWIDE=5 

NDEEP=18 

NPL=90 
ELSEIF( LAY4. EQ. '2' )THEN 

NWIDE=9 

NDEEP=10 

NPL=90 
ELSEIF( LAY4. EQ. '3' ) THEN 

NWIDE=8 

NDEEP=10 

NPL=80 
ELSEIF(LAY4. EQ. '4' )THEN 

NWIDE=13 

NDEEP=6 

NPL=78 
ENDIF 

GOTO 1241 


EPOXY CHARACTERISTICS 


CALL CLS 

WRITE(*, 8000) 

FORMAT ( / //, | eH ARE KARR IRR IKKE RRR RER REE 
PHAHKAAAKKAIAKKK AK I | / 
+! reredededesesededeicdede dese reser resend KER ENE POXY CHARACTERISTICS errricicttededededededed 
pre! | / 
$e! dedededededededede de dete RKTT IER RTI REITER ERR EERE RRIRRIRERE 
picickeaeraesed! 7 |) 


PROVIDE CORRECT UNIT ABBREVIATIONS 


IF(ANSN. EQ. 'S') THEN 
WRITE(*,5400) 
FORMAT(' ALL ENTRIES ARE IN SI NOTATION. ',/) 
UL='cm' 
UK='Watts/cm/C' 
Oi=" Ce 
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5401 


7 99 


UW=' gm' 
ELSEIF(ANSN. EQ. 'E') THEN 


WRITE(*, 5401) 
FORMAT(' ALL ENTRIES ARE IN ENGLISH NOTATION’ ,/) 
UL=' tan 
UK='Btu/hr/F' 
UT=' yr! 
UW='oz' 
ENDIF 


WRITE(* ,5402) UL 


FORMAT(/,' ENTER EPOXY LAYER LENGTH (',A2,'): ',2X, ) 

READ *,EL 

WRITE(*,5403) UL 

FORMAT(/,' ENTER EPOXY LAYER WIDTH (',A2,'): ‘',2X, ) 

READ *,EW 

WRITE(* ,6403) 

FORMAT(//,' THE EPOXY LAYER THICKNESS WILL DEFAULT TO 0.0625 in ( 
+0.15875 cm)',/,' DO YOU WANT TO CHANGE THE EPOXY THICKNESS? (Y OR 
TN 52k, ) 

READ(*,5799) ANSQ 

FORMAT(A1) 

IF(ANSQ. EQ. 'Y') THEN 

GOTO 5798 


ELSEIF(ANSQ. EQ. 'N'. AND. ANSN. EQ. 'S') THEN 
DELE=0. 15875 


GOTO 4338 
ELSEIF(ANSQ. EQ. 'N'. AND. ANSN. EQ. 'E') THEN 
DELE=0. 0625 
GOTO 4338 
ELSE 
GOTO 4337 
ENDIF 
5798 WRITE(*,5404) UL 
5404 FORMAT(/,' ENTER EPOXY LAYER THICKNESS (',A2,'): ',2X, ) 
READ *,DELE 
C 
4338 WRITE(*,5405) UK 
5405 FORMAT(/,' ENTER EPOXY LAYER THERMAL CONDUCTIVITY (',A10,'): ‘',2 
18, ) 
READ *,UKE 
5459 CALL CLS 
g 
C MAKE CHANGES OR CORRECTIONS TO EPOXY ENTRIES 
C 
WRITE(* , 5406) 
5406 FORMAT(///,' YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE EPOXY 
+ LAYER: ',/,) 
5407 WRITE(*,89343) EL,UL 
89343 FORMAT(/,' 1.) LENGTH: "MROCGw INA? , ) 
WRITE(* ,89344) EW,UL 
89344 FORMAT(/,' 2.) WIDTH: ' F9.4,1X,A2, ) 
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WRITE(*, 89345) DELE,UL 


89345 FORMAT(/,' 3.) THICKNESS: ',F9.4,1X,A2, ) 


WRITE(*,89346) UKE,UK 


89346 FORMAT(/,' lam) kk: " F9.4, IX SAl 15 


C 


5408 


5409 


5422 


C 
5440 


5442 


5443 


5444 


5445 


WRITE(* ,5408) 
FORMAT(/,' DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR YES 


+ AND N FOR NO: "SD 


READ(* ,5409) ANSE 
FORMAT(A1) 


IF(ANSE. EQ. 'Y') THEN 
CALL CLS 
WRITE(*,5410) EL,UL 
FORMAT(////,' | THE CURRENT ENTRY FOR LENGTH IS: ',F9.4,1X,A2) 
WRITE(*, 5411) 
FORMAT(/,' WOULD YOU LIKE TO CHANGE THE LENGTH? (Y OR N): 


+ 52k) 


READ(* ,5412) ANSEL 

FORMAT(A1) 

PRINT * 
IF( ANSEL. EQ. 'Y') THEN 
WRITE(* ,5422) UL 

FORMAT(/,' ENTER THE EPOXY LENGTH (° ,A2y =) Seon 

READ *,EL 
ELSEIF( ANSEL. EQ. 'N') THEN 
GOTO 5440 
ELSE 
GOTO 5441 
ENDIF 


CALL CLS 
WRITE(* 5442) EW,UL 
FORMAT(////,' THE CURRENT ENTRY FOR WIDTH IS: ',F9.4,1X,A2) 
WRITE(* , 5443) 
FORMAT(/,' WOULD YOU LIKE TO CHANGE THE WIDTH? (Y OR N): 


ge yD 


READ(* , 5444) ANSEW 

FORMAT(A1) 

PRINT * 
IF( ANSEW. EQ. 'Y') THEN 
WRITE(* ,5445) UL 

FORMAT(/,' ENTER THE EPOXY WIDTH (',A2,'): ‘',2X, ) 

READ *,EW 
ELSEIF( ANSEW. EQ. 'N') THEN 
GOTO 5446 
ELSE 
GOTO 5440 
ENDIF 


GALLACLS 
WRITE(*,5447) DELE,UL 
FORMAT(////,° THE CURRENT ENTRY FOR THICKNESS IS: ',F9.4,1X,A2) 
WRITE(* , 5448) 
FORMAT(/,' WOULD YOU LIKE TO CHANGE THE THICKNESS? (Y OR N 


58 


5 2K50) 
READ(* ,5449 ) ANSET 
5449 FORMAT(A1) 
PRINT * 
IF(ANSET. EQ. 'Y') THEN 
WRITE(* 5450) UL 
5450 FORMAT(/,' ENTER THE EPOXY THICKNESS (',A2,'): ',2X, ) 
READ *,DELE 
ELSEIF( ANSET. EQ. 'N’) THEN 
GOTO 5451 
ELSE 
GOTO 5446 
ENDIF 
C 
5451 GALES CLS 
WRITE(* ,5452) UKE,UK 
5452 FORMAT(////,' THE CURRENT ENTRY FOR THERMAL CONDUCTIVITY IS: ', 
+F9.4,1X,A10) 
WRITE(* ,5453) 


5453 FORMAT(/,' WOULD YOU LIKE TO CHANGE THE THERMAL CONDUCTIVIT 
+¥? (Y OR N): ',2X, ) 
READ(* ,5454) ANSEK 
5454 FORMAT(A1) 
PRINT * 


IF( ANSEK. EQ. 'Y') THEN 
WRITE(*,5455) UK 

5455 FORMAT(/,' ENTER THE THERMAL CONDUCTIVITY (',A10,'): ',2X, ) 
READ *,UKE 
ELSEIF( ANSEK. EQ. 'N’) THEN 
GOTO 5456 
ELSE 
GOTO 5451 
ENDIF 

5456 CAE CLS 
WRITE(* ,5457) 

5457 FORMAT(///,' YOU HAVE MADE THE FOLLOWING CORRECTIONS TO THE 

+ EPOXY ENTRIES: ',//) 
GOTO 5407 


ELSEIF(ANSE. EQ. 'N' ) THEN 
GOTO 5458 

ELSE 

GOTO 5459 

ENDIF 


CALL SUBROUTINE COPPER 


WmWanraIranrIinanana 


foo) OALL COPPER(CEW,EL,ANSN,UK,SELECT,T1L1,T2L1,T2L2,T3L1,T3L2 ,T3L3,T4L 
+1,T4L2 ,T4L3 ,T4L4,SPEVAL,UKC,CL,CW) 


C2 C2 
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CALL SUBROUTINE PCBS1 


1180 CALL PCBS1 


CALL SUBROUTINE PCBS2 


> Sp > Dw Ee ae ge ice [ge 


CALL PCBS2( THEAT , THPN , AHEAT , NHEAT ,NPL,NWIDE ,NDEEP ,HEAT,IH,JH,ANSN, 
+EL,EW) 


CALL SUBROUTINE PCBS3 


C2 ©) Garca c 


CALL PCBS3(EW,EL,CW,CL,UKE ,DELE,UKC,SELECT ,T1L1,12L1,121L2,T3L ee 
+2 ,T3L3,T4L1,T4L2 ,T4L3 ,T4L4 ,NWIDE,NDEEP,NPL, IH, JH,HEAT ,COEF ,IBT,UPR 
+T, LWRI RIG. Br ANSN) 

END 

CMa AAA a IRIN IRIE RIN KIKI IK IK KIKI KICKIN HIKARI IRIK KARIERRE IR 


SUBROUTINE INTRO 


GES MODEL BUILDER 

SUBROUT: INTRO 

AUTHOR: LT STEVE GLASER 

DATE: OF Ub Iga 

COMPILER: MICROSOFT VERSION 4.01 
LINKER: MICROSOFT VERSION 3.55 


THIS SUBROUTINE PROVIDES THE USER WITH AN OVERVIEW OF THE THERMAL 
ANALYZER MODEL BUILDER. 


Ca Cai Ca ee) Carta C2 Ca Co Ca C2. Cacca 


COMMON IBT,UPRT,LWRT,RT,LT,FT,BT,UT, ANSWER 


DEFINE ONE VARIABLE CHARACTER VALUES 


CCC) Gr) C2 


CHARACTER*1 ANSWER 


PROVIDE THE USER WITH A PROGRAM OVERVIEW 


Ci Ca Carca Cc) ©) 


60 


C 
C 


7002 FORMAT (Al) 


7003 FORMAT(////// 5" 


7004 


7005 


PALL CLS 

WRITE(*, 7003) 

Teveterastesteseg OVERVIEW KiekKR RT« 
THIS PROGRAM PERFORMS A NODAL ANALYSIS OF A PRINTED 
CIRCUIT BOARD CONTAINING UP TO FOUR COPPER LAYERS 
(WITH EPOXY LAYERS IN BETWEEN). THE OUTPUT CONSISTS 
OF UP TO 720 COEFFICIENTS THAT CONTRIBUTE TO THE 
DETERMINATION OF THE TEMPERATURE DISTRIBUTION 

OF THE PRINTED CIRCUIT BOARD (PCB) WHEN FED INTO THE 
THERMAL ANALYZER. ',//, 

THE FOLLOWING IS AN OUTLINE OF THE MAJOR y. 
SECTIONS OF THIS PROGRAM AND WHAT ENTRIES ARE 
REQUIRED OF THE USER. 


—= = = = == «= 


— == =s = 


PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS 


—= = 


Fete tteteteteteetet 


PRESS <ENTER> TO CONTINUE 
READ(* , 7002) ANSWER 


CALL CLS 
WRITE(*, 7004) 


FORMAT(////// 5" 
A. DATA OUTPUT FILE: 


Ke OVERVIEW wieiiedesedeseiek 


THIS PROGRAM GENERATES AN OUTPUT DATA FILE WHICH IS 
TO BE THE INPUT TO THE THERMAL ANALYZER. 

DURING THE COURSE OF THE PROGRAM, THE USER WILL BE 
ASKED TO PROVIDE A NAME FOR THE OUTPUT FILE. 


WHEN PROMPTED PLEASE ENTER THE NAME OF THE OUTPUT 
DATA FILE. THE DATA FILE NAME SHOULD BE NO LONGER 
THAN SIX LETTERS, AND MAY NOT HAVE ANY SPACES. 


— = = = 


PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS 


—_-— = 


tettteteteteteett 


PRESS <ENTER> TO CONTINUE 
READ(*, 7002) ANSWER 


CALL CLS 

WRITE(*,7005) 

FORMAT( //////,' fekieinieinieiek OVERVIEW tities 
B. STRUCTURE PHYSICAL CHARACTERISTICS 


—_- = 


1. YOU WILL BE ASKED TO SELECT UNIT TYPE 
(SI OR ENGLISH). 

2. THE PRINTED CIRCUIT BOARD IS MAINLY COMPOSED OF 
ALTERNATING COPPER AND EPOXY LAYERS. 
IS DESIGNED TO PROVIDE THE USER WITH A RANGE 
OF 1 TO 4 COPPER LAYERS. IT IS ASSUMED THAT EACH 

COPPER LAYER LIES BETWEEN EPOXY LAYERS. 


-_- = 


= 


tet tte teeter 


THE PROGRAM WILL ALLOW THE USER TO CHOOSE FROM 
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DEPENDING ON THE NUMBER OF COPPER LAYERS SELECTED, 


‘wal 
— 
“— 


-— = ~e we we ww w& >= 2A - = 
we we ewe A RL SR Stew Lew 
i ee 

we we wae 


w wo wo ) 


on 


of 


3 


t 
3 
f 
3 3 
t 
3 3 
f 
3 


— 
we 


THE PROGRAM’ ,/, 


we 
ie 


»/s 
ye 


3 3 


7006 


7007 


7008 


+! FOUR ALTERNATIVE NODAL ASPECT RATIOS. "Se 
ar ’ 
+' PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS './. 
at 
it 


/- 
' PRESS <ENTER> TO CONTINUE OX, 
READ(* , 7002 )ANSWER 
CALL GES 
WRITE(*, 7006) 
FORMAT(//////,' rik OVERVIEW wriikddeek ei 
" B. STRUCTURE PHYSICAL CHARACTERISTICS ( CONTINUED ) eg 


—=_-_ == 


3. AFTER SELECTING THE DESIRED NUMBER OF COPPER aad 
LAYERS AND NODAL ASPECT RATIO, THE PROGRAM WILL ‘yl 
THEN ASK THE USER TO PROVIDE LAYER CHARACTERISTICS. ',/, 
EPOXY AND COPPER LENGTH, THICKNESS, AND WIDTH, -) 
AS WELL AS CONSTANTS CONSTITUTE THE QUERIES. a 
THE PROGRAM WILL ALSO ASK THE USER TO PROVIDE ',/, 
THE PERCENT COVERAGE FOR EACH COPPER LAYER. il 


PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS ale 


—=-_- ss «5 


+++ +tt+eettttt 


—a 


PRESS <ENTER> TO CONTINUE 2k 
READ(* , 7002) ANSWER 

CALL CLS 

WRITEC* , 7007 ) 

FORMAT(//////,' KKK OVERVIEW wire Fe 
C. INTITIAL AND AMBIENT TEMPERATURES 


1 
b 
1. AFTER SELECTING THE DESIRED NUMBER OF COPPER 7s 
LAYERS AND NODAL ASPECT RATIO, THE PROGRAM WILL - 
THEN ASK THE USER TO PROVIDE TEMPERATURES FOR THE .' 
BOARD. 


—_-_ = s «= = 
-~-= 


PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS 


<a 


t+t+e ee teteetett 


PRESS <ENTER> TO CONTINUE 
READ(* , 7002) ANSWER 

CALL CLS 

WRITE(*, 7008) 

FORMAT(//////,' dere OVERVIEW wiciieiddidk Ye 
D. HEAT INEGR. ie 


—_-_ = « 
=a 


1. HEAT INJECTION OCCURS ONLY ON THE UPPER COPPER fe 
LAYER. THIS PROGRAM SUPPLIES THE USER FOR ‘y/ 
; 


-=-_— = 


ALTERNATIVE METHODS FOR ENTERING HEAT. 


A.) TOTAL HEAT OVER SURFACE. 
B. ) AVERAGE HEAT PER UNIT AREA ' 
C.) INPUT HEAT NODE BY NODE ugh. 
D.) NO HEAT INPUT | 


- «= 


PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS we 


Ww 
~ 
ww 


+++ teteeeeeteteretes+ 


— 


****THIS CONCLUDES THE PROGRAM OVERVIEW**** Se 
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+' PRESS <ENTER> TO CONTINUE mex, ) 
READ(* , 7002) ANSWER 
CALL CLS 
END 


SUBROUTINE PCBS1 
Rr Tetetedeveds RPK ERT HR RRR RTE KR Re 


ar cE: MODEL BUILDER 

SUBROUT: PCBS1 

AUTHOR: LT STEVE GLASER 

DATE: 09 MAY 1991 

COMPILER: MICROSOFT VERSION 4.01 
LINKER: MICROSOFT VERSION 3.55 


THIS SUBROUTINE HANDLES THE TEMPERATURE INPUTS TO THE PRINTED CIRCUIT 
BOARD. 


COMMON IBT,UPRT,LWRT,RT,LT,FT,BT,UT, ANSWER 


DEFINE REAL VARIABLES 


ReAGel BT SUPRT, LWRT,RI,LI,FT,BT 


DEFINE ONE CHARACTER VARIABLES 


aaqn maaan MOI aaaaN ao Giera 


CHARACTER*1 ANSTEM,ANSTEB , ANSTEU , ANSTEL, ANSTER, ANSLFT , ANSFRT , ANSBC 


on, UT 
C 
31180 CALL CLS 
WRITE(*,60000) 


60000 FORMAT ( ///, | exe tdeederete eK AR AR TERR RRR EE EERE RRR ER IRR RRR 
PIM HAAKAEREKEEERIE | S| ; 
| reikdededikddstai ike AMBIENT TEMPERATURE INPUTS tcccccccdcdcdvdedcdcdedededcdevedededs 
hivekdeddesek ! : 
fe! dedededededededeaede tease de Het de IRR RT TE RRERTIKRR RR EREREEERERERRERI 
irik! ff.) 


C 
@ 
WRITE(* ,60001) UT 
60001 FORMAT(/,' ENTER THE INITIAL BOARD TEMPERATURE (',Al,'): ',2X, 
+) 
READ *,IBT 
C 
@ 
WRITE(*,60002) UT 
60002 FORMAT(/,' ENTER THE UPPER SURFACE AMBIENT TEMPERATURE (',Al,') 
ee eX; +) 
READ *,UPRT 
C 
6 
C 
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WRITE(*,60003) UT 


60003 FORMAT(/, ' ENTER THE LOWER SURFACE AMBIENT TEMPERATURE (',Al,') 
+: 2 
READ *,LWRT 
Ee 
C 
WRITE(*,60004) UT 
60004 FORMAT( /,' ENTER THE RIGHT SURFACE AMBIENT TEMPERATURE (',Al1,') 
Te.) 
READ *,RT 
C 
C 
C 
WRITE(*,60005) UT 
60005 FORMAT(/, ' ENTER THE LEFT SURFACE AMBIENT TEMPERATURE (',Al,'): 
+ 3 3 
READ *,LT 
C 
G 
WRITE(* ,60006) UT 
60006 FORMAT(/ , ' ENTER THE FRONT SURFACE AMBIENT TEMPERATURE (',Al,') 
Tc) 
READ *,FT 
C 
C 
C 
WRITE(*,60007) UT 
60007 FORMAT(/,’ ENTER THE REAR SURFACE AMBIENT TEMPERATURE (',Al,'): 
2K) 
READ *,BT 


C 
C 
C REVIEW THE TEMPERATURE ENTRIES 
C 
6 


0052 CALL CLS 
WRITE(* ,60008) IBT,UT 


60008 FORMAT(////,° YOU HAVE MADE THE FOLLOWING AMBIENT TEMPERATURE E 
+NTRIES: ',////,' 1. INITIAL BOARD TEMPERATURE: ‘',1X,F9.3,1X,A1) 
WRITE(* ,60009) UPRT,UT 

60009 FORMATC(' 2. UPPER AMBIENT TEMPERATURE: ',1X,F9.3,1X,Al) 
WRITE(*,60010) LWRT,UT 

60010 FORMAT( ' 3. LOWER AMBIENT TEMPERATURE: ',1X,F9.3,1X,A1) 
WRITE(* ,60011) RT,UT 

60011 FORMATC( ' 4. RIGHT AMBIENT TEMPERATURE: ',1X,F9.3,1X,A1) 
WRITE(*,60012) LT,UT 

60012 FORMAT( ' 5. LEFT AMBIENT TEMPERATURE: ',1X,F9.3,1X,A1) 
WRITE(* ,60013) FT,UT 

60013 FORMATC( ' 6. FRONT AMBIENT TEMPERATURE: ',1X,F9.3,1X,A1) 
WRITE(*,60014) BT,UT 

60014 FORMAT(' 7. REAR AMBIENT TEMPERATURE: ',1X,F9.3,1X,A1,//) 

C 

G 

C 
WRITE(* ,60015) 

60015 FORMAT(/,' DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR YES A 


+ND N FOR NO: ',2X, ) 


READ(*, 60016) ANSTEM 
60016 FORMAT(A1) 


MAKE CORRECTIONS BOARD TEMPERATURES 


C2 @ c) 2 @ 


IF(ANSTEM. EQ. 'Y') THEN 
60022 CALL CLS 
WRITE(*,60017) IBT,UT 


60017. FORMAT(///,' THE INITIAL BOARD TEMPERATURE IS: ',F9.3,1X,A1) 
WRITE(*, 60018) 
60018 FORMAT(/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 
+X 


ae) 
READ(*,60019) ANSTEB 
60019 FORMAT(A1) 
IF(ANSTEB. EQ. 'Y') THEN 
WRITE(* ,60020) UT 
60020 FORMAT(/,' ENTER THE NEW VALUE (',Al1,'): ',2X, ) 
READ *,IBT 
ELSEIF(ANSTEB. EQ. 'N') THEN 
GOTO 60021 
ELSE 
GOTO 60022 
ENDIF 
C 
60021 CALLAGIS 
WRITE(*,70022) UPRT,UT 


70022 FORMAT(///,' THE UPPER AMBIENT TEMPERATURE IS: ',F9.3,1X,A1) 
WRITE(* , 60023) 
60023. FORMAT(/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 
+X 


$ 
READ(*,60024) ANSTEU ° 
60024 FORMAT(A1) 
IF(ANSTEU. EQ. 'Y') THEN 

WRITE(*,60025) UT 

60025 FORMAT(/,' ENTER THE NEW VALUE (',Al,'): ',2X, ) 
READ *,UPRT 

ELSEIF(ANSTEU. EQ. 'N') THEN 


GOTO 60026 
ELSE 
GOTO 60021 
ENDIF 
C 
C 


60026 CALL CLS 
WRITE(* ,60027) LWRT,UT 


60027. FORMAT(///,' THE LOWER AMBIENT TEMPERATURE IS: ',F9.3,1X,A1) 
WRITE(* , 60028) 
60028 FORMAT(/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 
+X 


READ(*,60029) ANSTEL 
60029 FORMAT(A1) 
IF(ANSTEL. EQ. 'Y') THEN 
WRITE(* ,60030) UT 
60030 FORMAT(/,' ENTER THE NEW VALUE (',Al,'): ',2X, ) 
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READ *,LWRT 
ELSEIF(ANSTEL. EQ. 'N') THEN 
GOTO 60031 
ELSE 
GOTO 60026 
ENDIF 


G 

C 

6 

60031 CALLLCLS 
WRITE(* ,60032) RT,UT 

60032  FORMAT(///,' THE RIGHT AMBIENT TEMPERATURE IS: ',F9.3,1X,A1) 
WRITE(* , 60033) 

60033. FORMAT(/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 

+X 


—) 
READ(*,60034) ANSTER 
60034  FORMAT(A1) 
IF(ANSTER. EQ. 'Y') THEN 
WRITE(*,60035) UT 
60035 FORMAT(/,' ENTER THE NEW VALUE (',Al,'): ',2X, ) 
READ *,RT 
ELSEIF( ANSTER. EQ. 'N') THEN 
GOTO 60036 
ELSE 
GOTO 60031 
ENDIF 


© 
C 
e 
c 
60036 CALL CES 
WRITE(* ,60037) LT,UT 
60037 FORMAT(///,' THE LEFT AMBIENT TEMPERATURE IS: ',F9.3,1X,A1l) 
WRITE(* ,60038) 
60038  FORMAT(/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 
+X, 
READ(* ,60039) ANSLFT 
60039 FORMAT(A1) 
IF( ANSLFT. EQ. 'Y') THEN 
WRITE(*,60040) UT 
60040 FORMAT(/,' ENTER THE NEW VALUE (',A1,'): ',2X, ) 
READ *,LT 
ELSEIF(ANSLFT. EQ. 'N') THEN 
GOTO 60041 
ELSE 


GOTO 60036 
ENDIF 


0041 CALL CLS 
WRITE(*,60042) FT,UT 

60042  FORMAT(///,' THE FRONT AMBIENT TEMPERATURE IS: ',F9.3,1X,A1) 
WRITE(*, 60043) 

60043 FORMAT(/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 
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+X, ) 
READ(*,60044) ANSFRT 
60044 FORMAT(A1) 
IF(ANSFRT. EQ. 'Y') THEN 
WRITE(* ,60045) UT 
60045 FORMAT(/,' ENTER THE NEW VALUE (',Al,'): ',2X, ) 
READ *,FT 
ELSEIF(ANSFRT. EQ. 'N') THEN 
GOTO 60046 
ELSE 
GOTO 60041 
ENDIF 


maaan 


60046 Ann, CLS 
WRITE(* ,60047) BT,UT 


60047. FORMAT(///,' THE REAR AMBIENT TEMPERATURE IS: ‘,F9.3,1X,A1) 
WRITE(*, 60048 ) 
60048 FORMAT(/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 
+X 


READ(* ,60049) ANSBCK 
60049 FORMAT(A1) 
IF(ANSBCK. EQ. 'Y') THEN 
WRITE(* ,60050) UT 
60050 FORMAT(/,' ENTER THE NEW VALUE (',Al1,'): ',2X, ) 
READ *,BT 
ELSEIF( ANSBCK. EQ. 'N') THEN 
CONTINUE 
ELSE 
GOTO 60046 
ENDIF 
ELSEIF(ANSTEM. EQ. 'N') THEN 
GOTO 60051 
ELSE 
GOTO 60052 
ENDIF 
60051 END 
SUBROUTINE PCBS2( THEAT, THPN, AHEAT , NHEAT ,NPL,NWIDE,NDEEP ,HEAT, IH, JH 
+,ANSN,EL, EW) 


TITLE: MODEL BUILDER 

SUBROUT: PCBS2 -- HEAT INPUT SUBROUTINE 
AUTHOR: LT STEVE GLASER 

DATE: OOSMAY 1991 

MICROSOFT VERSION 4.01 

LINKER: MICROSOFT VERSION 3.55 


THIS SUBROUTINE HANDLES THE HEAT INPUTS TO THE PRINTED CIRCUIT BOARD 


DEFINE REAL VARIABLES 


IAI’q_*~aanagaqaaaqaradanriany 
© 
O 
<< 
rd 
— ° 
on 
try 
Se 


REAL THEAT,THPN, AHEAT ,NHEAT, SL, SW 


7) 
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DEFINE INTEGERS 


©). @) o) 


INTEGER NWIDE ,NDEEP ,NPL,NUM,NUMA,CH,H,TOTNOD,NC,NN,I,J,1M 


DEFINE ONE CHARACTER VARIABLES 


CD esp ay Leg 


CHARACTER*1 SELH,ANSH,ANSN,ATH, ANSHA, AHN , DUMMY 
DEFINE CHARACTER VARIABLES 

CHARACTER UH*13,UAH*6 
DEFINE REAL MATRICES 

REAL HEAT( 100,100) ,COEFF(740,9) 
DEFINE INTEGER VARIABLES 

INTEGER IH( 1000) ,JH( 1000) , UMMY( 1000) 
PROVIDE A CORRELATION BETWEEN NODE NUMBERS AND MATRIX LOCATION 


mWaQanqan CasGac G2 C2 C2 Ga © C21 @& 


1 NUM=1 
DO 60 I=1,NPL/NWIDE 
DO 61 J=1,NWIDE 
JH(NUM)=J 
NUM=NUM+1 
61 CONTINUE 
60 CONTINUE 


NUMA=1 
CH=0 
H=1 
DO 62 I=1,NPL 
IH(NUMA)=H 
CH=CH+1 
IF(CH. EQ. NWIDE) THEN 
H=H+1 
CH=0 
ELSE 
CONTINUE 
ENDIF 
NUMA=NUMA+1 
62 CONTINUE 


C 
1241 CALL CLS 


WRITE(* , 8000) 
8000 FORMAT ( ///, | eiiaiaactae aa t A RRR RANA AN AIR EER RE ERR AAN KARR RETIN 


PIAKKIIKKAAIKKEARKRI | / 
+ debts aaa aaK aaa Ka we KAKKKEEHEAT INPUTS tere ivaedcacacieddeaeaeae aeaetedrdede aed ier 
PAIK! 
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HHA KAAAKRKEKKKKKKKEKAKEKEKERERKEEK REE RERKE REE ERERERERKEKKE 
eae wie 


HEAT INPUT TO THE PCB OCCURS ONLY ON THE UPPER COPPER',/, 
LAYER. HEAT INPUT IS ACCOMPLISHED BY ONE OF THE ',/, 
FOLLOWING METHODS: ',//, 

1. ENTER AS TOTAL HEAT APPLIED TO THE PCB',/, 

2. ENTER AS AVERAGE HEAT PER UNIT AREA ',/, 

3. ENTER HEAT NODE BY NODE ',/, 

4. NO HEAT INPUT',//, 

PLEASE SELECT NUMBER 1 THROUGH 4: ',2X, ) 


++tteeett 


READ(*,5799) SELH 
5799 FORMAT(A1) 

IF(SELH. EQ. '1'. OR. SELH. EQ. '2'. OR. SELH. EQ. '3'. OR. SELH. EQ. '4') THEN 
444  WRITE(*,291) SELH 


291 FORMAT(//,' YOU HAVE SELECTED NUMER ',Al,' OF FOUR ALTERNATIVES 
1, 1S THIS THE, DESIRED SELECTMON@(Y OR N)2 42x, ) 
ELSE 
GOTO 124] 
ENDIF 


READ(* , 5409) ANSH 
5409 FORMAT(A1) 


IF( ANSH. EQ. 'Y') THEN 
5441 GOTO 63 
ELSEIF(ANSH. EQ. 'N') THEN 
GOTO 51 
ELSE 
CALL CLS 
WRITE(* , 8000) 
GOTO 444 
3 ENDIF 


DETERMINE UNIT FOR HEAT INPUT 


aad a a 


IF ((ANSN. EQ. jE ANDeCSELH. EG. 2° )i) THEN 
UH = ‘Btu/(hr*in 2 
UAH = 'Btu/hr' 

ELSEIF ((ANSN. EQ. 'E'). AND. (SELH. EQ. '1'. OR. SELH. EQ. '3')) THEN 
UH = 'Btu/hr' 

ELSEIF ((ANSN. EQ. 's'). AND. (SELH. EQ. '2' )) THEN 
UH = ‘Watts/(cm ip 
UAH = 'Watts' 

ELSEIF ((ANSN. EQ. 'S'). AND. (SELH. EQ. '1'. OR. SELH. EQ. '3')) THEN 
UH = ‘Watts’ 


ALLOW FOR RE-SELECTION OF HEAT INPUT METHOD OR CONTINUE WITH 
INITIAL SELECTION. 


CHOICE #1 


aqaaadagaiqa 
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989 
988 


996 


coco C2 C2 C2 


998 
997 


*C 


CALEOCLS 
IF (SELH. EQ. '1'. AND. ANSN. EQ. 'E’) THEN 
WRITE(*,988) UH 


FORMAT(///,' YOU HAVE SELECTED TO INPUT HEAT AS A TOTAL HEAT 
ae es APPLIED TO THE SURFACE. ',//, 
+" ENTER TOTAL HEAT APPLIED TO THE SURFACE (',A6,'): ',2X, ) 


READ *,THEAT 

WRITE(*,991) 

FORMAT(/,' IS THIS THE CORRECT ENTRY? (Y OR N): ‘',2X, ) 
READ(* ,992)ATH 

FORMAT( Al) 


MAKE ENTRY AND ALLOW FOR CORRECTIONS 
if (CATH, BOs YeetHEN 


THPN=THEAT/NPL 
WRITE(* 993) THPN WH 


FORMAT(/,' TOTAL HEAT PER NODE IS: ',F9.4,1X,A13) 
WRITE(* , 9323) 

FORMAT(//,' <PRESS ENTER TO CONTINUE>' ) 
READ(* , 9324 )DUMMY 

FORMAT(A1) 


FILL HEAT MATRIX WITH DESIRED VALUES 


DO 994 I=1,NDEEP 
DO 995 J=1,NWIDE 
HEAT(I,J)=THPN 
CONTINUE 
CONTINUE 


ELSEIF( ATH. EQ. 'N') THEN 
CALGACLS 
GOTO 989 

ELSE 
CALL CLS 
WRITE(*,996) THEAT,UH 
FORMAT(////,' TOTAL HEAT APPLIED TO THE SURFACE IS: ',F9 

+.4,1X,A13, ) 
GOTO 990 
ENDIF 


CHOICE #2 


ELSEIF(SELH. EQ. '2') THEN 
WRITE(*,997) UH 


FORMAT(////,° YOU HAVE SELECTED TO ENTER THE AVERAGE HEA 
+T OVER THE',/,' UPPER PCB SURFACE.',//,' ENTER THE DESIR 
+ED HEAT INPUT: @ ,AlSy) 2) 32) 

READ *,AHEAT 
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G 

C 
1000 
he 


1001 


1006 


Cy C2 C17 C3 


1200 
1201 


Ioaaaaqaraiannay 


H202 


MAKE ENTRY AND ALLOW FOR CORRECTION 


WRITE( *,999 ) 

FORMAT(/,' Secs THE CORRECT ENTRY? a0 Y OR Njageoe, ) 
READ(*, 1001) ANSHA 

FORMAT(A1) 


IF (ANSHA. EQ. 'Y') THEN 
THPN=AHEAT*EL*EW /NPL 
WRITE(* ,1002) THPN,UAH 


FORMAT(///,' TOTAL HEAT PER NODE IS: ',F9.4,1X,A6) 

WRITE(* ,4323) 

FORMAT(//,' <PRESS ENTER TO CONTINUE>' ) 
READ(* ,4324) DUMMY 

FORMAT(A1) 


FILL HEAT MATRIX WITH DESIRED VALUES 


DO 1003 I=1,NDEEP 
DO 1004 J=1,NWIDE 
HEAT(I,J)=THPN 
CONTINUE 
CONTINUE 


ELSEIF(ANSHA. EQ. 'N') THEN 
CALL CLS 
GOTO 998 

ELSE 
GALE CLS 
WRITE(*,1006) AHEAT,UH 
FORMAT(////,' AVERAGE HEAT OVER PCB SURFACE IS: ',F9.4,1 

+X,A13, ) 
GOTO 1000 
ENDIF 


CHOICE #3 


ELSEIF(SELH. EQ. '3') THEN 


WRITE(*, 1201) 

FORMAT(////,' YOU HAVE SELECTED TO ENTER THE HEAT NODALLY' 
t,//,5" ENTER THE TOTAL NUMBER OF NODES DESIGNATED FOR HEAT INP 
mor: 6 5 2k) 


READ * ,TOTNOD 


THIS IS DONE NODE BY NODE. GET NUMBER OF ENTRIES AND THEN LOOP UNTIL 
ALL ENTRIES HAVE BEEN MADE. 
TELL USER MAXIMUM NUMBER OF ENTRIES POSSIBLE 

IF (TOTNOD. GT. NPL) THEN 


WRITE(*,1202) NPL 
FORMAT(///,' THE MAXIMUM ENTRY IS: ',I4) 


al 


2239 


2240 


Go ccf 


1204 


1205 


5345 


8032 


i Oa 
1302 
11204 
11203 


1303 


1304 


1305 


WRITEG*| 2239) 


FORMAT(///,' PLEASE PRESS <ENTER> TO CONTINUE a 
READ(*,2240) DUMMY 
FORMAT( Al) 
GALL CLS 
GOTO 1200 
ENDIF 
MAKE ENTRIES 
DO 1203 I=1,TOTNOD 
NC=I 
CALL CLS 
WRITE(* , 1205) NC,TOTNOD 
FORMAT( ///,' THIS IS NUMBER ',I3,' OF ',13,' ENTRIES') 
WRITE(* ,5345) 
FORMAT(//,' ENTER THE NODE NUMBER FOR HEAT INPUT: ',2 
hae) 
READ (*,8032) NN 
FORMAT( 14) 
UMMY( I)=NN 
IF (NN. EQ. 0. OR. NN. GT. NPL) THEN 
GOTO 1204 
ENDIF 
WRITE(*, 1206) UH 
FORMAT(/,' ENTER THE HEAT INPUT (',A6,'): ',2X, ) 
READ * ,NHEAT 
HEAT(IH(NN) , JH(NN)) = NHEAT 
CONTINUE 
PROVIDE OPPORTUNITY TO MAKE CORRECTIONS 


CALEBGIS 
WRITE(*,1302) TOTNOD 
FORMAT(///,' YOU HAVE MADE ',I3,' NODAL ENTRIES: ',/) 


DO 11203 IM=1,TOTNOD 
WRITEC* ,11204) UMMYC IM) ,HEAT( IHC UMMY( IM) ) | JHCUMMY Gite 


FORMAT(1X,' NODE NUMBER ',14,': ',1X,F9.4,1X,A6) 

CONTINUE 

WRITE(*, 1303) 

FORMAT(/,' DO YOU WISH TO MAKE ANY MORE ENTRIES OR CORRECTION 


+S? (Y OR N): ',2X, ) 


READ(*,1304) AHN 
FORMAT(A1) 
IF( AHN. EQ. 'Y') THEN 
CALI ells 
GOTO 1200 
ELSEIF( AHN. EQ. 'N') THEN 
GOTO 1305 
ELSE 
GOTO 1301 
ENDIF 
ELSEIF(SELH. EQ. '4') THEN 
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GOTO 1306 
1306 ENDIF 
END 
SUBROUTINE PCBS3( XEW, XEL, XCW,XCL, XUKE, XDELE, XUKC ,SELECT,X1L1,X2L1, 
+X2L2,43L1,43L2, X3L3 , X4L1 , X4L2,X4L3 ,X4L4 , XWIDE , XDEEP, XPL, XIH, XJH,X 
fil, KOBE, AIBI ,XPRI,XWRT, XRT ,XLT,XFT,XBT,XNITS) 


CRAKEAR KEHK KKK KKK KKK KK KKK KE KEE ERIK EKEKRERKEEKEKERERERE 


C 
oe «TITLE: MODEL BUILDER 
C SUBROUTINE: PCBS3 
C AUTHOR: LT STEVE GLASER 
C ‘DATE: 09 JUL 1991 
C COMPILER: MICROSOFT VERSION 4.01 
C ‘LINKER: MICROSOFT VERSION 3.55 
C 
C THIS SUBROUTINE IS THE "HEART" OF THE MODEL BUILDER. THE 
c NODE COEFFICIENTS ARE CALCULATED IN THIS SUBROUTINE AND 
C PLACED IN MATRICES. PCBS3 IN TURN ALSO CALLS OTHER 
C SUBROUTINES: S1, S2, S3, AND S4. THESE SUBROUTINES CALLED 
C BY PCBS3 GENERATE THE OUTPUT DATA FILE FOR THE THERMAL 
C ANALYZER. 
C 
C 
% 
C DEFINE REAL VARIABLES 
C 
REAL XEL,XEW,XDELE, XUKE ,XCL,XCW,X1L1,X2L1,X2L2,X3L1,X3L2,X3L3,X4L1 
+,X4L2,X4L3,X4L4 , XUKC,EZB 
REAL EYLR,EFXB,EYY,EZC11,£2C21,E2022,£2C31,E2C32,£2C33,EZC41,E2C42 
+ ,EZC43 ,EZC44 , EXX , DELX ,DELY 
REAL CYLR11,CYLR21,CYLR22 ,CYLR31,CYLR32,CYLR33,CYLR41,CYLR42 , CYLR4 
+3, CYLR44 
REAL CFXB11,CFXB21,CFXB22,CFXB31,CFXB32,CFXB33,CFXB41,CFXB42 ,CFXB4 
+3, CFXB44 
REAL CYY11,CYY21,CYY22 ,CYY31,CYY32 ,CYY33,CYY41,CYY42,CYY43,CYY44 
REAL CXX11,CXX21,CXX22 ,CXX31,CXX32 ,CXX33 , CXX41 , CXX42 ,CXX43, CXX44 
REAL CZE11,CZE21,CZE22 ,CZE31,CZE32 ,CZE33 ,CZE41,CZE42 ,CZE43 ,CZE44 
REAL PC11,PC21,PC22,PC31,PC32,PC33,PC41,PC42,PC43,PC44 
REAL XEAT( 100,100) ,XOEF(1000,12) 
C 
C DEFINE CHARACTER VARIABLES 
€ 
C 
CHARACTER*1 SELECT,ANS,XNITS 
C 
C DEFINE INTEGER VARIABLES 
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INTEGER XAWIDE, XDEEP,XPL,NUMC,I,N,IB,USEL 


INTEGER XIH( 1000) ,XJHC 1000) 
C 
CIR TIT TTR TTC TTT FREER TREE TEER EE ERR ERER ER ERREEREEERE REREERERERE 
Crete sdek eek kkk cee DETERMINE COEFFICIENTS ¥ i vx vk eavcedstecek tcteksernak ick iv 
CREE EEE ERE REERERERERRERE ERR ERR AUERRRERERRERREREREREERREEEEEREERE 


C 


C 
C DETERMINE INCREMENTAL MEASUREMENTS IN THE X AND Y DIRECTIONS 
C 
C 
C XNITS IS A CODE USED IN THE OUTPUT DATA FILE. 1 MEANS ENGLISH UNITS, 
C AND 2 MEANS SI UNITS 
C 
C 
IF (XNITS. EQ. 'E') THEN 
USEL=1 
ELSE 
USEL=2 
ENDIF 
C 
C DETERMINE INCREMENTS IN THE ARRAY 
C 


DELX = XEL/XDEEP 
DELY = XEW/XWIDE 
111 CONTINUE 


1h 
C 
C 
C COPPER COVERAGE, ONE COPPER LAYER CASE 
C 
C 


IFCSELECT, EQ. >) THEN 
931 CALL CLS 


ENTER PERCENT COPPER COVERAGE 


WRITE(*, 3941) 
3941 FORMAT(////,' ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR 

+ LAYER 1: ',2X, 

READ *,PC11 

IF (PC11.GT. 100) THEN 

Goer 931 

ENDIF 

PC11=PC11/100 


COPPER COVERAGE, TWO COPPER LAYERS CASE 


Gece co CC) C2 Co C 
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ELSEIF( SELECT. EQ. '2') THEN 
e332 CALL CLS 
WRITE(*, 3942) 
3942 FORMAT(////,/ ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR 
+ LAYER 1: ',2X, 
READ *,PC21 
IF (PC21.GT.100) THEN 
GOTO 932 
ENDIF 
E@21=PC21/100 
933. WRITE(*,3943) 
3943 FORMAT(/,' ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR LA 
+YER 2: ',2X, ) 
READ *,PC22 
IF (PC22.GT.100) THEN 
GOTO 933 
ENDIF 
PC22=PC22/100 


COPPER COVERAGE, THREE COPPER LAYERS CASE 


qNanran 


ELSEIF( SELECT. EQ. '3') THEN 
eo) CALL CLS 
WRITE(* , 4942) 
4942 FORMAT(////,' ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR 
+ LAYER 1: ',2X, 
READ *,PC31 
IF (PC31.GT. 100) THEN 
GOTO 942 
ENDIF 
PC31=PC31/100 
433 WRITE(*,4943) 
4943 FORMAT(/,' ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR LA 
+YER 2: ',2X, ) 
READ *,PC32 
IF (PC32.GT.100) THEN 
GOTO 433 
ENDIF 
PC32=PC32/100 
633 WRITE(*,8943) 
8943 FORMAT(/,' ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR LA 
+YER 3: ',2X, ) 
READ *,PC33 
IF (PC33.GT. 100) THEN 
GOTO 633 
ENDIF 
PC33=PC33/100 


COPPER COVERAGE, FOUR COPPER LAYERS CASE 


C2 C7 <2 2 


ELSEIF( SELECT. EQ. '4') THEN 
Ze2)) CALL CLS 
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WRITE(* , 2942) 
2942 FORMAT(////,' ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR 
+ LAYER 1: ',2X, ) 
READ *,PC41 
IF (PC41.GT. 100) THEN 
GOTO 242 
ENDIF 
PC41=PC41/100 
233. WRITE(*, 2943) 
2943 FORMAT(/,' ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR LA 
+YER 2: ',2X, ) 
READ *,PC42 
IF (PC42.GT. 100) THEN 
GOTO 233 
ENDIF 
PC42=PC42/100 
133. WRITE(*,1943) 
1943 FORMAT(/,' ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR LA 
+YER 3: ',2X, ) 
READ *,PC43 
IF (PC43.GT.100) THEN 
GOTO 133 
ENDIF 
PC43=PC43/100 
Be WRITE(* ,43) 
43 FORMAT(/,' ENTER THE PERCENT COPPER COVERAGE (i.e. 67) FOR LA 
+YER 4: ',2X, ) 
READ *,PC44 
IF (PC44.GT.100) THEN 
GOTO 33 
ENDIF 
PC44=PC44/100 
ENDIF 
CALLaCES 


VERIFY THAT THE SELECTED ENTRIES ARE CORRECT 


Cac) C) Ca Ca C2 Ca 


IF( SELECT. EQ.'1') THEN 
305 WRITE(*,300) PC11*100 
300 FORMAT(/////' YOU HAVE SELECTED ',F6.2,'% COVERAGE.',//,' 
+ IS THIS YOUR DESIRED ENTRY? (Y OR N): ',1X, 
READ (*,301)ANS 
301 #FORMAT(A1) 
IF( ANS. EQ. 'Y') THEN 
GOTO 333 
ELSEIF( ANS. EQ. 'N') THEN 
Geno iis 
ELSE 
CALL AGES 
GOTO 305 
ENDIF 
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306 
307 


308 


310 
Shall 


312 


520 
321 


B22 


W 
W 


qQaqraaqawiaanaqnad Ww 


ELSEIF( SELECT. EQ. '2') THEN 

WRITE(* ,307) PC21*100,PC22*100 

FORMAT( /////' YOU HAVE SELECTED: ',F6.2,'% COVERAGE, LAYER 1' 
ta 25%,F6.2, % COVERAGE, LAYER 2.',//,. IS THIS YOUR DESIRED 
+ ENTRY? (Y ORN): ',1X, ) 

READ (*,308)ANS 

FORMAT(A1) 

IF( ANS. EQ. 'Y') THEN 

GOTO 333 

ELSEIF( ANS. EQ. 'N' ) THEN 

GOTO 1111 

ELSE 

GALL CLS 

GOTO 306 

ENDIF 

ELSEIF( SELECT. EQ. '3') THEN 

WRITE(*,311) PC31*100,PC32*100,PC33*100 

FORMAT( /////' YOU HAVE SELECTED: ',F6.2,'% COVERAGE, LAYER 1' 
+,/,25X,F6.2,'% COVERAGE, LAYER 2',/,25X,F6.2,'% COVERAGE, LAYER 3. 
ee | IS THIS YOUR DESIRED ENTRY? (Y OR N): ',1X, ) 

READ (*,312)ANS 

FORMAT(A1) 

IF(ANS. EQ. 'Y') THEN 

GOTO 333 

ELSEIF(ANS. EQ. 'N') THEN 

COTO 1111 

ELSE 

CALL CLS 

GOTO 310 

ENDIF 

ELSEIF( SELECT. EQ. '4') THEN 

WRITE(*,321) PC41*100,PC42*100, PC43*100 , PC44*100 


FORMAT(/////' YOU HAVE SELECTED: ',F6.2,'% COVERAGE, LAYER 1' 
+,/,25X,F6.2,'% COVERAGE, LAYER 2',/,25X,F6.2,'% COVERAGE, LAYER 3° 
+,/,25X,F6.2,'% COVERAGE, LAYER 4.',//,' IS THIS YOUR DESIRED 


+ ENTRY? (Y OR N): ',1X, ) 
READ (*,322)ANS 
FORMAT(A1) 
IF(ANS. EQ. 'Y') THEN 
GOTO 333 
ELSEIF( ANS. EQ. 'N') THEN 
GOTO 1111 
ELSE 
PALL CLS 
GOTO 320 
ENDIF 
ENDIF 
CONTINUE 


GENERATE CONSTANTS FOR THE EPOXY LAYERS 


LEFT OF RIGHT EDGE TO OUTSIDE 


77 


SP ip Sp Vi ee Ua a PE a a Tal Sp Ie >) > A Ge Ea) C2 0) C2 Caio (2° C2) C2 Gea C2 C2) Goa 


c3°G) @) Gi C) C20) 


Circa ca, C2 C2 C2 C2 C2 


EYLR = 2*XUKE*DELX*XDELE/DELY 
FRONT OR BACK TO OUTSIDE 


EFXB = 2*XUKE*DELY*XDELE /DELX 

INNER MATRIX MOVEMENT IN THE Y DIRECTION 
EYY = XUKE*DELX*XDELE/DELY 

INNER MATRIX MOVEMENT IN THE X DIRECTION 
EXX = XUKE*XDELE*DELY/DELX 


EPOXY TORGOPE:R 


ONE COPPER LAYER CASE 


IF(SELECT. EQ. '1') THEN 
EZC11= 2*DELX*DELY/((XDELE/XUKE )+(X1L1*PC11/XUKC)) 


TWO COPPER LAYERS CASE 


ELSEIF( SELECT. EQ. '2') THEN 
EZC21 = 2*DELX*DELY/((XDELE/XUKE )+(X2L1*PC21/XUKC) ) 
EZC22 = 2*DELX*DELY/((XDELE/XUKE)+( X2L2*PC22/XUKC)) 


THREE COPPER LAYERS CASE 


ELSEIF(SELECT. EQ. '3') THEN 
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EZC31 = 2*DELX*DELY/((XDELE/XUKE)+(X3L1*PC31/XUKC)) 
EZC32 = 2*DELX*DELY/((XDELE/XUKE )+(X3L2*PC32/XUKC) ) 
EZC33 = 2*DELX*DELY/((XDELE/XUKE )+( X3L3*PC33/XUKC) ) 

c 

C 

C FOUR COPPER LAYERS CASE 

C 

® 

c 

c 
ELSEIF(SELECT. EQ. '4') THEN 
EZC41 = 2*DELX*DELY/((XDELE/XUKE )+(X4L1*PC41/XUKC)) 
EZC42 = 2*DELX*DELY/((XDELE/XUKE )+(X4L2*PC42 /XUKC)) 
EZC43 = 2*DELX*DELY/((XDELE/XUKE )+( X4L3*PC43/XUKC)) 
EZC44 = 2*DELX*DELY/((XDELE/XUKE)+( X4L4*PC44/XUKC) ) 
ELSE 
CONTINUE 
ENDIF 

G 

C 

c 

@ EPOXY TO BOTTOM OUTER EDGE 

C 


EZB = 2*XUKE*DELX*DELY/XDELE 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


eG 

C GENERATE CONSTANTS FOR THE COPPER LAYERS 

@ 

€ 

C 

C LEFT OF RIGHT EDGE TO OUTSIDE 

C 

C 

C ONE COPPER LAYER CASE 

C 

€ 

C 

C 

C 
IF(SELECT. EQ. '1') THEN 
CYLR11= 2*XUKC*X1L1*PC11/DELY 

C 

C 

c TWO COPPER LAYERS CASE 

@ 

C 

6 

@ 


ELSEIF( SELECT. EQ. '2') THEN 
CYLR21= 2*XUKC*X2L1*PC21/DELY 
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De hw elem El ik a 


ee CP 


CYLR22= 2*XUKC*X2L2*PC22/DELY 


THREE COPPER LAYERS CASE 


ELSEIF( SELECT. EQ. '3') THEN 

CYLR31= 2*XUKC*X3L1*PC31/DELY 
CYLR32= 2*XUKC*X3L2*PC32/DELY 
CYLR33= 2*XUKC*X3L3*PC33/DELY 


FOUR COPPER LAYERS CASE 


ELSEIF(SELECT. EQ. '4') THEN 

CYLR41= 2*XUKC*X4L1*PC41/DELY 
CYLR42= 2*XUKC*X4L2*PC42/DELY 
CYLR43= 2*XUKC*X4L3*PC43/DELY 
CYLR44= 2*XUKC*X4L4*PC44/DELY 


ELSE 
CONTINUE 
ENDIF 


FRONT OR BACK TO OUTSIDE 


ONE COPPER LAYER CASE 


IF(SELECT. EQ. '1') THEN 
CFXBl1= 2*XUKC*DELY*X1L1*PC11/DELX 


TWO COPPER LAYERS CASE 


ELSEIF( SELECT. EQ. '2') THEN 
CFXB21= 2*XUKC*DELY*X2L1*PC21/DELX 
CFXB22= 2*XUKC*DELY*X2L2*PC22/DELX 


THREE COPPER LAYERS CASE 
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aqadianq 


Bs Sl el ep 1 ae | ap i a 


IqIaqaaanrananyn C2 C2 Conga? OC) Gc) 


Gaec)? ©) ClaGaaearC) 


ELSEIF( SELECT. EQ. '3') THEN 

CFXB31= 2*XUKC*DELY*X3L1*PC31/DELX 
CFXB32= 2*XUKC*DELY*X3L2*PC32 /DELX 
CFXB33= 2*XUKC*DELY*X3L3*PC33/DELX 


FOUR COPPER LAYERS CASE 


ELSEIF( SELECT. EQ. '4') THEN 

CFXB41= 2*XUKC*DELY*X4L1*PC41/DELX 
CFXB42= 2*XUKC*DELY*X4L2*PC42 / DELX 
CFXB43= 2*XUKC*DELY*X4L3*PC43 /DELX 
CFXB44= 2*XUKC*“DELY*X4L4*PC44 /DELX 


ELSE 
CONTINUE 
ENDIF 


INNER MATRIX MOVEMENT IN THE Y DIRECTION 


ONE COPPER LAYER CASE 


IF(SELECT. EQ. '1') THEN 


CYY11= XUKC*DELX*X1L1*PC11/DELY 


TWO COPPER LAYERS CASE 


ELSEIF(SELECT. EQ. '2') THEN 


CYY21= XUKC*DELX*X2L1*PC21/DELY 
CYY22= XUKC*DELX*X2L2*PC22/DELY 


THREE COPPER LAYERS CASE 


ELSEIF( SELECT. EQ. '3') THEN 
CYY31= XUKC*DELX*X3L1*PC31/DELY 
CYY32= XUKC*DELX*X3L2*PC32/DELY 


8] 


a I ap ie ap Ble Je oP hsp le, 


DS Me ep eS al ep ae | CG) CC) CC) eC) Gs Ga 


Ceca Gc C2 C) 


C2 C7 €2 


CYY33= XUKC*DELX*X3L3*PC33 /DELY 


FOUR COPPER LAYERS CASE 


ELSEIF( SELECT. EQ. '4')} THEN 

CYY41= XUKC*DELX*X4L1*PC41/DELY 
CYY42= XUKC*DELX*X4L2*PC42/DELY 
CYY43= XUKC*DELX*X4L3*PC43/DELY 
CYY44= XUKC*DELX*X4L4*PC44/DELY 


ELSE 
CONTINUE 
ENDIF 


INNER MATRIX MOVEMENT IN THE X DIRECTION 


ONE COPPER LAYER CASE 


IF( SELECT. EQ. '1') THEN 


CXX11= XUKC*X1L1*PC11*DELY/DELX 


TWO COPPER LAYERS CASE 


ELSEIF(SELECT. EQ. '2') THEN 


CXX21= XUKC*X2L1*PC21*DELY/DELX 
CXX22= XUKC*X2L2*PC22*DELY /DELX 


THREE COPPER LAYERS CASE 
ELSEIF( SELECT. EQ. '3')} THEN 
CXX31= XUKC*X3L1*PC31*DELY/DELX 


CXX32= XUKC*X3L2*PC32*DELY /DELX 
CXX33= XUKC*X3L3*PC33*DELY/DELX 


FOUR COPPER LAYERS CASE 
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aqaanaacn 
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ELSEIF( SELECT. EQ. '4') THEN 

CXX41= XUKC*X4L1*PC41*DELY/DELX 
CXX42= XUKC*X4L2*PC42*DELY /DELX 
CXX43= XUKC*X4L3*PC43*DELY /DELX 
CXX44= XUKC*X4L4*PC44*DELY /DELX 


ELSE 
CONTINUE 
ENDIF 


COPPER TO EPOXY (OR AIR) 


ONE COPPER LAYER CASE 


IF( SELECT. EQ. '1') THEN 
CZE11= 2*XUKC*DELX*DELY/(X1L1*PC11) 


TWO COPPER LAYERS CASE 


ELSEIF( SELECT. EQ. '2') THEN 


CZE21 2*DELX*DELY/( (XDELE/XUKE )+( X2L1*PC21/XUKC) ) 
CZE22 = 2*XUKC*DELX*DELY/(X2L2*PC22) 


THREE COPPER LAYER CASE 


ELSEIF( SELECT. EQ. '3') THEN 


CZE31 = 2*DELX*DELY/(( XDELE/XUKE )+( X3L1*PC31/XUKC) ) 
CZE32 = 2*DELX*DELY/( (XDELE/XUKE )+( X3L2*PC32/XUKC ) ) 
CZE33 = 2*XUKC*DELX*DELY/(PC33*X3L3) 


FOUR COPPER LAYER CASE 


ELSEIF( SELECT. EQ. '4') THEN 
CZE41 = 2*DELX*DELY/( (XDELE/XUKE)+(X4L1*PC41/XUKC) ) 
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CZE42 = 2*DELX*DELY/((XDELE/XUKE)+(X4L2*PC42/XUKC) ) 
CZE43 = 2*DELX*DELY/((XDELE/XUKE)+( X4L3*PC43/XUKC) ) 
CZE44 = 2*XUKC*DELX*DELY/(PC44*X4L4) 

ELSE 

CONTINUE 

ENDIF 


CALL THE OUTPUT DATA FILE SUBROUTINES 


ONE COPPER LAYER CASE 


IF( SELECT. EQ. '1') THEN 

CALL S1(EYLR,EFXB,EYY,EXX,E2C11,E2B,CYLR11 ,CFXB11,CYY11, Cxx eee 
+1,XPL,XIH,XJH,XEAT, XWIDE , XDEEP, XOEF ,XIBT, XPRT, XWRT, XRT, XLT, XFT, XBT 
+,USEL) 


TWO COPPER LAYER CASE 


ELSEIF(SELECT. EQ. '2') THEN 

CALL S2(EYLR,EFXB,EYY ,EXX,EZC21,EZC22 ,EZB,CYLR21,CYLR22,CFXB2I9eue 
+B22 ,CYY21,CYY22,CXX21,CXX22 ,CZE21,CZE22 ,XPL, XIH, XJH, XEAT, AW IDE 
+EP,XOEF,XIBT,XPRT,XWRT, XRT, XLT, XFT,XBT,USEL) 


THREE COPPER LAYER CASE 


ELSEIF(SELECT. EG. 3°) THEN 

CALL S3(EYLR,EFXB ,EYY ,EXX,E2031,EZC32,EZC33 ,EZB,CYLR31,\CYUR32 ee 
+33 ,CFXB31,CFXB32 , CFXB33 ,CYY31,CYY32,CYY33 , CXXS I Cy 3200 on eee 
+CZE32 ,CZE33 , XPL,XIH,XJH, XEAT, XWIDE, XDEEP, XOEF (XIBI (APRIAWR i An ee 
LT AFT, XBT SUSE 


FOUR COPPER LAYER CASE 


ELSEIF(SELECT. EQ. '4') THEN 
CALL S4(EYLR,EFXB,EYY,EXX,EZC41,E2Z042 E2043, E2044, 628.0) bal, CY en 
+2 ,CYLR43,CYLR44 ,CFXB41,CFXB42 , CFXB43 , CFXB44 , CYY41,CYY42 ,CYY43 ,CYY4 
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+4 ,CXX41 ,CXX42 ,CXX43 , CXX44 ,CZE41 ,CZE42 ,CZE43 ,CZE44 ,XPL, XIH, XJH, XEAT 
fo LOE, ADEEP. XOEF ,XIBT ,XPRI, XWRT, XRT, XR XPT, XBT ,USEL) 

ELSE 

CONTINUE 

ENDIF 

END 

SUBROUTINE COPPER(CEW,EL,ANSN,UK,SELECT,T1L1,T2L1,T2L2,T3L1,T3L2,T3 
+L3 ,T4L1,T4L2 ,T4L3 ,T4L4,SPEVAL,UKC ,CL,CW) 


TITLE: MODEL BUILDER 

AUTHOR: LT STEVE GLASER 

DATE: 09 JUL 1991 

COMPILER: MICROSOFT VERSION 4.01 
LINKER: MICROSOFT VERSION 3.55 


DEFINE REAL VARIABLES 


SUBROUTINE COPPER GATHERS ALL THE INFORMATION REQUIRED FOR THE COPPER 
LAYERS. 

REAL EL,EW,UN,T1L1,T2L1,T2L2,T3L1,T3L2,T3L3,T4L1,T4L2,T4L3 ,T4L4, Vo 
+L1,VOL2,VOL3 , VOL4 ,T1W1,T2W1,T2W2,T3W1,T3W2,T3W3 ,T4W1,T4W2,T4W3 , T4W 
+4 ,UKC 


Gemc 2 C1 C2) C2 CC) Co Co C2 Gaca 


INTEGER NWIDE,NDEEP ,NPL 
C DEFINE ALL ONE CHARACTER VARIABLES 


CHARACTER*1 ANSN, SELECT ,WORL, SPEVAL, ANSE1,ANST1,ANSK1 
C CHARACTER VARIABLES OF MORE THAN ONE POSITION 


CHARACTER LOCO*2,UK*10,UT*1 


C 

C DEFINE MATRICES 

C 

C 

C 

C 

C 

C COPPER LAYER CHARACTERISTICS 
C 

C 


5458 CALL CLS 


31171 WRITE(*, 18000) 
18000 FORMAT ( / // , | tev dedededederede sede rede Fete de de He HH HHH II IIH HIRAI IKK INK IK HEHE HIKER 


PieaeyerewIAkHIHRRK IHW KIK KI | | / 

+! kickin KKKKAKEEKCOPPER LAYER CHARACTERISTICS Xi x irrdesedededered deserter sever 
piik HHH! | / 

4! detesededededededededenededede dedededete KAI I I I CRITI ICI KIRK IHKKIET KKK KI 
pH! |) 


PROVIDE CORRECT UNIT ABBREVIATIONS 


ap Pap el ae 


AREA=EW*EL 
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IF( ANSN. EQ. 'S') THEN 
UN=12*2. 54*12%2. 54*, 00134374*2. 54 


WRITE(*, 15400) 
15400 FORMAT(’ ALL ENTRIES ARE IN SI NOTATION. ',/) 
LOCO=' cm' 
C UK='Watts/cm/C' 
C Ci — ey 


ELSEIF(ANSN. EQ. 'E') THEN 
UN=12*12*0. 00134374 


WRITE(*,15401) 
15401 FORMAT(' ALL ENTRIES ARE IN ENGLISH NOTATION’ ,/) 
1GOCO= nm 
C UX= Btu nee 
c Ui=. 
ENDIF 


Se Be eS 


WRITE(*,15402) EL,LOCO 
15402 FORMAT(/,' THE COPPER LAYER LENGTH IS THE SAME AS THE EPOXY LAYER 
tit EO at XN. 2 Xe 
CleEu 
C 
WRITE(*,15403) EW,LOCO 
15403 FORMAT(//,' THE COPPER LAYER WIDTH IS THE SAME AS THE EPOXY LAYER 
+: ' FO. 4, 1X42 2 
CW=EW 
14337 WRITE(*,16403) SELECT 
16403 FORMAT(//,' YOU SELECTED ',Al,' COPPER LAYER(S) FOR THE PCB.',/,' 
+ YOU NOW WILL BE ASKED TO ENTER THE COPPER LAYER THICKNESS’ ,2X, ) 
99999 WRITE(*, 31112) 
31112 FORMAT(//,' DO YOU WANT TO SPECIFY THICKNESS BY LENGTH OR WEIGHT? 
+',/,' ENTER L IF YOU WANT TO ENTER LENGTH, W IF YOU WANT TO ENTER 
+ WEIGHT: ',2X, 
READ(*,31113) WORL 
31113 FORMAT(A1) 
IFCWORL. EQ. 'L’. ANDUSELECT. BOQ. 1) )) SHER 
CALL “Gigs 
WRITE(*,31114) LOCO 
31114  FORMAT(////,' ENTER THE THICKNESS FOR LAYER 1 (',A2,'): ',2X, ) 
READ *,T1L1 
COE@ns11 7.0 
ELSEIF(WORL. EQ. 'L'. AND. SELECT. EQ. '2' )THEN 
CALL GELS 
WRITE(* , 31115) LOCO 
31115  FORMAT(////,' ENTER THE THICKNESS FOR LAYER 1 (',A2,'): ',2X, ) 
READ *,T2L1 
WRITE(*,31116) LOCO 
31116 FORMAT(/,' ENTER THE THICKNESS FOR LAYER 2 (',A2,'): ',2X, ) 
READ *,T2L2 
GOTO 31170 
ELSEIF(WORL. EQ. 'L'. AND. SELECT. EQ. '3' ) THEN 
CALIACES 
WRITE(*,31117) LOCO 
31117 FORMAT(////,' ENTER THE THICKNESS FOR LAYER 1 (',A2,'): ',2X, ) 
READ *,T3L1 
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WRITE(*,31118) LOCO 
31118 FORMAT(/,' ENTER THE THICKNESS FOR LAYER 2 (',A2,'): ',2X, ) 
READ *,T3L2 
WRITE(*,31119) LOCO 
31119  FORMAT(/,' ENTER THE THICKNESS FOR LAYER 3 (',A2,'): ',2X, ) 
READ *,T3L3 
GOTO 31170 
ELSEIF(WORL. EQ. 'L'. AND. SELECT. EQ. '4' )THEN 
CALL CLS 
WRITE(*,31120) LOCO 
31120 FORMAT(////,' ENTER THE THICKNESS FOR LAYER 1 (',A2,'): ',2X, ) 
READ *,T4L1 
WRITE(*,31121) LOCO 
31121 FORMAT(/,' ENTER THE THICKNESS FOR LAYER 2 (',A2,'): ',2X, ) 
READ *,T4L2 
WRITE(*,31122) LOCO 
31122 FORMAT(/,' ENTER THE THICKNESS FOR LAYER 3 (',A2,'): ',2X, ) 
READ *,T41L3 
WRITE(*,31123) LOCO 
31123. FORMAT(/,' ENTER THE THICKNESS FOR LAYER 4 (',A2,'): ',2X, ) 
READ *,T4L4 
GOTO 31170 
ELSEIF(WORL. EQ. 'W'. AND. SELECT. EQ. '1') THEN 
CALL CLS 
WRITE(* , 31142) 
31142  FORMAT(////,' ENTER THE WEIGHT FOR LAYER 1 (oz): ',2X, ) 
IF(UW. EQ. 'gm') THEN 
READ *,T1W1 
T1W1=T1W1*62. 5 
ELSE 
READ *,T1W1 
ENDIF 
VOL1=UN*T1W1 
T1L1=VOL1/AREA 
GOTO 31170 
ELSEIF(WORL. EQ. 'W'. AND. SELECT. EQ. '2' )THEN 
CALL CLS 
WRITE(*, 31143) 
31143. FORMAT(////,' ENTER THE WEIGHT FOR LAYER 1 (oz): ',2X, ) 
IF(UW. EQ. 'gm') THEN 
READ *,T2W1 
T2W1=T2W1*62. 5 
ELSE 
READ *,T2W1 
ENDIF 
VOL1=UN*T2W1 
T2L1=VOL1/AREA 
WRITE(*, 31144) 
31144 FORMAT(/,' ENTER THE WEIGHT FOR LAYER 2 (oz): ',2X, ) 
IF( UW. EQ. 'gm') THEN 
READ *,T2W2 
T2W2=T2W2*62. 5 
ELSE 
READ * ,T2W2 
ENDIF 
VOL2=UN*T2W2 
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31145 


31146 


31147 


31148 


31149 


T2L2=VOL2/AREA 
GOONS 76 
ELSEIF(WORL. EQ. 'W'. AND. SELECT. EQ. '3' )THEN 
CALICCLS 
WRITE(* , 31145) 


FORMAT(////,' ENTER THE WEIGHT FOR LAYER 1 (oz): ',2X, ) 


IF(UW. EQ. 'gm') THEN 
READ *,T3W1 
T3W1=T3W1*62. 5 
ELSE 
READ *,T3W1 
ENDIF 
VOL1=UN*T3W1 
T3L1=VOL1/AREA 
WRITE(*,31146) 
FORMAT(/,' ENTER THE WEIGHT FOR LAYER 2 (oz): 
IF(UW. EQ. 'gm') THEN 
READ *,T3W2 
T3W2=T3W2*62. 5 
ELSE 
READ *,T3W2 
ENDIF 
VOL2=UN*T3W2 
T3L2=VOL2/AREA 
WRITE(*, 31147) 
FORMAT(/,' ENTER THE WEIGHT FOR LAYER 3 (02): 
IF( UW. EQ. 'gm') THEN 
READ *,T3W3 
T3W3=T3W3*62. 5 
ELSE 
READ *,T3W3 
ENDIF 
VOL3=UN*T3W3 
T3L3=VOL3/AREA 

GOTO 31170 
ELSEIF(WORL. EQ. 'W'. AND. SELECT. EQ. '4' )THEN 
CALIMCES. 
WRITE(*, 31148) 


» 92K, ) 


FORMAT(////,' ENTER THE WEIGHT FOR LAYER 1 (oz): ',2X, ) 


IF(UW. EQ. 'gm') THEN 

READ *,T4W1 

T4W1=T4W1*62. 5 

ELSE 

READ *,T4W1 

ENDIF 

VOL1=UN*T4W1 

T4L1=VOL1/ AREA 
WRITE(* ,31149) 

FORMAT(/,' ENTER THE WEIGHT FOR LAYER 2 (oz): 

IF(UW. EQ. 'gm') THEN 

READ *,T4W2 

T4W2=T4W2*62. 5 

ELSE 

READ *,T4W2 

ENDIF 

VOL2=UN*T4W2 
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ey ) 


T4L2=VOL2/AREA 
WRITE(*, 31150) 
31150 FORMAT(/, ’ ENTER THE WEIGHT FOR LAYER 3 (0z): ',2X, ) 
IF( UW. EQ. 'gm') THEN 
READ *,T4W3 
T4W3=T4W3*62. 5 
ELSE 
READ * ,T4W3 
ENDIF 
VOL3=UN*T4W3 
T4L3=VOL3 / AREA 
WRITE(* , 31151) 
BEi51 FORMAT(/, - ENTER THE WEIGHT FOR LAYER 4 (0z): ',2X, ) 
IF(UW. EQ. 'gm') THEN 
READ *,T4W4 
T4W4=T4W4*62. 5 
BLE 
READ *,T4W4 
ENDIF 
VOL4=UN*T4W4 
T4L4=VOL4/AREA 
GOrOn31170 
CALL CLS 
GOTOy3 11.70 
ELSE 
CALL CLS 
Gea@ «31171 
ENDIF 
31170 IF(SPEVAL. EQ. 'B') THEN 
GOTO 41186 
ELSE 
CONTINUE 
ENDIF 
WRITE (31172) UK 
31172 FORMAT(/,' ENTER COPPER LAYER THERMAL CONDUCTIVITY (',A10,'): °, 
os, + 
READ *,UKC 


C 
C MAKE CHANGES OR CORRECTIONS TO COPPER LAYER ENTRIES 
C 
§) 


1193 IF(SELECT. EQ. '1' )THEN 


CALL. CLS 
WRITE(*,31173) CL, LOCO 
31173 FORMAT(///,' YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE COPPE 
+R LAYER(S). ee) , 
+' em) LENGTH: se ey et Cy Wei), «) 
WRITE(*,73) CW,LOCO 
73 FORMAT(/,' 2.) WIDTH: oF9)4, 1X ene 2k e) 
WRITE(*,74) T1L1, LOCO 
74 FORMAT(/,° 3.) THICKNESS LAYER 1: ',F9.4,1X,A2,2X, ) 
WRITE(*,75) UKC,UK 
75 FORMAT(/,' d. )_k: PFS 248 1X A105 2X9) 
ELSEIF( SELECT. EQ. '2')THEN 
GALL CLS 
WRITE(*, 311) CL,LOCO 
311  FORMAT(///,' YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE COPPE 
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+R SLAYER( S eae, (7 


+! 1.) LENGTH: " pO eGel erAee X 
WRITE(*,312) CW,LOCO 
312 FORMAT(/,' 2.) WIDTH: " ,F9.4,1X,A202 Xe 
WRITE(*,314) T2L1, LOCO 
314  FORMAT(/,' 3.) THICKNESS LAYER 1:  ,F9.4, WAR 2 
WRITE(* (315) T2L2 eco 
315 FORMAT(/,' THICKNESS LAYER 2: ',F9.4,1X,A2,2X, ) 
WRITE(*,316) UKC,UK 
316 FORMAT(/,' 4. ) k: ' ,F9. 4, 1X JA10 2 
ELSEIF( SELECT. EQ. '3' )THEN 
CALI CLS 
WRITE(*,317) CL,LOCO 
suly) FORMAT(///, ' YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE COPPE 
+R LAYER(S).',///, 
+! 1.) LENGTH: ' F9.4,1%,A2, 25 
WRITE(*,318) CW,LOCO 
318 FORMAT(/,' 2.) WIDTH: ' F9.4,1X%,A2,2% 
WRITE(* 599) [3L1,L0co 
319 FORMAT(/,' 3.) THICKNESS LAYER 1: * ,F9. 4 ikeAg= 2 
WRILE(* 520 132 snoOee 
320 FORMAT(/,' THICKNESS LAYER 2: ~ ,F9.4,1X,A297x 
WRITE(*,321) T3L3,LOCO 
321 FORMAT(/,' THICKNESS LAYER 3: ',F9.4,1X,A2,2X, ) 
WRITE(*, 322) UKC,UK 
322 FORMAT(/,' 4.) k: " ,F9.4,1X,A105 2x) 
ELSEIF( SELECT. EQ. '4' )THEN 
CAEEACES 
WRITE(*, 323) CL, LOCO 
323 FORMAT(/// 5 | YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE COPPE 
+R LAYER(S).',///, 
+! 1.) LENGTH: "FS. 4.1%, 42-2 
WRITE(*,324) CW, LOCO 
324 FORMAT(/,' 2.) WIDTH: ',F9. 49x, A2, 20a 
WRITEC* 325) 141 nOco 
325 FORMAT(/,' 3.) THICKNESS LAYER 1: ',F9.4,1X,A2,2X, ) 
WRITE(*,326) T4L2,LOCO 
326 FORMAT(/,' THICKNESS LAYER 2: ',F9.4,1X,A2,2X, ) 
WRITE(* ,327) T4L3, LOCO 
327 (FORMATC THICKNESS LAYER 3: ',F9.4,1X,A2,2X, ) 
WRITE(*,328) T4L4, LOCO 
328 FORMAT(/,' THICKNESS LAYER 4: ',F9.4,1X,A2,2X, ) 
WRITE(*, 329) UKC,UK 
329  FORMAT(/,' AM kes ' ,F9. 4,1X,A10,2 
ENDIF 
C 
31194 WRITE(*,71180) 
71180 FORMAT(///,' DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR Y 
+ES AND N FOR NO: ',2X, ) 
READ(*,31181)ANSE1 


31181 FORMAT(A1) 
IF(ANSE1.EQ.'Y') THEN 
56342  WRITE(*,31184) 
31184  FORMAT(/," WOULD YOU LIKE TO CHANGE THE THICKNESS? (Y OR N 
+e 
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READ(*,31185)ANST1 
31185 FORMAT( Al) 
PRINT * 
IF(ANST1. EQ. 'Y') THEN 
SPEVAL='B' 
CAEL CLS 
WRITE(* ,99998) 
99998 FORMAT(////) 
GOTO 99999 
ELSEIF(ANST1. EQ. 'N') THEN 
CONTINUE 
ELSE 
Cc0i.CLS 
WRITE(* , 56343) 
56343 FORMAT( ////) 
GOTO 56342 
ENDIF 
C 
41186 GALL CLS 
WRITE(*, 31187) 
31187 FORMAT(////) 
mee THE CURRENT ENTRY FOR THERMAL CONDUCTIVITY IS ',UK 
+C, , UK 
32186 WRITE(*, 31188) 
31188 FORMAT(/,' WOULD YOU LIKE TO CHANGE THE THERMAL CONDUCTIVIT 
+Y¥? (Y OR N): ',2X, ) 
READ(*,31189)ANSK1 
31189 FORMAT( Al) 
PRINT * 
IF(ANSK1. EQ. 'Y') THEN 
WRITE(*,31190) UK 
31190 FORMAT(/,' ENTER THE THERMAL CONDUCTIVITY (',A10,'): ',2X, ) 
READ *,UKC 
ELSEIF(ANSK1. EQ. 'N') THEN 
GOTO 31191 
ELSE 
CALL CLS 
WRITE(* , 56344) 
56344 FORMAT( ////) 
GOTO 32186 
ENDIF 
31191 GOTO 31193 


ELSEIF(ANSE1. EQ. 'N') THEN 
GOTO 31180 
ELSE 
CALL. CLS 
GOTO 31194 
ENDIF 
31180 END 
Crete tedededete dete dete dete dete tek ete de KPI III ILD TELLER IRE TERI RRR IRR REITER ERR ERE EERE EE 


C 


C pies: MODEL BUILDER 
C SUBROUTINE: S1 
C DATE: GF JUL 91 
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NIOaaIaIantian 


898 


803 


AUTHOR: LT STEVE GLASER 
COMPILER: MICROSOFT VERSION 4.01 
LINKER: MICROSOFT VERSION 3.55 


CALLED FROM SUBROUTINE PCBS3. THIS SUBROUTINE GENERATES THE OUTPUT 
DATA FILE FOR THE ONE COPPER LAYER CASE 


SUBROUTINE S1 (GYLR,GFXB,GYY,GXX,GZC11,GZB ,GYLR11,GFXB11,GYY11, GXX 
+11,GZE11,GPL,GIH,GJH,GEAT, GWIDE ,GDEEP,GOEF ,GIBT,GPRT ,GWRT,GRT,GLT, 
+GFT,GBT,GSEL) 


REAL GYLR,GFXB ,GYY ,GXX,GZC11,GZB ,GYLR11,GFXB11,GYY11,GXX11,GZE11,A 
+CC , DAMP ,CONFAC ,GIBT ,GPRT ,GWRT ,GRT ,GLT ,GFT, GBT 
REAL GEAT( 100,100) ,GOEF( 1000, 12) 


CHARACTER*1 SELECT, ANS 
CHARACTER DATAF*79 , NAME*6 


INTEGER GPL,I,N,IB,GWIDE,GDEEP , COUNT , CONTEMP , ZER ,GSEL,NMAX , TMAX , HT 
+RS,D1,D2,D3,D4,D5,D6,D7 ,MAXIT, LOCVAR 


INTEGER GIH( 1000) ,GJH( 1000) ,GCON(1000,12) 
CONTEMP=6 
ZER=0 
NMAX=750 
TMAX=50 
HTRS=6 

D1=2 

D2=4 

D3=6 

D4=0 

D5=0 

D6=0 

D7=0 

ACC=0. 05 
DAMP=0. 666667 
MAXIT=12 
CONFAC=0. 8 
COUNT=2*GPL 
LOCVAR=COUNT 


CALL CLS 


WRITE(*, 803) 
FORMAT(///,' THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY 
+ INTO THE',/, 


+! EXISTING THERMAL ANALYZER, FURTHERMORE, THIS PROGRAM DOES',/, 
+" NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE’ ,/, 
ft THE USER WILL NAME THE DATA FILE FOR EACH RUN OF THIS we 
2uM PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS, AND’ ,/, 
4 SHOULD NOT HAVE ANY SPACES. "FJ fee 
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+! PLEASE ENTER THE DESIRED DATA FILE NAME: ',2X ) 
READ(* ,804) NAME 

804  FORMAT(A6) 

552 WRITE(*,910)NAME 

D:D FORMAT(///// ,' YOU SELECTED ',A6,' FOR YOUR DATA FILE NAME 
+ 


a / 
911 WRITE(* , 811) 
811 FORMAT( IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N 
+FOR NO: 7) 
READ(*,812) ANS 
812 FORMAT(A1) 


IF(ANS. EQ. 'N') THEN 
GOTO 898 

ELSE 
CONTINUE 

ENDIF 

IF(ANS. EQ. 'Y') THEN 
GOTO 897 

ELSE 
CALE CLS 
GOTO 552 

897 ENDIF 


C ALLOW THE USER TO PROVIDE TITLE LINE FOR DATA FILE 


e535) GALL CLS 
WRITE(* ,805) 
805  FORMAT(////,' ENTER THE DESIRED TITLE TO BE PLACED ON LINE 
ee / _ NUMBER ONE OF THE OUTPUT DATA BILE: 9/7 /, 
+ 2, 
READ(*,806) DATAF 
806 FORMAT(A79) 
835 WRITE(*,831) 
831 FORMAT(////,' DO YOU WISH TO CHANGE THE TITLE OF YOUR OUTPU 
+T DATA FILE?',/,' ENTER Y FOR YES AND N FOR NO: WaaXx, + 
READ(*,832) ANS 
832 FORMAT(A1) 


IF(ANS. EQ. 'Y') THEN 
GOTO 833 

ELSE 
CONTINUE 

ENDIF 

IF(ANS. EQ. 'N') THEN 
GOTO 834 

ELSE 
GOTO 835 

834 ENDIF 


Creterevededede de vee dese te te Fe Fete Fede Fete FoF KCK HR RIK RRR RIKER RRR RE RRER EEREREREREEREREREER 


Citiicinineieeet COEFFICIENTS FOR EPOXY AND COPPER LAYERS*iheddditiicininninininins 
Crete deter de eve de te Fete te Riese te HF TIL ERIC REPRE TEPPER TEER TERRE LER TELERIK REE RRR EREEE ERE 


é 
DO 90 I=1,GPL 


eS) 


N=1 
IB=GPL+I 
C 
C 
CHARA KRRIKEREHREKRKKREEEEEEK CORNERS RARER KHAKI 
C 
IF ((GIH(1I). EQ. 1. OR. GIH(1). EQ. GDEEP). AND. (GJH(1). EQ. 1. OR. GJH(1). EQ 
+. GWIDE)) THEN 


C DETERMINE COEFFICIENTS FOR TOP LAYER 
IF (GEATCGIH(U1) ,GJH(1}>). EG: 0 Oe. 
GCONCI,N) = 6 


ELSE 
GCON(I,N) = 7 
ENDIF 
C 
C CONNECTIONS FOR EPOXY LAYER 
GCON(IB,N) = 6 
C 
é 
C LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
6 
C LEFT EDGE 
IF (GJH(1).EQ.1) THEN 
C LEFT COEFFICIENT 
GOEF(I,N) = GYLR11 
GOEF(IB,N) = GYLR 
N=N+1 
GCON(I,N) = 7551 
GCON(IB,N) = 7551 
C 


6 RIGHT COEFFICIENT 
GOEF(I,N) = GYY11 
GOEF(IB,N) = GYY 
N=N+1 
GCON(I,N) = 10 * (I+1) + 1 
GCON(IB,N) = ‘10* (IB+1) + 1 


C RIGHT EDGE 
ELSEIF (GJH(1I).EQ.GWIDE) THEN 


C 
C LEFT COEFFICIENT 
C 
GOEF(I,N) = GYY11 
GOEF(IB,N) = GYY 
N=N+1 
GCON(I,N) = 10*(I-1)+1 
GCON(IB,N) = 10*(IB-1)+1 
C 
C RIGHT COEFFICIENT 
C 


GOEF(I,N) = GYLR11 
GOEF(IB,N) = GYLR 
=N+1 

GCON(I,N) = 7541 
GCON(IB,N) = 7541 
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7) qoaoqgagy 


eS | 


ENDIF 


FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 


FRONT EDGE 
me CGinGl}. EQ, 1) THEN 
FRONT COEFFICIENT 
GOEF( I ,N)=GFXB11 
GOEF(IB,N)=GFXB 
N=N+1 
GCONCI,N) =7521 
GCONCIB,N) = 7521 


BACK COEFFICIENT 
GOEF(I,N)=GXX11 
GOEF( IB ,N)=GXX 
N=N+1 
GCON(I,N) =10*( I+GWIDE)+1 
GCONCIB,N) = 10*( IB+GWIDE)+1 


BACK EDGE 
ELSEIF (GIH(1).EQ.GDEEP) THEN 
FRONT COEFFICIENT 
GOEF(I,N)=GXX11 
GOEF( IB ,N)=GXX 
=N+1 
GCON(I,N) =10*( 1-GWIDE)+1 
GCON(CIB,N) = 10*( IB-GWIDE)+1 


BACK COEFFICIENT 
GOEF(I,N)=GFXB11 
GOEF(IB,N)=GFXB 
N=N+1 
GCON(CI,N) = 7511 
GCON(IB,N) = 7511 
ENDIF 


TOP COEFFICIENT 
GOEF(I,N)=GZE11 
GOEF( IB ,N)=GZC11 
N=N+1 
GCONCI,N) = 7511 
GCONCIB,N) = 10*I+1 


BOTTOM COEFFICIENT 


GOEF(I,N) = GZC11 
GOEFCIB,N) = GZB 

=N+1 

GCON(I,N) =10*( I+GPL)+1 
GCONCIB,N) = 10*( IB+GPL)+1 


HEAT INPUT 
IF(GEAT(GIH(1) ,GJH(I)).NE. 0.0) THEN 
GOEF(I,N) = GEAT(GIH(1) ,GJH(I)) 
N=N+1 
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aqIanan aQao 


Cl} G2a°Garc3 


© 


GCON(1,N)=9991 
ENDIF 


ELSEIF((GIH(I). EQ. 1. OR. GIH(1). EQ. GDEEP). AND. (GJH( 1). NEW ORSGoeee 


+. NE. GWIDE)) THEN 


DETERMINE NUMBER OF CONNECTIONS FOR COPPER AND EPOXY LAYERS 


IF (GEAT(GIH(1),GJHC1I)).EQ. 0.0) THEN 
GCONCI,N) = 6 

ELSE 
GCONCI,N) = 7 

ENO IE 


CONNECTIONS FOR EPOXY LAYER 
GCON(CIB,N) = 6 


LEFT AND RIGHT COEFFICIENTS 
LEFT COEFFICIENT 

GOEF(I,N) = GYY11 
GOEF(IB,N) = GYY 

N=N+1 

GCON(I,N) = 10*(I-1)+1 
GCON(IB,N) = 10%*(IB-1)+1 


RIGHT COEFFICIENT 
GOEF(I,N) = GYY11 
GOEF(IB,N) = GYY 
N=N+1 
GCON(I,N) = 10 * (I+1) + 1 
GCON(IB,N) = 10* (IB+1) + 1 


FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 


FRONT EDGE 

IF (GIH(I).EQ.1) THEN 
FRONT COEFFICIENT 
GOEF(1,N)=GFXB11 
GOEF( 1B,N)=GFXB 

N=N+1 

GCON(I,N) =7521 
GCON(IB,N) = 7521 


BACK COEFFICIENT 

GOEF( I,N)=GXX11 

GOEF( 1IB,N)=GXxX 

N=N+1 

GCONCI,N) =10*( I+GWIDE)+1 
GCONCIB,N) = 10*( IB+GWIDE)+1 


BACK EDGE 

ELSEIF (GIH(1). EQ. GDEEF) itm 
FRONT COEFFICIENT 
GOEF(1I,N)=GXX11 
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GOEF( IB ,N)=GXX 

=N+1 

GCON(I,N) =10*( I-GWIDE)+1 
GCON(IB,N) = 10*( IB-GWIDE)+1 


C BACK COEFFICIENT 
GOEF(I,N)=GFXB11 
GOEF( IB ,N)=GFXB 

=N+1 

GCON(I,N) = 7511 
GCON(IB,N) = 7511 
ENDIF 


C TOP COEFFICIENT 
GOEF(I,N)=GZE11 
GOEF( IB ,N)=GZC11 

=N+1 
GCON(I,N) = 7511 
GCONCIB,N) = 10*I+1 


C 
C BOTTOM COEFFICIENT 
C 
GOEF(I,N) = GZC11 
GOEF(IB,N) = GZB 
N=N+1 
GCON(I,N) =10*(I+GPL)+1 
GCON(IB,N) = 10%( IB+GPL)+1 
€ 


C HEAT INPUT 
IF(GEAT(GIH(I),GJH(I)).NE. 0.0) THEN 
GOEF(I,N) = GEAT(GIH(1I),GJH(I)) 
N=N+1 
GCON(I,N)=9991 
ENDIF 


CXCXCXC 


LEFT AND RIGHT EDGES EXCLUDING CORNERS 


aqaaqa 


meee eeCGIHC) 2 EQ. 1. OR. GIH( 1). EQ. GDEEP). AND. (GJH(1). EQ. 1. OR. GJH( 1) 
+.EQ.GWIDE)) THEN 


ao? 


DETERMINE COEFFICIENTS FOR TOP LAYER 
IF (GEAT(GIH(I),GJH(1)).NE.0.0) THEN 
GCON(I,N) = 7 
ELSE 
GCON(I,N) = 6 
ENDIF 


6 CONNECTIONS FOR EPOXY LAYER 
GCON(IB,N) = 6 


aqaag 


LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
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Ge lS C37 C2 © 


LP jal a> 


LEFT EDGE 

cP CGI LD PEG ye tHe 
LEFT COEFFICIENT 
GOEF(I,N) = GYLR11 
GOEF(IB,N) = GYLR 
N=N+1 

GCONCI,N) = 7551 
GCONCIB,N) = 7551 


RIGHT COEFFICIENT 
GOEF(I,N) = GYY11 
GOEF(IB,N) = GYY 
N=N+1 
GCON(CI,N) = 10 * (I+1) + 1 
GCONCIB,N) = 10* (IB+1) + 1 


RIGHT EDGE 
ELSEIF (GJH(1).EQ.GWIDE) THEN 


LEFT COEFF ICIENDL 


GOEF(I,N) = GYY11 
GOEF(IB,N) = GYY 

N=N+1 

GCON(I,N) = 10*(I-1)+41 
GCON(IB,N) = 10*(IB-1)+1 


RIGHT COEFFICIENT 


GOEF(I,N) = GYLR11 
GOEF(IB,N) = GYLR 
N=N+1 

GCON(I,N) = 7541 
GCON(IB,N) = 7541 
ENDIF 


FRONT COEFFICIENT 
GOEF(1,N)=GXX11 
GOEF( IB ,N)=GXX 
N=N+1 


GCONCI,N) =10*( I-GWIDE)+1 
GCON(CIB,N) = 10*( IB-GWIDE)+1 


BACK COEFFICIENT 
GOEF(I,N)=GXX11 
GOEF( IB ,N)=GXX 
N=N+1 
GCON(I,N) =10%*( I+GWIDE)+1 
GCON(IB,N) = 10*( IB+GWIDE)+1 


TOP COEFFICIENT 
GOEF(1I,N)=GZE11 
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GOEF(IB,N)=GZC11 
N=N+1 

GCON(I,N) = 7511 
GCON(IB,N) = 10*I+1 


C 
c BOTTOM COEFFICIENT 
C 
GOEF(I,N) = GZC11 
GOEF(IB,N) = GZB 
N=N+1 
GCON(I,N) =10*(I+GPL)+1 
GCON(IB,N) = 10*( IB+GPL)+1 
C 


C HEAT INPUT 
IF(GEAT(GIH(I),GJH(I)).NE. 0.0) THEN 
GOEF(I,N) = GEAT(GIH(I),GJH(I)) 
N=N+1 
GCON(I,N)=9991 
ENDIF 


KKEEDETERMINE COEFFICIENTS FOR ALL NODES NOT TOUCHING AN EDGE* xix 


(> 1 > > ja ap 


Peet C(GIHCT)]NE. 1. OR. GIH(I).NE. GDEEP). AND. (GJH(1). NE. 1. OR. GJH(I) 
+. NE.GWIDE)) THEN 


Pe) 


DETERMINE CONNECTIONS FOR TOP LAYER 
IF (GEAT(GIH(I) ,GJH(I)).NE.0.0) THEN 
GCON(I,N) = 7 
ELSE 
GCON(I,N) = 6 
ENDIF 


Pep] 


CONNECTIONS FOR EPOXY LAYER 
GCONCIB,N) = 6 


Sp ha > Ui ie a a> 


LEFT COEFFICIENT 
GOEF(I,N) = GYY11 
GOEF(IB,N) = GYY 
N=N+1 
GEONCE SN) = 10*(I-1)+1 
GCON(CIB,N) = 10*(IB-1)+1 


C RIGHT COEFFICIENT 
GOEF(I,N) = GYY11 
GOEF(IB,N) = GYY 
N=N+1 
GCON(I,N) = 10 * (I+1) + 1 
GCON(IB,N) = 10% (IB+1) + 1 
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C FRONT COEFFICIENT 
GOEF(I,N)=GXX11 
GOEF( IB ,N)=GXX 
N=N+1 
GCON(I,N) =10*( I-GWIDE)+1 
GCONCIB,N) = 10*( IB-GWIDE)+1 


C BACK COEFFICIENT 
GOEF(I,N)=GXX11 
GOEF(IB,N)=GXX 
N=N+1 
GCON(I,N) =10*( I+GWIDE)+1 
GCONCIB,N) = 10*( IB+GWIDE)+1 


C TOP COEFFICIENT 
GOEF( I ,N)=GZE11 
GOEF(IB,N)=GZC11 
N=N+1 
GCONGI,N) =a7511 
GCONCIB,N) = 10*I+1 


C 
C BOTTOM COEFFICIENT 
C 
GOEF(I,N) = G2ZCl1l1 
GOEF(IB,N) = GZB 
N=N+1 
GCON(CI,N) =10*( I+GPL)+1 
GCONCIB,N) = 10*( IB+GPL)+1 
C 


G HEAT INPUT 
IF(GEAT(GIH(1I) ,GJH(I)).NE.0.0) THEN 
GOEF(I,N) = GEAT(GIH(1I),GJH(I)) 
N=N+1 
GCON(I,N)=999 
ENDIF 7 
ENDIF 

0 CONTINUE 


9 
C 
C GENERATE DATA FILE 
C 


OPEN (3,FILE=NAME , FORM='FORMATTED' ,ACCESS=' DIRECT’ , RECL=108 , STATUS 

+='NEW' ) 

WRITE( 3,909) DATAF 
909  FORMAT(1X,A79) 

WRITE( 3,908) COUNT,CONTEMP, ZER, ZER, ZER, ZER, ZER, ZER, GSEL 
908  FORMAT(2X,9(13,5X)) 

WRITE( 3,907) ZER,ZER, ZER 
907 FORMAT(2X,3(13,5X)) 

WRITE(3,9081) NMAX,TMAX,HTRS,D1,D2,D3,D4,D5,D6,D7 
9081 FORMAT(2X,9(13,5X)) 

WRITE(3,905) ACC,DAMP,MAXIT,CONFAC ,GIBT 
905  FORMAT(1X,2(F9.7,1X),14X,12,1X, 59) 7 ee 

WRITE( 3,906) GPRT,GFT,GBT,GRT,GLT,GWRT 
906  FORMAT(1X,6(F12. 4,1X)) 
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DO 112 I=1,LOCVAR 
WRITE(3,9100) GCON(I,1),GCON(I,2),GCON(I,3),GCON(I,4),GCON(I,5),GC 
+ON(1,6),GCON(I,7),GCON(I,8) 

9100 FORMAT(14,3X,7(15,7X)) 
WRITE(3,9110) GOEF(I,1),GOEF(I,2),GOEF(1,3),GOEF(1,4),GOEF(I,5),GO 
+EF(1,6),GOEF(I,7) 

9110 FORMAT(7(F9. 3,3X)) 

112 CONTINUE 


CLOSE (3) 
CALL CLS 
WRITE(*,999) NAME 
999 FORMAT(////,' THE OUTPUT DATA HAS BEEN PLACED IN A FILE 
+NAMED ',A6,/////,° <PRESS ENTER TO CONTINUE>' ) 


READ(*,5912) ANS 
5912 FORMAT(A1) 


Cre I I He IH eI Fe Fe Fe FH Fe IH FHI IH FHI TCR FCI TCR TET KE TCR FRICKE TRE IR KEKE ERE 


S2 IS CALLED FROM PCBS3. S2 IS THE OUTPUT DATA FILE GENERATOR FOR THE 
THE COPPER LAYER CASE. 


C 

C TITLE: MODEL BUILDER 

C SUBROUTINE: S2 

C DATE: 09 JUL 91 

C AUTHOR: LT STEVE GLASER 

C COMPILER: MICROSOFT VERSION 4.01 
C LINKER: MICROSOFT VERSION 3.55 
C 

C 

C 

C 


SUBROUTINE S2 (GYLR,GFXB,GYY,GXX ,GZC21,GZC22 ,GZB,GYLR21,GYLR22 ,GFX 
+B21,GFXB22 ,GYY21,GYY22,GXX21,GXX22 ,GZE21,GZE22,GPL,GIH,GJH,GEAT,GW 
+IDE ,GDEEP ,GOEF ,GIBT ,GPRT,GWRT,GRT,GLT,GFT,GBT,GSEL) 


REAL GYLR,GFXB,GYY ,GXX ,GZC21,GZB ,GYLR21,GFXB21,GYY21,GXX21,GZE21,A 
+CC , DAMP , CONFAC ,GIBT,GPRT,GWRT,GRT,GLT,GFT,GBT,GYLR22 ,GFXB22,GYY22, 
ESO ONGZE22 .GZC22 


REAL GEAT( 100,100) ,GOEF( 1000, 12) 


CHARACTER*1 SELECT, ANS 
CHARACTER DATAF*79 , NAME*6 


INTEGER GPL,I,N,1B,GWIDE,GDEEP , COUNT, CONTEMP , ZER,GSEL,NMAX , TMAX, HT 
+RS ,D1,D2,D3,D4,D5 ,D6,D7 ,MAXIT, LOCVAR,IC, ID 


INTEGER GIH( 1000) ,GJH( 1000) ,GCON(1000,12) 
CONTEMP=6 

ZER=0 

NMAX=750 

TMAX=50 

HTRS=6 

D1=2 

D2=4 

D3=6 

D4=0 
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805 


806 


D5=0 

D6=0 

D7=0 

ACC=0. 05 
DAMP=0. 666667 
MAXIT=12 
CONFAC=0. 8 
COUNT=4*GPL 
LOCVAR=COUNT 


CALL CLS 


WRITE(*, 803) 
FORMAT(///,' 
+ INITIO THE ye 


— 


THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY 


EXISTING THERMAL ANALYZER, FURTHERMORE, THIS PROGRAM DOES',/, 


NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE',/, 


: THE USER WI 


SHOULD NOT 
PLEASE ENTE 
READ(*,804) NAME 
FORMAT( A6) 
WRITE(* ,910)NAME 
FORMAT(/////,' 

+ araie) 
WRITE(* , 
FORMAT( 

+FOR NO: Oe. 
READ(*,812) ANS 
FORMAT(A1) 


811) 


LL NAME THE DATA FILE FOR EACH RUN OF THIS o. 


HAVE ANY SPACES. ' Je 
R THE DESIRED DATA FILE NAME: ',2X ) 


ie 
ae, 
ay, 
ui PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS, AND',/, 
ae 
a 


YOU SELECTED ',A6,' FOR YOUR DATA FILE NAME 


IS THIS THE DESIRED SELECTION? ENIER Y FOR YEouanva 


) 


IF(ANS. EQ. 'N') THEN 


GOTO 898 
ELSE 

CONTINUE 
ENDIF 


IF(ANS. EQ. 'Y') THEN 


GOTO 897 
ELSE 
CALL CLS 
GOTO 552 
ENDIF 


ALLOW THE USER TO PROVIDE TITLE LINE FOR DATA FILE 


GALL CLS 

WRITE(*, 805) 

FORMAT(////,' ENTER THE DESIRED TITLE TO BE PLACED ON LINE 
uae NUMBER ONE OF THE OUTPUT DATA FILE:',///, 
BE eae) 


READ(*,806) DATAF 


FORMAT(A79 ) 
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835 WRITE(*,831) 
831 FORMAT(////,' DO YOU WISH TO CHANGE THE TITLE OF YOUR OUTPU 
fe DATA BILE?’ ./,' ENTER Y FOR YES AND N FOR NO: mex, ) 
READ(* ,832) ANS 
832 FORMAT(A1) 
C 
IF( ANS. EQ. 'Y') THEN 
GOTO 833 
ELSE 
CONTINUE 
ENDIF 
IF(ANS. EQ. 'N') THEN 
GOTO 834 
ELSE 
GOTO 835 
834 ENDIF 


C rede rede se dese seve He ve Fe He Fe se Fe Fee Fe Fe WF He Fe Fe He Fe VEN Fe He Fe Fe He Fe He Fe Fe HE Fe NE We TE He FE Fee Te HE I TENET Fe IE I Fe ICH Ie HE I Fe He FE IH 


Cavietedcsedesesesede eee COEFFICIENTS FOR EPOXY AND COPPER LAYER Sveitvestiedevedcdedcsedeseaesedese 
Caksbikdededederededesedededeiesedeicdedeacded sedesederede dered sede aede seeder escseseiedededeseve sede dedcsedeaesevederede seve de severe dew ye 


C 
DO 90 I[=1,GPL 
N=1 
IB=GPL+tI 
IC=2*GPL+I 
ID=3*GPL+I 

C 


CHR IKKAKKIKKK  KKIKE R  E R  K CORNER SHRINK IIR IIH KIKI IKK IH KI HII 
C 
ie (Giant). EQ. 1. OR. GIA( 1). EQ. GDEEP). AND. (GJH(1). EQ. 1. OR. GJH(I). EQ 
+. GWIDE)) THEN 


C DETERMINE COEFFICIENTS FOR TOP COPPER 
IF (GEAT(GIH(I),GJHC(I)).EQ.0.0) THEN 
GCON(I,N) = 6 


ELSE 
GCON(I,N) = 7 
ENDIF 
6 
C CONNECTIONS FOR EPOXY LAYER 
GCON(IB,N) = 6 
GCON(ID,N) = 6 
GCON(IC,N) = 6 
¢ 
C 
C LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
c 
C LEFT EDGE 
IF (GJH(I).EQ. 1) THEN 
C LEFT COEFFICIENT 
GOEF(I,N) = GYLR21 
GOEF(IC,N) = GYLR22 
GOEF(IB,N) = GYLR 
GOEF(ID,N) = GYLR 
N=N+1 
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oo 


gqon 


C 


C2 c) ©) 


‘?) 


GCON(I,N) = 7551 


GCONCIB,N) = 7551 
GCONCIC,N) = 7551 
GCON(CID,N) = 7551 


RIGHT COEFFICIENT 
GOEF(I,N) = GYY21 


GOEF(IC,N) = GYY22 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
N=N+1 


GCONCI,N) = 10 (ieee 


GCON(IB,N) = 10* (IB+1) + 1 

GCON(IC,N) = 10 * (IC+1) +1 

GCON(ID,N) = 10 * (ID+1) + 1 
RIGHT EDGE 


ELSEIF (GJH(I).EQ. GWIDE) THEN 
LEFT COEFFICIENT 
GOEFCI,N)) = Gr Y2i 


GOEF(IC,N) = GYY22 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 

N=N+1 

GCON(I,N) = 10*(I-1)+1 
GCON(IB,N) = 10*(IB-1)+1 
GCON(IC,N) = 10*(IC-1)+1 
GCON(ID,N) = 10*(ID-1)+1 


RIGHT COEFFICIENT 


GOEF(I,N) = GYLR21 


GOEF(IC,N) = GYLR22 
GOEF(IB,N) = GYLR 
GOEF(ID,N) = GYLR 
N=N+1 

GCON(I,N) = 7541 
GCON(IB,N) = 7541 
GCON(IC,N) = 7541 
GCON(ID,N) = 7541 
ENDIF 


FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 


FRONT EDGE 
IF (GIH(I).EQ.1) THEN 
FRONT COEFFICIENT 
GOEF(1I,N)=GFXB21 
GOEF( IC,N)=GFXB22 
GOEF( IB ,N)=GFXB 
GOEF( ID ,N)=GFXB 
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N=N+1 
GCON(CI,N) =7521 


iSGON( TBGN) = 7521 
GCONCIC,N) = 7521 
GCONCID,N) = 7521 


C BACK COEFFICIENT 
GOEF(I,N)=GXX21 
GOEF( IC ,N)=GXX22 
GOEF( 1B ,N)=GXX 
GOEF( ID ,N)=GXX 
N=N+1 
GCON(I,N) =10*( I+GWIDE)+1 
GCON(IB,N) = 10*( IB+GWIDE)+1 
GCON(IC,N) = 10*( IC+GWIDE)+1 
GCON(ID,N) = 10*( ID+GWIDE)+1 


C BACK EDGE 
ELSEIF (GIH(1I).EQ.GDEEP) THEN 
C FRONT COEFFICIENT 
GOEF(1I,N)=GXX21 
GOEF( IC ,N)=GXX22 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
N=N+1 
GCON(I,N) =10*( I-GWIDE)+1 


GCONCIB,N) = 10*( IB-GWIDE)+1 
GCONCIC,N) = 10*( IC-GWIDE)+1 
GCONCID,N) = 10*( ID-GWIDE)+1 


C BACK COEFFICIENT 
GOEF(1I,N)=GFXB21 
GOEFC IC ,N)=GFXB22 
GOEF( IB ,N)=GFXB 
GOEF(ID,N)=GFXB 


N=N+1 
GCON(I,N) = 7511 
GCON(IB,N) = 7511 
GCON(IC,N) = 7511 
GCON(ID,N) = 7511 
ENDIF 


C TOP COEFFICIENT 
GOEF(I,N)=GZE21 
GOEF( IC ,N)=GZE22 
GOEF(IB,N)=GZC21 
GOEF(ID,N)=G2C22 


N=N+1 

GCON(I,N) = 7511 
GCONCIC,N) = 7511 
GCONCIB,N) = 10*I+i 
GCONCID,N) = 10*IC+1 
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C 


qo 


1 > it Ge a a 


C 


BOTTOM COEFFICIENT 
GOEF(I,N) = GZC21 


GOEFCIC,N) = GZC22 
GOEF(IB,N) = GZE22 
GOEF(ID,N) = GZB 
N=N+1 


GCONCI,N) =10*(I+GPL)+1 
GCONCIB,N) = 10*( IB+GPL)+1 


GCON(IC,N) = 10*(IC+GPL)+1 | 
GCON(ID,N) = 10*(ID+GPL)+1 
HEAT INPUT : 


IF(GEAT(GIH(1I),GJH(I)).NE.0.0) THEN 
GOEF(I,N) = GEAT(GIH(I) ,GJH(I)) 
N=N+1 

GCON(I,N)=9991 

ENDIF 


ELSEIF((GIH(I). EQ. 1. OR. GIH( I). EQ. GDEEP). AND. (GJH(1). NE. 1) ORNGoiia 
+. NE; GWIDE) ) THEN 


DETERMINE NUMBER OF CONNECTIONS FOR COPPER AND EPOXY LAYERS 
IF (GEAT(GIH(1I),GJH(I)).EQ.0.0) THEN 
GCON(I,N) = 6 
ELSE 
GCON(I,N) = 7 
ENDIF 


CONNECTIONS 


HO 
GCON(IB,N) = 


R EPOXY LAYER 

6 
GCON(ID,N) = 6 
GCON(IC,N) = 6 


LEFT AND RIGHT COEFFICIENTS 
LEFT COEFFICIENT 
GOEF(I,N) = GYY21 


GOEFCIC,N) = GYY22 
GOEF(IB,N) = GYY 
GOEPCIDEN) = Gia 

N=N+1 

GCONCI,N) = 10*(I-1)+1 
GCONCIB,N) = 10*( TB=1) =i 
GCONCIC,N) = 10*(IC-1)+1 
GCONCID,N) = 10*(ID-1)+1 


RIGHT COEFFICIENT 
GOEF(I,N) = GYY21 
GOEF(CIC,N) Grea 2 
GOEF(IB,N) GYY 
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GOEF(ID,N) = GYY 
N=N+1 
EGON( TSN) = 10 * (I+1) + 1 


GCONCIB,N) = 10* (IB+1) + 1 
GCONCIC,N) = 10 * (IC+1) + 1 
GECONCIDSN) = 10 * (ID+1) + 1 


FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 


FRONT EDGE 
IF (GIH(I).EQ.1) THEN 
FRONT COEFFICIENT 
GOEF(I,N)=GFXB21 
GOEF( IC,N)=GFXB22 
GOEF( IB,N)=GFXB 
GOEF( ID,N)=GFXB 
N=N+1 
GCON(I,N) =7521 
GCON(IB,N) = 7521 
GEON(IC,N) = 7521 
EsONCID,N) = 7521 


>) Be ee al ae Jal ae 


@ BACK COEFFICIENT 
GOEF(1I,N)=GXxX21 
GOEF( IC ,N)=GXX22 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 

=N+1 

GCON(I,N) =10*( I+GWIDE)+1 
GCON(IB,N) = 10*( IB+GWIDE)+1 
GCON(IC,N) = 10*(IC+GWIDE)+1 
GCON(ID,N) = 10*( ID+GWIDE)+1 


C BACK EDGE 
ELSEIF (GIH(1).EQ.GDEEP) THEN 
C FRONT COEFFICIENT 
GOEF(I,N)=GXX21 
GOEF( IC ,N)=GXX22 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
=N+1 
GCON(I,N) =10*(I-GWIDE)+1 
GCON(IB,N) = 10*( IB-GWIDE)+1 
GCON(IC,N) = 10*(IC-GWIDE)+1 
GCON(ID,N) = 10*(ID-GWIDE)+1 


C BACK COEFFICIENT 
GOEF(I,N)=GFXB21 
GOEF(CIC,N)=GFXB22 
GOEF(IB,N)=GFXB 
GOEF( ID,N)=GFXB 


N=N+1 

GCON(I,N) = 7511 
GCON(IB,N) = 7511 
GCON(IC,N) = 7511 
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GCONCID,N) = 7511 
ENDIF 


C TOP COEFFICIENT 
GOEF(I,N)=GZE21 
GOEF(C IC ,N)=GZE22 
GOEF( IB ,N)=GZC21 
GOEF( ID,N)=GZC22 


N=N+1 
GCON(I,N) = 7511 
GCONCIC,N) = 7511 
GCONCIB,N) = 10*I+1 
GCONCID,N) = 10*IC+1 

C 

C BOTTOM COEFFICIENT 

C 
GOEF(I,N) = GZC21 
GOEFCIC,N) = GZC22 
GOEF(IB,N) = GZB 
GOEFCID,N) = GZB 
N=N+1 
GCON(I,N) =10*( I+GPL)+1 
GCON(CIB,N) = 10*( IB+GPL)+1 
GCONCIC,N) = 10*( IC+GPL)+1 
GCONCID,N) = 10*( ID+GPL)+1 

C 


C HEAT INPUT 
IF(GEAT(GIH(1),GJH(I)).NE. 0.0) THEN 
GOEF(I,N) = GEAT(GIH(1),GJH(I)) 


N=N+1 
GCON(I,N)=9991 
ENDIF 
CXEXOXG 
CXCXCXC 
6 
C 
C LEFT AND RIGHT EDGES EXCLUDING CORNERS 
c 
ELSEIF((GIH(1). EQ. 1. OR. GIH(1). EQ. GDEEP). AND. (GJH(1). EQ. 1. OR. GJH(I) 
+.EQ.GWIDE)) THEN 
C 
C DETERMINE COEFFICIENTS FOR TOP LAYER 
IF (GEAT(GIH(I),GJH(1)).NE.0.0) THEN 
GCON(I,N) = 7 
ELSE 
GCON(I,N) = 6 
ENDIF 
C 
C CONNECTIONS FOR EPOXY LAYER 
GCON(IB,N) = 6 
GCON(IC,N) = 6 
GCON(ID,N) = 6 
é 
C 
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oe , eae 


C LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 

C 

C LEFT EDGE 
PCC Int). Oot) THEN 

C LEFT COEFFICIENT 
GOEF(I,N) = GYLR21 
GOEFCIC,N) = GYLR22 
GOEF(IB,N) = GYLR 
GOEF(ID,N) = GYLR 
N=N+1 
GCON(I,N) = 7551 
GCONCIB,N) = 7551 
GCONCIC,N) = 7551 
GCONCID,N) = 7551 

C 


C RIGHT COEFFICIENT 
GOEF(I,N) = GYY21 


GOEFCIC,N) = GYY22 
GOEFCIB,N) = GYY 
GOEF(ID,N) = GYY 
N=N+1 


GCONCI,N) = 10 * (I+1) + 1 


GCON(IB,N) = 10* (IB+1) + 1 
GCON(IC,N) = 10 * (IC+1) + 1 
GCON(ID,N) = 10 * (ID+1) + 1 
C 
C RIGHT EDGE 
ELSEIF (GJH(1I).EQ.GWIDE) THEN 
6 
6 LEFT COEFFICIENT 
GOEF(I,N) = GYY21 
BSEF(IC.N) = GYY22 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
N=N+1 
GCON(I,N) = 10*(I-1)+1 
GCON(IB,N) = 10*(IB-1)+1 
GCON(IC,N) = 10*(IC-1)41 
GCON(ID,N) = 10*(ID-1)+1 
@ 
C RIGHT COEFFICIENT 
C 
GOEF(I,N) = GYLR21 
GOEF(IC,N) = GYLR22 
GOEF(IB,N) = GYLR 
GOEF(ID,N) = GYLR 
N=N+1 
GCON(I,N) = 7541 
GCON(IB,N) = 7541 
GCON(IC,N) = 7541 
GCON(ID,N) = 7541 
ENDIF 
c 
C 


C FRONT COEFFICIENT 
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GOEF(I,N)=GXX21 
GOEF( IC ,N)=GXX22 
GOEF( IB ,N)=GXX 

GOEF( ID,N)=GXX 

N=N+1 

GCON(I,N) =10*( I-GWIDE)+1 
GCON(IB,N) = 10*(IB-GWIDE)+1 
GCON(IC,N) = 10*( IC-GWIDE)+1 
GCON(ID,N) = 10*( ID-GWIDE)+1 


C BACK COEFFICIENT 
GOEF(1I,N)=GXX21 
GOEF( IC ,N)=GXX22 
GOEF(IB,N)=GXX 
GOEF( ID ,N)=GXX 
N=N+1 
GCON(CI,N) =10*( I+GWIDE)+1 
GCONCIB,N) = 10*( IB+GWIDE)+1 
GCONCIC,N) = 10*( IC+GWIDE)+1 
GCONCID,N) = 10*( ID+GWIDE)+1 


C TOE CORFFICIENT 
GOEF(1I,N)=GZE21 
GOEFC IC ,N)=GZE22 
GOEF(IB,N)=GZC21 
GOEF( ID,N)=GZC22 


N=N+1 
GCON(I,N) = 7511 
GCON(IC,N) = 7511 
GCON(IB,N) = 10*I+1 
GCON(ID,N) = 10*IC+1 

C 

C BOTTOM COEFFICIENT 

C 
GOEF(I,N) = GZC21 
GOEF(IC,N) = GZC22 
GOEF(IB,N) = GZE22 
GOEF(ID,N) = GZB 
N=N+1 
GCON(1,N) =10*(I+GPL)+1 
GCON(IB,N) = 10*( IB+GPL)+1 
GCON(IC,N) = 10*( IC+GPL)+1 
GCON(ID,N) = 10*( ID+GPL)+1 

© 


C HEAT INPUT 
IF(GEAT(GIH(1I) ,GJH(1)).NE.0.0) THEN 
GOEF(I,N) = GEAT(GIH(1),GJH(I)) 
N=N+1 
GCON( I ,N)=9991 
ENDIF 


*eEXDETERMINE COEFFICIENTS FOR ALL NODES NOT TOUCHING AN EDGE**X%ri% 


G2) GC) ©) ©) 
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Seta GGliCl) Ne. OR. GIH( 1). NE. GDEEP). AND. (GJHC 1). NE. 1. OR. GJH(T) 
+.NE.GWIDE)) THEN 


C DETERMINE CONNECTIONS FOR TOP LAYER 
IF (GEAT(GIH(I),GJHCI)).NE.0.0) THEN 
GCONCI,N) = 7 


ELSE 
GCONCI,N) = 6 

ENDIF 

C 

C CONNECTIONS FOR EPOXY LAYER 
GCONCIB,N) = 6 
GCONCIC,N) = 6 
GCONCID,N) = 6 

C 

C 

C 

C LEFT COEFFICIENT 
Geprcti,N) = GYY21 
GOEFCIC,N) = GYY22 
GOEF(IB,N) = GYY 
GOEFCID,N) = GYY 
N=N+1 


ECONC TSN) = 10 * (I-11) + 1 


GCONCIB,N) = 10* (IB-1) + 1 
GCONCIC,N) = 10 * (IC-1) + 1 
GCONCID,N) = 10 * (ID-1) + 1 


C RIGHT COEFFICIENT 
GOEF(I,N) = GYY21 


GOEFCIC,N) = GYY22 
GOEFCIB,N) = GYY 
GOEF(ID,N) = GYY 
N=N+1 


GCON(I,N) = 10 * (I+1) + 1 


GCON(IB,N) = 10* (IB+1) + 1 
GCON(IC,N) = 10 * (IC+1) +1 
GCON(ID,N) = 10 * (ID+1) + 1 


FRONT COEFFICIENT 
GOEF(I ,N)=GXX21 
GOEF( IC ,N)=GXX22 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
=N+1 
GCON(I,N) =10*( I-GWIDE)+1 


GCONCIB,N) = 10*( IB-GWIDE)+1 
GCONCIC,N) = 10*(IC-GWIDE)+1 
GCONCID,N) = 10*( ID-GWIDE)+1 


C BACK COEFFICIENT 


It 


C 


GOEF( I ,N)=GXX21 
GOEF( IC ,N)=GXX22 
GOEF( IB ,N)=GXX 

GOEF( ID,N)=GXxX 

N=N+1 

GCONCI,N) =10*( I+GWIDE)+1 
GCONCIB,N) = 10*( IB+GWIDE)+1 
GCONCIC,N) = 10*( IC+GWIDE)+1 
GCONCID,N) = 10*( ID+GWIDE)+1 


TOP COEFFICIENT 
GOEF(I,N)=GZE21 
GOEF( IC,N)=GZE22 
GOEF( IB ,N)=GZC21 
GOEF(ID,N)=GZC22 


N=N+1 

GCON(CI,N) = 7511 
GCONCIB,N) = 10*I+1 
GCONCIC,N) = 7511 
GCONCID,N) = 10*IC+1 


BOTTOM COEFFICIENT 
GOEF(I,N) = GZC21 


GOEF(IC,N) = GZC22 
GOEF(IB,N) = GZB 
GOEFCID,N) = GZE22 
N=N+1 


GCON(I,N) =10*(I+GPL)+1 

GCON(IB,N) = 10*(IB+GPL)+1 
GCON(IC,N) = 10*(IC+GPL)+1 
GCON(ID,N) = 10*( ID+GPL)+1 


HEAT INPUT 


[FCGEAT(GIH( 1) ,GJHC1)). NE. 0. 0) SEN 


GOEFCI,N) = GEATCGIN( 1) (GING)? 


N=N+1 
GCON(I,N)=9991 
ENDIF 

ENDIF 

CONTINUE 


GENERATE DATA FILE 


OPEN (3,FILE=NAME ,FORM=' FORMATTED’ , ACCESS='DIRECT’ , RECL=108, STATUS 


+='NEW') 
WRITE(3,909) DATAF 
FORMAT( 1X, A79) 


WRITE( 3,908) COUNT ,CONTEMP, ZER, ZER, ZER, ZER, ZER, ZER, GSEL 


FORMAT( 2X,9(13,5X)) 
WRITE(3,907) ZER,ZER, ZER 
FORMAT( 2X, 3(13,5X)) 


112 


9081 


905 


906 


9100 


2110 
2 


999 


Sia 


CWC? COO as Ss a a 


WRITE(3,9081) NMAX,TMAX,HTRS,D1,D2,D3,D4,D5,D6,D7 

FORMAT( 2X,9(13,5X)) 

WRITE(3,905) ACC,DAMP,MAXIT,CONFAC,GIBT 

FORMAT( 1X,2(F9. 7,1X),14X,12,1X,F9. 7,1X,F9. 5) 

WRITE(3,906) GPRT,GFT,GBT,GRT,GLT,GWRT 

FORMAT(1X,6(F12. 3,1X)) 

DO 112 I=1,LOCVAR 

WRITE(3,9100) GCON(I,1),GCON(I,2),GCON(I,3),GCON(I,4),GCON(I,5),GC 
+0N(I,6),GCON(I,7),GCON(I,8) 

FORMAT( 14, 3X,7(15,7X)) 

WRITE(3,9110) GOEF(I,1),GOEF(1I,2),GOEF(1I,3),GOEF(I,4),GOEF(1I,5),GO 
+EF(I,6),GOEF(I,7) 


FORMAT( 7(F9. 3,3X)) 

CONTINUE 

CLOSE (3) 

CAE, CLS 

WRITE(* ,999) NAME 

FORMAT(////,' THE OUTPUT DATA HAS BEEN PLACED IN A FILE 
+NAMED ',A6,/////,' <PRESS ENTER TO CONTINUE>') 
READ(*,5912) ANS 

FORMAT(A1) 

END 
TITLE: MODEL BUILDER 

SUBROUTINE: $3 
DATE: 09 JUL 91 
AUTHOR: LT STEVE GLASER 

COMPILER: MICROSOFT VERSION 4.01 
LINKER: MICROSOFT VERSION 3.55 

SUBROUTINE S3 IS CALLED FROM SUBROUTINE PCBS3. $3 GENERATES THE 
OUTPUT DATA FILE FOR THE THREE COPPER LAYER CASE. 


SUBROUTINE S3 (GYLR,GFXB,GYY,GXX,GZC31,GZC32 ,GZC33,GZB,GYLR31,GYLR 
+32 ,GYLR33,GFXB31,GFXB32,GFXB33 ,GYY31,GYY32 ,GYY33, GXX31,GXX32 ,GXX33 
+,GZE31,GZE32 ,GZE33,GPL,GIH,GJH,GEAT,GWIDE,GDEEP , GOEF ,GIBT, GPRT,GWR 
+T,GRT,GLT,GFT,GBT,GSEL) 


REAL GYLR,GFXB ,GYY,GXX,GZC31,GZB,GYLR31,GFXB31 ,GYY31,GXX31,GZE31,A 
+CC , DAMP ,CONFAC ,GIBT ,GPRT ,GWRT ,GRT ,GLT,GFT,GBT,GYLR32 ,GFXB32 ,GYY32, 
+GXX32 ,GZE32 ,GZC32 ,GYLR33 , GFXB33 , GYY33 , GXX33 ,GZE33 ,GZC33 

REAL GEAT( 100,100) ,GOEF(1000,12) 

CHARACTER*1 SELECT, ANS 

CHARACTER DATAF*79 , NAME*6 


INTEGER GPL,I,N,IB,GWIDE,GDEEP, COUNT, CONTEMP , ZER, GSEL, NMAX, TMAX, HT 
+RS,D1,D2,D3,D4,D5,D6,D7, MAXIT, LOCVAR,IC,ID,IE, IF 


INTEGER GIH( 1000) ,GJH( 1000) ,GCON(1000,12) 


CONTEMP=6 
ZER=0 
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897 


NMAX=750 
TMAX=50 
HTRS=6 

D1=2 

D2=4 

D3=6 

D4=0 

D5=0 

D6=0 

D7=0 
ACC=0. 05 
DAMP=0. 666667 
MAXIT=12 
CONFAC=0. 8 
COUNT=6*GPL 
LOCVAR=COUNT 


CALL CLS 


WRITE(*, 803) 
FORMAT(///,' 
+ INTO THE',/, 


= 


THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY 


EXISTING THERMAL ANALYZER, FURTHERMORE, THIS PROGRAM DOES' ,/, 


NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE’ ,/, 


: THE USER WI 


SHOULD NOT 
PLEASE ENTE 
READ(* ,804) NAME 
FORMAT( A6) 
WRITE(* ,910)NAME 
FORMAT(/////," 

a 3 
WRITE(*,811) 


LL NAME THE DATA FILE FOR EACH RUN OF THIS " ie 


HAVE ANY SPACES. lite 
R THE DESIRED DATA FILE NAME: ‘',2X ) 


+ 
+ 
+ 
+' PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS, AND’ ,/, 
“|. 
+ 


YOU SELECTED ',A6,' FOR YOUR DATA FILE NAME 


FORMAT( IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N 


+FOR NO: Oe 
READ(*,812) ANS 
FORMAT(A1) 


) 


IF(ANS. EQ. 'N') THEN 


GOTO 898 
ELSE 

CONTINUE 
ENDIF 


IFC ANS. EQ. 'Y') THEN 


GOTO 897 


ALLOW THE USER TO PROVIDE TITLE LINE FOR DATA FILE 
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goo) CALL CLS 
WRITE(*, 805) 
805  FORMAT(////,' ENTER THE DESIRED TITLE TO BE PLACED ON LINE 
ais > _ .NUMBER ONE OF THE OUTPUT DATA FILE: ',///, 
a »2X, ) 
READ(*,806) DATAF 
806 FORMAT(A79) 
835 WRITE(*,831) 
831 FORMAT(////,' DO YOU WISH TO CHANGE THE TITLE OF YOUR OUTPU 
+T DATA FILE?',/,' ENTER Y FOR YES AND N FOR NO: m2 xX) 
READ(*,832) ANS 
832 FORMAT(A1) 


IF(ANS. EQ. 'Y') THEN 
GOTO 833 
ELSE 
CONTINUE 
ENDIF 
IF( ANS. EQ. 'N') THEN 
GOTO 834 
ELSE 
GOTO 835 
834 ENDIF 
C 


CARR RIK exe eM I I IK TTT EKER ERE RE REREREREEERERRRERERERERE REE 


Crrdkekieieekek COEFFICIENTS FOR EPOXY AND COPPER LAYERS witisediedeitdesedcdedcdededcy 
Cre tededeetede de Wee deed eM K THI IN TT IIT IIT RITE I KIER ITE TERE TET ITE INTER ITER III KI III 


C 
DO 90 I=1,GPL 
N=] 
IB=GPL+I 
IC=2*GPLt+I 
ID=3*GPL+I 
TE=4*GPLt+I 
IF=5*GPL+I 


C 
CHR KIKI IRIE IE I CORNERS Hie ee aI III I IIIT IIR RII 
C 
IF ((GIH(I). EQ. 1. OR. GIH( 1). EQ. GDEEP). AND. (GJH( 1). EQ. 1. OR. GJH( 1). EQ 
+. GWIDE)) THEN 


C DETERMINE COEFFICIENTS FOR TOP COPPER 
IF (GEAT(GIH(1),GJH(1)).EQ. 0.0) THEN 
GCON(I,N) = 6 
ELSE 
GCON(I,N) = 7 
ENDIF 
C CONNECTIONS FOR EPOXY LAYER 
GCON(IB,N) = 
GCON(ID,N) = 
GCON(IC,N) = 
GCON(ID,N) = 
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€2 02 C2 Gira 


>) 


Gc) 


oa Val Se i ap 


GCON( IE,N) 
GCON( IF ,N) 


ON OV 


LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 


REPT EDGE 

TE S(GJHCI) 280: 1) THEN 
LEFI COEFFICIENT 
GOEF(I,N) = GYLR31 


GOEF(IC,N) = GYLR32 
GOEF(IE,N) = GYLR33 
GOEF(IB,N) = GYLR 
GOEF(ID,N) = GYLR 
GOEF(IF,N) = GYLR 
N=N+1 

GCON(I,N) = 7551 
GCON(IB,N) = 7551 
GCON(IC,N) = 7551 
GCON(ID,N) = 7551 
GCON(IE,N) = 7551 
GCON(IF,N) = 7551 


RIGHT COEFFICIENT 
GOEF(I,N) = GYY31 


GOEF(IC,N) = GYY32 
GOEFCIESN) = Giese 
GOEF(IB,N) = GYY 
GOEF(CID,N) = GYY 
GOEFCIF,N) = GYY 
=N+1 


GCON(CI,N) = 10 * (I+1) + 1 


GCONCIB,N) = 10* (IB+1) + 1 

GCON(CIC,N) = 10 * (IC+1) + 1 

GGONCIBEN) = 10°* Crit) a 

GEONCIE ,.N)-= 100% (inet es 

GCONCIF,N) = 10 * (IF+1) + 1 
RIGHT EDGE 


ELSEIF (GIH(1). EG2GwitE) THEN 
LEFT COEFFICIENT 


GOEF(I,N) = GYY31 


GOEF(IC,N) = GYY32 
GOEF(IE,N) = GYY33 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 

N=N+1 

GCON(I,N) = 10*(I-1)+1 
GCON(IB,N) = 10*(IB-1)+1 
GCON(IC,N) = 10*(IC-1)+1 
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GCON(C ID,N) 20*C1D-1)+1 


GCON(IE,N) = 10*(IE-1)41 
GCON(IF,N) = 10*(IF-1)+1 

c 

C RIGHT COEFFICIENT 

c 
GOEF(I,N) = GYLR31 
GOEF(IC,N) = GYLR32 
GOEF(IE,N) = GYLR33 
GOEF(IB,N) = GYLR 
GOEF(ID,N) = GYLR 
GOEF(IF,N) = GYLR 
N=N+1 
GCON(I,N) = 7541 
GCON(IB,N) = 7541 
GCON(IC,N) = 7541 
GCON(ID,N) = 7541 
GCON(IE,N) = 7541 
GCON(IF,N) = 7541 
ENDIF 

C 

c FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 

c 

C FRONT EDGE 
IF (GIH(I).EQ.1) THEN 

c FRONT COEFFICIENT 


GOEF( I ,N)=GFXB31 
GOEF( IC ,N)=GFXB32 
GOEF( IE,N)=GFXB33 
GOEF(IB,N)=GFXB 
GOEF( ID ,N)=GFXB 
GOEF( EF ,N)=GFXB 


N=N+1 

GCON(I,N) = 7521 
GCON(IB,N) = 7521 
GCON(IC,N) = 7521 
GCON(ID,N) = 7521 
GCON(IE,N) = 7521 


GCON(IF,N) = 7521 

© BACK COEFFICIENT 
GOEF( I ,N)=GXX31 
GOEF( IC ,N)=GXX32 
GOEF(IE,N)=GXX33 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
GOEF( IF ,N)=GXX 
=N+1 
GCON(I,N) =10*(I+GWIDE)+1 
GCON(IB,N) = 10*( IB+GWIDE)+1 


GCON(IC,N) = 10*( IC+GWIDE)+1 
GCON(ID,N) = 10*( ID+GWIDE)+1 
GCON(IE,N) = 10*( IE+GWIDE)+1 
GCON(IF,N) = 10*(IF+GWIDE)+1 
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C BACK EDGE 
ELSEIF (GIH(1I).EQ.GDEEP) THEN 
E FRONT COEFFICIENT 
GOEF( I ,N)=GXX31 
GOEF( IC ,N)=GXX32 
GOEF( IE ,N)=GXX33 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
GOEF( IF ,N)=GXX 


N=N+1 

GCON(CI,N) = 10*( I-GWIDE)+1 

GCON(CIB,N) = 10*( IB-GWIDE)+1 
GCONCIC,N) = 10*( IC-GWIDE)+1 
GCONCID,N) = 10*( ID-GWIDE)+1 
GCONCIE,N) = 10*( IE-GWIDE)+1 
GCON(IF,N) = 10*( IF-GWIDE)+1 


C BACK COEFFICIENT 
GOEF( 1,N)=GFXB31 
GOEF( IC ,N)=GFXB32 
GOEF( IE,N)=GFXB33 
GOEF( IB ,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEF( IF ,N)=GFXB 


N=N+1 

GCON(1,N) = 7511 
GCON(IB,N) = 7511 
GCON(IC,N) = 7511 
GCON(ID,N) = 7511 
GCON(IE,N) = 7511 
GCON(IF,N) = 7511 
ENDIF 


C TOP COEFFICIENT 
GOEF(I,N)=GZE31 
GOEF(IC,N)=GZE32 
GOEF( IE,N)=GZE33 
GOEF(IB,N)=GZC31 
GOEF( ID ,N)=GZC32 
GOEF( IF ,N)=GZC33 


N=N+1 
GCON(I,N) = 7511 
GCON(IC,N) = 7511 
GCON(IE,N) = 7511 
GCON(IB,N) = 10*I+1 
GCON(ID,N) = 10*IC+1 
GCON(IF,N) = 10*IE+1 

C 

® BOTTOM COEFFICIENT 

C 
GOEF(I,N) = GZE31 
GOEF(IC,N) = GZE32 
GOEF(IE,N) = GZE33 
GOEF(IB,N) = GZC31 
GOEF(ID,N) = GZC32 
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a aa 


GOEFCIF,N) = GZB 

N=N+1 

GCON(CI,N) =10*( I+GPL)+1 
GCONCIB,N) = 10*( IB+GPL)+1 
GCONCIC,N) 10*( IC+GPL)+1 
GCONCID,N) 10*( ID+GPL)+1 


C HEAT INPUT 
IF(GEAT(GIH(I),GJH(1)).NE. 0.0) THEN 
GOEF(I,N) = GEAT(GIH(I),GJHCI)) 


N=N+1 
GCON(I,N)=9991 
ENDIF 
C 
iobwmeeGlH( l).E@. 1. OR. GIH(1). EQ. GDEEP). AND. (GJHC1). NE. 1. OR. GJH(1) 
+. NE. GWIDE)) THEN 
C 


C DETERMINE NUMBER OF CONNECTIONS FOR COPPER AND EPOXY LAYERS 
Pee cGEaucGIH(l),GJH(1)).EQ.0.0) THEN 
GCON(I,N) = 6 


ELSE 
GCON(I,N) = 7 

ENDIF 

é 

C CONNECTIONS FOR EPOXY LAYER 
GCON(IB,N) = 6 
GCON(ID,N) = 6 
GCON(IC,N) = 6 
GCON(IE,N) = 6 
GCON(IF,N) = 6 

© 

C 

C LEFT AND RIGHT COEFFICIENTS 

C LEFT COEFFICIENT 
GOEF(I,N) = GYY31 
GOEF(IC,N) = GYY32 
GOEF(IE,N) = GYY33 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 
N=N+1 
GCON(I,N) = 10*(I-1)+1 
GCON(IB,N) = 10*(IB-1)+1 
GCON(IC,N) = 10*(IC-1)+1 
GCON(ID,N) = 10*(ID-1)+1 
GCON(IE,N) = 10*(IE-1)+1 
GCON(IF,N) = 10*(IF-1)+1 

C 


C RIGHT COEFFICIENT 
GOEF(I,N) = GYY31 


GOEFC(IC,N) = GYY32 
GOEFCIE,N) = GYY33 
GOEFCIB,N) = GYY 


me 


GOEF(ID,N) = GYY 
GOEPCIE ND: = Gry 
N=N+1 


GCONCI,N) = 10 * (I+1) + 1 


GCON(IB,N) = 10* (IB+1) + 1 

GCON(IC,N) = 10 * (IC+1) +1 
GCON(ID,N) = 10 * (ID+1) + 1 
GCON(IE,N) = 10 * (IE+1) + 1 
GCON(IF,N) = 10 * (IF+1) +1 


FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 


FRONT EDGE 

IF (GIH(I).EQ.1) THEN 
FRONT COEFFICIENT 
GOEF(1I,N)=GFXB31 
GOEF( IC ,N)=GFXB32 
GOEF( IE,N)=GFXB33 
GOEF(IB,N)=GFXB 

GOEF( ID,N)=GFXB 
GOEF( IF ,N)=GFXB 


© C2 C2 Carta 


N=N+1 

GCON(I,N) =7521 
GCON(IB,N) = 7521 
GCON(IC,N) = 7521 
GCON(ID,N) = 7521 
GCON(IE,N) = 7521 
GCON(IF,N) = 7521 


é BACK COEFFICIENT 
GOEF( I ,N)=GXX31 
GOEF( IC ,N)=Gxx32 
GOEF( IE ,N)=GXX33 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
GOEF( IF ,N)=GXX 


N=N+1 

GCONCI,N) =10*( I+GWIDE)+1 
GCONCIB,N) = 10*( IB+GWIDE)+1 
GCONCIC,N) = 10*( IC+GWIDE)+1 
GCONCID,N) = 10*(C ID+GWIDE)+1 
GCONCIE,N) = 10*( IE+GWIDE)+1 
GCONCIF,N) = 10*( IF+GWIDE)+1 


C BACK EDGE 
ELSEIF (GIH(1).EQ.GDEEP) THEN 
C FRONT COEFFICIENT 
GOEF(I,N)=GXX31 
GOEF( IC ,N)=GXX32 
GOEF( IE ,N)=GXX33 
GOEF(C IB ,N)=GXX 
GOEF (ID ,N)=GXX 
GOEF(C IF ,N)=GXX 
N=N+1 
GCON(I,N) =10*( I-GWIDE)+1 
GCONCIB,N) = 10*( IB-GWIDE)+1 
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GCONC IC ,N) 10*( IC-GWIDE)+1 


GCON(ID,N) = 10*( ID-GWIDE)+1 
GCONCIE,N) = 10*( IE-GWIDE)+1 
GCONCIF,N) = 10*( IF-GWIDE)+1 


C BACK COEFFICIENT 
GOEF(1I,N)=GFXB31 
GOEF(IC,N)=GFXB32 
GOEF( IE,N)=GFXB33 
GOEF(IB,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEF( IF ,N)=GFXB 


N=N+1 

GCON(I,N) = 7511 
GCON(IB,N) = 7511 
GCON(IC,N) = 7511 
GCON(ID,N) = 7511 
GCON(IE,N) = 7511 
GCON(IF,N) = 7511 
ENDIF 


C TOP COEFFICIENT 
GOEF(I,N)=GZE31 
GOEF(IC,N)=GZE32 
GOEF(IE,N)=G2E33 
GOEF(IB,N)=GZC31 
GOEF(ID,N)=GZC32 
GOEF( IF ,N)=GZC33 


N=N+1 
GCONC(I,N) = 7511 
GCON(IC,N) = 7511 
GCONCIB,N) = 10*I+1 
GCON(CID,N) = 10*IC+1 
GCON(IE,N) = 7511 
GCONCIF,N) = 10*IE+1 

| c 

| E BOTTOM COEFFICIENT 

C 

GOEF(I,N) = GZE31 
GOEF(IC,N) = GZE32 
GOEF(IE,N) = GZE33 
GOEF(IB,N) = GZC31 
GOEF(ID,N) = GZC32 
GOEF(IF,N) = GZB 
N=N+1 
GCON(I,N) =10*( I+GPL)+1 
GCON(IB,N) = 10*( IB+GPL)+1 
GCON(IC,N) = 10*( IC+GPL)+1 
GCONCID,N) = 10*( ID+GPL)+1 
GCON(IE,N) = 10*( IE+GPL)+1 


12] 





GCONCIF,N) = 10*( IF+GPL)+1 


C 

C HEAT INPUT 
IF(GEAT(GIH(1),GJH(I)).NE. 0.0) THEN 
GOEF(I,N) = GEAT(GIH(I),GJH(I)) 


N=N+1 
GCON(I,N)=9991 
ENDIF 
CXCXCXC 
CXCXCXG 
€ 
E 
C LEFT AND RIGHT EDGES EXCLUDING CORNERS 
C 
ELSEIF((GIH(I). EQ. 1. OR. GIH(I). EQ. GDEEP). AND. (GJH(I). EQ. 1. OR. GJH(I) 
+.EQ.GWIDE)) THEN 
c 
¢ DETERMINE COEFFICIENTS FOR TOP LAYER 
IF (GEAT(GIH(1I),GJH(I)).NE. 0.0) THEN 
GCONGH Ihe 7 
ELSE 
GCON(I,N) = 6 
ENDIF 
é 
é CONNECTIONS FOR EPOXY LAYER 
GCON(IB,N) = 6 
GCON(IC,N) = 6 
GCON(ID,N) = 6 
GCON(IE,N) = 6 
GCONCIF,N) = 6 
C 
E 
c LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
c 
C LEFT EDGE 
IF (GJH(I).EQ.1) THEN 
C LEFT COEFFICIENT 
GOEF(I,N) = GYLR31 
GOEF(IC,N) = GYLR32 
GOEF(IE,N) = GYLR33 
GOEF(IB,N) = GYLR 
GOEF(ID,N) = GYLR 
GOEF(IF,N) = GYLR 
=N+1 
GCON(I,N) = 7551 
GCON(IB,N) = 7551 
GCON(IC,N) = 7551 
GCON(ID,N) = 7551 
GCONCIE,N) = 7551 
GCON(IF,N) = 7551 
C 


@ RIGHT COEFFICIENT 
GOEF(I,N) = GYY31 
GOEF(IC,N) = GYY32 
GOEF(IE,N) = GYY33 
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GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 
N=N+1 


GCON(CI,N) = 10 * (I+1) + 1 


GCON(IB,N) = 10* (IB+1) + 1 
GCON(IC,N) = 10 * (IC+1) + 1 
GCON(ID,N) = 10 * (ID+1) +1 
GCON(IE,N) = 10 * (IE+1) + 1 
GCON(IF,N) = 10 * (IF+1) + 1 
% 
C RIGHT EDGE 
ELSEIF (GJH(1I).EQ.GWIDE) THEN 
C 
€ LEFT COEFFICIENT 
GOEF(I,N) = GYY31 
GOEF(IC,N) = GYY32 
GOEF(IE,N) = GYY33 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 
N=N+1 
GCON(I,N) = 10*(I-1)+41 
GCON(IB,N) = 10*(IB-1)+1 
GCON(IC,N) = 10*(IC-1)41 
GCON(ID,N) = 10*(ID-1)+1 
GCON(IE,N) = 10*(IE-1)+1 
GCON(IF,N) = 10*(IF-1)+1 
C 
¢ RIGHT COEFFICIENT 
C 
GOEF(I,N) = GYLR31 
GOEF(IC,N) = GYLR32 
GOEF(IE,N) = GYLR33 
GOEF(IB,N) = GYLR 
GOEF(ID,N) = GYLR 
GOEF(IF,N) = GYLR 
N=N+1 
GCON(I,N) = 7541 
GCON(IB,N) = 7541 
GCON(IC,N) = 7541 
GCON(ID,N) = 7541 
GCON(IE,N) = 7541 
GCON(IF,N) = 7541 
ENDIF 
c 


C FRONT COEFFICIENT 
GOEF(I,N) =GXX31 
GOEF( IC ,N)=GXX32 
GOEF( IE,N)=GXX33 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
GOEF( IF ,N)=GXX 
=N+1 
GCON(I,N) =10*(I-GWIDE)+1 
GCON(IB,N) = 10*( IB-GWIDE)+1 
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GCON(CIC,N) = 10*(C IC-GWIDE)+1 
GCON(CID,N) = 10*( ID-GWIDE)+1 
GCONCIE,N) = 10*( IE-GWIDE)+1 
GCONCIF,N) = 10*( IF-GWIDE)+1 

C BACK COEFFICIENT 
GOEF( I ,N)=GXX31 
GOEFCIC,N)=GXXK32 
GOEFCIE,N)=GXX33 
GOEF( IB ,N)=GXX 
GOEF( ID,N)=GXX 
GOEF( IF ,N)=GXX 
N=N+1 
GCON(I,N) =10*( I+GWIDE)+1 
GCON(CIB,N) = 10*( IB+GWIDE)+1 


GCONCIC,N) = 10*( IC+GWIDE)+1 
GCONCID,N) = 10*( ID+GWIDE)+1 
GCON(CIE,N) = 10*( IE+GWIDE)+1 
GCONCIF,N) = 10*( IF+GWIDE)+1 


C DOP CORPPICIE® 
GOEF(I,N)=GZE31 
GOEF( IC ,N)=GZE32 
GOEF( ID,N)=GZE33 
GOEF(IB,N)=GZC31 
GOEF( ID,N)=GZC32 
GOEF(C IF ,N)=GZC33 


=N+1 
GCON(I,N) = 7511 
GeenciGnNy = 7541 
GCON(IB,N) = 10*I+1 
GCON(ID,N) = 10*IC+1 
GCON(IE,N) = 7511 
GCON(IF,N) = 10*IE+1 

C 

C BOTTOM COEFFICIENT 

C 
GOEF(I,N) = GZE31 
GOEF(IC,N) = GZE32 
GOEF(IE,N) = GZE33 
GOEF(IB,N) = GZC31 
GOEF(ID,N) = GZC32 
GOEF(IF,N) = EZB 
N=N+1 
GCON(I,N) = 10*(I+GPL)+1 
GCON(IB,N) = 10*(IB+GPL)+1 
GCON(IC,N) = 10*( IC+GPL)+1 
GCON(ID,N) = 10*( ID+GPL)+1 
GCON(IE,N) = 10%( IE+GPL)+1 
GCON(IF,N) = 10*(IF+GPL)+1 

c 


C HEAT INPUT 
IF(GEAT(GIH(1I) ,GJH(1I)).NE.0.0) THEN 
GOEF(I,N) = GEAT(GIH(1),GJH(I)) 
N=N+1 
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GCON(1I,N)=9991 


ENDIF 
% 
e 
c 
Cy**DETERMINE COEFFICIENTS FOR ALL NODES NOT TOUCHING AN EDGE? tii 
C 
ELSEIF((GIH(I).NE. 1. OR. GIH(1). NE. GDEEP). AND. (GJH(1). NE. 1. OR. GJH(I) 
+.NE.GWIDE)) THEN 
c 
C DETERMINE CONNECTIONS FOR TOP LAYER 
IF (GEAT(GIH(1),GJH(I)).NE.0.0) THEN 
GCON(I,N) = 7 
ELSE 
GCON(I,N) = 6 
ENDIF 
C 
C CONNECTIONS FOR EPOXY LAYER 
GCON(IB,N) = 6 
GCON(IC,N) = 6 
GCON(ID,N) = 6 
GCON(IE,N) = 6 
GCON(IF,N) = 6 
C 
e 
C 
C LEFT COEFFICIENT 
GOEF(I,N) = GYY31 
GOEF(IC,N) = GYY32 
GOEF(IE,N) = GYY33 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 
N=N+1 


ESONCI,N) = 10 * (1-1) + 1 


GCON(IB,N) = 10* (IB-1) + 1 

GCON(IC,N) = 10 * (IC-1) +1 
GCON(ID,N) = 10 * (ID-1) + 1 
GCON(IE,N) = 10 * (IE-1) +1 
GCON(IF,N) = 10 * (IF-1) + 1 


C RIGHT COEFFICIENT 
GOEF(I,N) = GYY31 


GOEF(IC,N) = GYY32 
GOEF(IE,N) = GYY33 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 
N=N+1 


GCON(I,N) = 10 * (I+1) + 1 


GCONCIB,N) = 10* (IB+1) + 1 
GCONCIC,N) = 10 * (IC+1) + 1 
GCONCID,N) = 10 * (ID+1) + 1 
GCONCIE,N) = 10 * (IE+1) + 1 
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GCONCIF,N) = 10 * (IF+1) + 1 


>) 


C FRONT COEFFICIENT 
GOEF(1,N)=GXX31 
GOEF( IC ,N)=GXX32 
GOEF( IE,N)=GXX33 
GOEF( IB ,N)=GXX 
GOEF( ID,N)=GXX 
GOEF( IF ,N)=GXX 
N=N+1 
GCON(I,N) =10*(I-GWIDE)+1 


GCON(IB,N) = 10*( IB-GWIDE)+1 
GCON(IC,N) = 10*( IC-GWIDE)+1 
GCON(ID,N) = 10*( ID-GWIDE)+1 
GCON(IE,N) = 10*( IE-GWIDE)+1 
GCON(IF,N) = 10*( IF-GWIDE)+1 


C BACK COEFFICIENT 
GOEF(I,N)=GXX31 
GOEF( IC ,N)=GXX32 
GOEF(C IE ,N)=GXX33 
GOEF(IB,N)=GXX 
GOEF( ID,N)=GxXxX 
GOEFCIF,N)=GXX 
N=N+1 
GCON(CI,N) =10*( I+GWIDE)+1 
GCONCIB,N) = 10*( IB+GWIDE)+1 


GCONCIC,N) = 10*( IC+GWIDE)+1 
GCONCID,N) = 10*( ID+GWIDE)+1 
GCONCIE,N) = 10*( IE+GWIDE)+1 
GCONCIF,N) = 10*( IF+GWIDE)+1 


C TOP COEFFICIENT 
GOEF(1,N)=GZE31 
GOEF( IC ,N)=GZE32 
GOEF( IE,N)=GZE33 
GOEF(1IB,N)=GZC31 
GOEF( ID,N)=GZ2C32 
GOEF( IF ,N)=GZC33 


N=N+1 
GCONCI,N) = 7511 
GCONCIB,N) = 10*I+1 
GCONCIC,N) = 7511 
GCONCID,N) = 10*IC+1 
GCONCIE,N) = 7511 
GCONCIF,N) = 10*IE+1 

C 

C BOTTOM COEFFICIENT 

C 
GOEF(I,N) = GZE31 
GOEFCIC,N) = GZE32 
GOEFCIE,N) = GZE33 
GOEF(IB,N) = GZC31 
GOEF(CID,N) = GZC32 
GOEFCIF,N) = GZC33 
N=N+1 


GCON(I,N) =10*(I+GPL)+1 
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GCON(IB,N) = 10*( IB+GPL)+1 
GCON(IC,N) = 10*(IC+GPL)+1 
GCON(ID,N) = 10*( ID+GPL)+1 
GCON(IE,N) = 10*( IE+GPL)+1 
GCON(IF,N) = 10*( IF+GPL)+1 


C HEAT INPUT 
IF(GEAT(GIH(1I) ,GJH(I)).NE.0.0) THEN 
GOEF(I,N) = GEAT(GIH(I),GJH(I)) 
N=N+1 
GCON(1I,N)=9991 
ENDIF 
ENDIF 

0 CONTINUE 


9 
C 
C GENERATE DATA FILE 
C 


OPEN (3,FILE=NAME ,FORM='FORMATTED' , ACCESS='DIRECT' , RECL=108 , STATUS 

+='NEW' ) 

WRITE(3,909) DATAF 
909 FORMAT(1X,A79) 

WRITE(3,908) COUNT,CONTEMP, ZER, ZER, ZER, ZER, ZER, ZER,GSEL 
908 FORMAT(2X,9(13,5X)) 

WRITE(3,907) ZER,ZER,ZER 
907. FORMAT(2X,3(13,5X)) 

WRITE(3,9081) NMAX,TMAX,HTRS ,D1,D2,D3,D4,D5 ,D6,D7 
9081 FORMAT(2X,9(13,5X)) 

WRITE(3,905) ACC,DAMP,MAXIT,CONFAC,GIBT 
905 FORMAT(1X,2(F9. 7,1X),14X,12,1X,F9. 7,1X,F9. 5) 

WRITE( 3,906) GPRT,GFT,GBT,GRT,GLT,GWRT 
906  FORMAT(1X,6(F12.3,1X)) 

DO 112 I=1,LOCVAR 

WRITE(3,9100) GCON(I,1),GCON(1I,2),GCON(I,3),GCON(I,4),GCON(1,5),GC 

+ON(1I,6) ,GCON(I,7),GCON(I,8) 
9100 FORMAT(14,3X,7(15,7X)) 

WRITE(3,9110) GOEF(I,1),GOEF(I,2),GOEF(1,3),GOEF(1,4) ,GOEF(I,5),GO 

+EF(1,6),GOEF(I,7) 
9110 FORMAT(7(F9.3,3X)) 
112 CONTINUE 


CLOSE (3) 
CALL CLS 
WRITE(*,999) NAME 
999 FORMAT(////,' THE OUTPUT DATA HAS BEEN PLACED IN A FILE 
+NAMED ',A6,/////,' <PRESS ENTER TO CONTINUE>' ) 


READ(*,5912) ANS 
5912 FORMAT(A1) 


END 
Cretetetetteateak tek teed tet fete Rit KH IC RICHI RIT FRI TTERIREKRE RE RERERERREREEE 


C 

C PIT LE: MODEL BUILDER 

C SUBROUTINE: S4 

C DATE: Os JUL 91 

C AUTHOR: LT STEVE GLASER 

C COMPILER: MICROSOFT VERSION 4.01 

C LINKER: MICROSOFT VERSION 3.55 
C SUBROUTINE S4 IS CALLED FROM SUBROUTINE PCBS3. S84 GENERATES THE 


ey 


C2°G? @) 


898 


803 


OUTPUT DATA FILE FOR THE FOUR COPPER LAYER CASE. 


SUBROUTINE S4 (GYLR,GFXB,GYY,GXX,GZC41,GZC42 ,GZC43 ,GZC44 ,GZB ,GYLR4 
+1,GYLR42 ,GYLR43 ,GYLR44,GFXB41 ,GFXB42 ,GFXB43 ,GFXB44 ,GYY41 ,GYY42,GYY 
+43 ,GYY44 , GXX41 , GXX42 , GXX43 ,GXX44, GZE41 ,GZE42 ,GZE43 ,GZE44,GPL,GIH,G 
+JH,GEAT ,GWIDE ,GDEEP,GOEF ,GIBT,GPRT,GWRT,GRT,GLT,GFT,GBT,GSEL) 


REAL GYLR,GFXB,GYY ,GXX ,GZC41,GZB ,GYLR41,GFXB41,GYY41,GXX41,GZE41,A 
+CC , DAMP , CONFAC ,GIBT ,GPRT ,GWRT,GRT,GLT,GFT,GBT,GYLR42 ,GFXB42 ,GYY42, 
+GXX42 ,GZE42 ,GZC42 , GYLR43 , GFXB43 , GYY43 , GXX43 ,GZE43,GZC43 , GYLR44,GFX 
+B44,GYY44, GXX44,GZE44 ,GZC44 


REAL GEAT( 100,100) ,GOEF( 1000, 12) 


CHARACTER*1 SELECT,ANS 
CHARACTER DATAF*79 , NAME*6 


INTEGER GPL,I,N,1IB,GWIDE ,GDEEP , COUNT , CONTEMP , ZER,GSEL,NMAX , TMAX , HT 
+RS ,D1,D2 ,D3,D4,D5,D6,D7 ,MAXIT ,LOCVAR, IC ID JEV iy tGa, 


INTEGER GIH( 1000) ,GJH( 1000) ,GCON( 1000, 12) 
CONTEMP=6 
ZER=0 
NMAX=750 
TMAX=50 
HTRS=6 

D1=2 

D2=4 

D3=6 

D4=0 

D5=0 

Dé6=0 

D7=0 

ACC=0. 05 
DAMP=0. 666667 
MAXIT=12 
CONFAC=0. 8 
COUNT=8*GPL 
LOCVAR=COUNT 


CALL CLS 


WRITE(* , 803) 

FORMAT(///, ' THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY 
+ INTO THE',/, 

EXISTING THERMAL ANALYZER, FURTHERMORE, THIS PROGRAM DOES' ,/, 
NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE Jn 
THE USER WILL NAME THE DATA FILE FOR EACH RUN OF THIS ’ Je 
PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS, AND ' yy 
SHOULD NOT HAVE ANY SPACES. " Jif 
PLEASE ENTER THE DESIRED DATA FILE NAME: ',2X ) 
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834 
CG 


+ 


+T DATA FILE?',/,' 


READ(* ,804) NAME 
FORMAT( A6 ) 
WRITE(*,910)NAME 


»//) 
WRITE(* 
FORMAT( 


»811) 


READ(*,812) ANS 
FORMAT(A1) 


IF( ANS. EQ. 'N') THEN 
GOTO 898 
ELSE 
CONTINUE 
ENDIF 
IF(ANS. EQ. 'Y') THEN 
GOTO 897 
ELSE 
CATt, CLS 
GOTO 552 
ENDIF 


ALLOW THE USER TO PROVIDE TITLE LINE FOR DATA FILE 


CALL CLS 
WRITE(*, 805) 


FORMAT(///// YOU SELECTED ',A6,' FOR YOUR DATA FILE NAME 


IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N 
+FOR NO: ' poten) 


FORMAT(////," ENTER THE DESIRED TITLE TO BE PLACED ON LINE 


+ 3/; NUMBER ONE OF THE OUTPUT DATA FILE: ',///, 
+ yen) 


READ(*,806) DATAF 
FORMAT( A79 ) 
WRITE(* ,831) 


FORMAT(////,' DO YOU WISH TO CHANGE THE TITLE OF YOUR OUTPU 


READ(*,832) ANS 
FORMAT(A1) 


IF( ANS. EQ. 'Y') THEN 
GOTO 833 

ELSE 
CONTINUE 

ENDIF 

IF( ANS. EQ. 'N') THEN 
GOTO 834 

ELSE 
GOTO 835 
ENDIF 


ENTER Y FOR YES AND N FOR NO: 


' 2X, ) 


Crete dededededede dete Fete ve Fete Fe Fete Fe Fe Fe Fe Fe Fe Fe FE Fe Te FC TCL TE FC TEFEN TC TE PEF E TFET TEFEN TER FOF FETE FRET RR ERREREEEE 


Crkikkireee ke COEFFICIENTS FOR EPOXY AND COPPER LAYERS *tvevediediiiid reir 
Cae dere dete tere tere re HHA IIIT I TI THIN IK RITTER IIR IIR IIIT TI TITIAN 


C 


DO 90 I=1,GPL 


9 


Ciicinkinivinkiokiokiekiiiiekiekicieieiek CORNER Sie serercse tesa aad dese svedesede de ved se dese desea tere seve dese ak a 


C 


ome) 


OOO Go) 


>) 


N=1 
IB=GPL+I 
IC=2*GPL+I 
ID=3*GPL+I 
IE=4*GPL+I 
IF=5*GPL+I 
IG=6*GPL+I 
TH=/7*GPL+I 


IF ((GIH(I). EQ. 1. OR. GIH(I). EQ. GDEEP). AND. (GJH( 1). EQ0 1. OR] Gai Gee 
+, GHIDE)  aHEN 


GCONCI,N) = 6 


ELSE 


GCONCI,N) = 7 


ENDIF 


CONNECTIONS 


GCONCIB,N) 
GCONCID,N) 
GCONCIC,N) 
GCON(CID,N) 
GCONCIE,N) 
GCONC IF ,N) 
GCON(C IG,N) 
GCON( IH,N) 


LEFT EDGE 
IP CGJHe) £OQ. 1) THEN 
LEFT COEERICIENT 


bub wot wt 


R 
6 
6 
6 
6 
6 
6 
6 
6 


GOEF(I,N) = GYLR41 
GOEF(CIC,N) = GYLR42 
GOEF(IE,N) = GYLR43 
GOEF(CIG,N) = GYLR44 
GOEF(IB,N) = GYLR 
GOEF(ID,N) = GYLR 
GOEF(IF,N) = GYLR 
GOEF(IG,N) = GYLR 
GOEF(IH,N) = GYLR 
N=N+1 

GCON(CI,N) = 7551 
GCONCIB,N) = 7551 
GCONCIC,N) = 7551 
GCONCID,N) = 7551 
GCONCIE,N) = 7551 
GGONC TFN) = 7551 
GCONCIG,N) = 7551 


EPOXY LAYER 


DETERMINE COEFFICIENTS FOR TOP COPPER 
IF (GEATC(CGIH(I) ,GJHCI)).EQ. 0.0) THEN 


LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
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C 


gqan oe) 


C 


GCONCIH,N) = 7551 


RIGHT COEFFICIENT 
GOEF(I,N) = GYY41 


GOEF(IC,N) 
GOEF(IE,N) 
GOEF(IG,N) 
GOEF(IB,N) 
GOEF(ID,N) 
GOEF( IF ,N) 
GOEF( IH,N) 
N=N+1 
GCON(I,N) 
GCON( IB,N) 
GCON(IC,N) 
GCON( ID,N) 
GCON( IE,N) 
GCON(IF,N) 
GCON(IG,N) 
GCON( IH,N) 


RIGHT EDGE 


GYY42 
GYY43 
GYY44 
GYY 
GYY 
Gr 
GYY 


OReilly 1 


Oley 1 


10 * (IC+1) 
10 * (ID+1) 
10 * (IE+1) 
10 * (IF+1) 
10 * (IG+1) 
10 * (IH+1) 


ELSEIF (GJH( 1). EQ. GWIDE) 


LEFT COEFFICIENT 


GOEF(I,N) = GYyY41 


GOEF(IC,N) 
GOEF(IE,N) 
GOEF(IG,N) 
GOEF(IB,N) 
GOEF(ID,N) 
GOEF( IF ,N) 
GOEF(IH,N) 
N=N+1 
GCON(I,N) 
GCON( IB,N) 
GCON(IC,N) 
GCON( ID,N) 
GCON(IE,N) 
GCON( IF ,N) 
GCON(IG,N) 
GCON( IH,N) 


GYY42 
GYY43 
GYY44 
ery 
GYY 
Coe 
eNO 


10*(I-1)+1 
10*(IB-1)+1 
10*(IC-1)+1 
10*(ID-1)+1 
10*( IE-1)+1 
10*( IF-1)+1 
10*( IG-1)+1 
10*( IH-1)+1 


RIGHT COEFFICIENT 


GOEF(I,N) 
GOEF(IC,N) 
GOEF(IE,N) 
GOEF(IG,N) 
GOEF(1IB,N) 
GOEF( ID ,N) 
GOEF( IF ,N) 
GOEF(IH,N) 


GYLR41 
GYLR42 
GYLR43 
GYLR44 
GYLR 
GYLR 
GYLR 
GYLR 


t++t+ett 
a oe oe oe oll 


THEN 
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N=N+1 
GCON(I,N) = 7541 


GCON(IB,N) = 7541 
GCON(IC,N) = 7541 
GCON(ID,N) = 7541 
GCON(IE,N) = 7541 
GCON(IF,N) = 7541 
GCON(IG,N) = 7541 
GCON(IH,N) = 7541 
ENDIF 


FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 


FRONT EDGE 
IF (GIH(1).EQ. 1) THEN 

FRONT COEFFICIENT 
GOEF(1,N)=GFXB41 
GOEF( IC ,N)=GFXB42 
GOEF( IE,N)=GFXB43 
GOEF(1G,N)=GFXB44 
GOEF(1B,N)=GFXB 
GOEF( ID ,N)=GFXB 
GOEF( IF ,N)=GFXB 
GOEF( 1H ,N)=GFXB 


>) anaan 


N=N+1 

GCONCI,N) = 7521 
GCONCIB,N) = 7521 
GCONCIC,N) = 7521 
GCONCID,N) = 7521 
GCONCIE,N) = 7521 
GCONCIF,N) = 7521 
GCONCIG,N) = 7521 
GCONCIH,N) = 7521 


C BACK COEFFICIENT 
GOEF( 1,N)=GXX41 
GOEF( IC ,N)=GXX42 
GOEF( IE ,N)=GXX43 
GOEF( IG,N)=GXX44 
GOEF( 1B,N)=GXX 
GOEF( ID,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF( IH,N)=GXX 
N=N+1 
GCON(1,N) =10*( I+GWIDE)+1 
GCON(IB,N) = 10*( IB+GWIDE)+1 


GCONCIC,N) = 10*( IC+GWIDE)+1 
GCONCID,N) = 10*( ID+GWIDE)+1 
GCONCIE,N) = 10*( IE+GWIDE)+1 
GCONCIF,N) = 10*(C IF+GWIDE)+1 


c BACK EDGE 
ELSEIF (GIH(1).EQ.GDEEP) THEN 
C FRONT COEFFICIENT 
GOEF(1I,N)=GXX41 
GOEF(1IC,N)=GXX42 


13z 


C 


C 


GOEF( IE ,N)=GXX43 
GOEF( IG ,N)=GXX44 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF( IH,N)=GXX 
=N+1 
GCON(I,N) 
GCON( IB ,N) 
GCON(IC,N) 


GCONCIE,N) 


GCONC IF ,N) 
GCONCIG,N) 


GOEF(I,N)=GFXB41 


GOEF( IC ,N)=GFXB42 
GOEF( IE ,N)=GFXB43 
GOEF( IG,N)=GFXB44 


GOEF( IB,N)=GFXB 
GOEF( ID,N)=GFXB 
GOEF( IF ,N)=GFXB 
GOEF( IH,N)=GFXB 
N=N+1 

GCON(I,N) = 7511 


GCON(IB,N) = 7511 
GCON(IC,N) = 7511 
GCON(ID,N) = 7511 
GCON(IE,N) = 7511 
GCON(IF,N) = 7511 
GCON(IG,N) = 7511 
GCON(IH,N) = 7511 
ENDIF 


TOP COEFFICIENT 


GOEF(1I,N)=GZE41 

GOEF(IC,N)=GZE42 
GOEF( IE ,N)=GZE43 
GOEF(IG,N)=GZE44 
GOEF( IB ,N)=GZC41 
GOEF( ID,N)=G2C42 
GOEF( IF ,N)=G2C43 
GOEF(IH,N)=GZC44 
N=N+1 

GCON(I,N) = 7511 


GCON(IC,N) = 7511 
GCON(IE,N) = 7511 
GCON(IB,N) = 10*I+1 
GCON(ID,N) = 
GCON(IF,N) = 
GCON(IG,N) = 7511 
GCON(IH,N) = 


BOTTOM COEFFICIENT 


= 10*( I-GWIDE)+1 
= 10*( IB-GWIDE)+1 
= 10*( IC-GWIDE)+1 
GCONCID,N) = 10*( ID-GWIDE)+1 
= 10*( IE-GWIDE)+1 
= 10*( IF-GWIDE)+1 
= 10*( IG-GWIDE)+1 
GCONCIH,N) = 10*( IH-GWIDE)+1 
BACK COEFFICIENT 


1071 C+ 1 
PO TE+ 


10*IF+1 
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SP Jal Sp I 


GOEFCI,N) = GZE41 


GOEF(IC,N) = GZE42 
GOEF(IE,N) = GZE43 
GOEF(IG,N) = GZE44 
GOEF(IB,N) = GZC41 
GOEF(ID,N) = GZC42 
GOEF(IF,N) = GZC43 
GOEF(IH,N) = GZB 
N=N+1 


GCON(I,N) =10*(1+GPL)+1 
GCON(IB,N) = 10%*( IB+GPL)+1 


GCON(IC,N) = 10*( IC+GPL)+1 

GCON(ID,N) = 10*( ID+GPL)+1 

GCON(IE,N) = 10*( IE+GPL)+1 

GCON(IF,N) = 10*(IF+GPL)+1 

GCON(IG,N) = 10*( IG+GPL)+1 

GCON(IH,N) = 10*( IH+GPL)+1 
HEAT INPUT 


IF(GEAT(GIH(I),GJH(1)).NE. 0.0) THEN 
GOEF(I,N) = GEAT(GIH(1I),GJH(I)) 
N=N+1 

GCON(I,N)=9991 

ENDIF 


ELSEIF((GIH(1). EQ. 1. OR. GIH(1). EQ. GDEEP). AND. (GJH(1). NE. 1) OR) GUHGas 
+. NE. GWIDE)) THEN 


DETERMINE NUMBER OF CONNECTIONS FOR COPPER AND EPOXY LAYERS 
IF (GEAT(GIH(1),GJH(I)).EQ. 0.0) THEN 
GCON(I,N) = 6 
ELSE 
GCON(I,N) = 7 
ENDIF 


CONNECTIONS 
GCON( IB ,N) 
GCON( ID,N) 
GCON(IC,N) 
GCON( IE,N) 
GCON( IF ,N) 
GCON(IG,N) 
GCON( IH,N) 


OR EPOXKY “LAYER 


now uu wm wou 


DHDHAAAN DH DW 


LEFT AND RIGHT COEFFICIENTS 
LEFT COEFFICIENT 


GOEF(I,N) = GYY41 
GOEF(IC,N) = GYY42 
GOEF(IE,N) = GYY43 
GOEF(IG,N) = GYY44 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 
GOEF(IH,N) = GYY 
N=N+1 


GCONCI,N) = 10*(I-1)+1 
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C 


>) amaaa 


C 


GCON( IB,N) 
GCON(IC,N) 
GCON(ID,N) 
GCON( IE,N) 
GCON( IF ,N) 
GCON( IG,N) 
GCON( IH,N) 


GOEF(I,N) 
GOEF(IC,N) 
GOEF(IE,N) 
GOEF(IG,N) 
GOEF( IB,N) 
GOEF(ID,N) 
GOEF(IF,N) 
GOEF(IH,N) 
N=N+1 
GCON(I,N) 
GCON( IB ,N) 
GCON(IC,N) 
GCON( ID,N) 
GCON( IE,N) 
GCON( IF,N) 
GCON( IG,N) 
GCON( IH,N) 


FRONT EDGE 


rund nu it de a 


— 
— 


nud du eu 


pe GiB =] )+1 
Oreo = 1) +1 
0s CID-1)+1 
Oe Ee = 1) +1 
EO“CIr=1 )-+1 
HOrC1G-1)+1 
10*( IH-1)+1 


RIGHT COEFFICIENT 


GYY41 
GYY42 
GYY43 
GYY44 
Grr 
Gry 
GYY 
Giy 


Ome itl) + 1 


hob] ) + 1 


10 * (IC+1) 
10 * (ID+1) 
10 * (IE+1) 
10 * (IF+1) 
10 * (IG+1) 
10 * (IH+1) 


mee (GIHCI).EQ. 1) THEN 


FRONT COEFFICIENT 


GOEF(I,N)=GFXB41 
GOEF( IC ,N)=GFXB42 
GOEF( IE ,N)=GFXB43 
GOEF( IG,N)=GFXB44 
GOEF( IB ,N)=GFXB 
GOEF( ID ,N)=GFXB 
GOEF( IF ,N)=GFXB 
GOEF( IH,N)=GFXB 


N=N+1 


GCON(I,N) =7521 


GCON( IB,N) 
GCON(IC,N) 
GCON(ID,N) 
GCON( IE,N) 
GCON( IF ,N) 
GCON(IG,N) 
GCON( IH,N) 


yal 
USall 
fo2l 
2 
fs Pall 
LS) 
1 Saal 


BACK COEFFICIENT 


GOEF( I ,N)=GXX41 
GOEF( IC ,N)=GXX42 
GOEF( IE ,N)=GXX43 
GOEF(IG,N)=GXX44 
GOEF( IB,N)=GXX 


t+++tt 
a ee 


FRONT AND BACK COEFFICIENTS 


DEPENDING ON WHICH EDGE 


15 


C 
C 


C 


GOEF( 1ID,N)=GXX 
GOEFC IF ,N)=GXX 
GOEF( IH,N)=GXX 


N=N+1 

GCONCI,N) =10*( I+GWIDE)+1 
GCONCIB,N) = 10*( IB+GWIDE)+1 
GCONCIC,N) = 10*( IC+GWIDE)+1 
GCONCID,N) = 10*( ID+GWIDE)+1 
GCONCIE,N) = 10*( IE+GWIDE)+1 
GCONCIF,N) = 10*( IF+GWIDE)+1 
GCONCIG,N) = 10*( IG+GWIDE)+1 
GCONCIH,N) = 10*( IH+GWIDE)+1 
BACK EDGE 


ELSEIF (GIH(1).EQ.GDEEP) THEN 


FRONT COEFFICIENT 


GOEF(CI,N)=GXX41 
GOEFCIC,N)=GXX42 
GOEF( IE ,N)=GXX43 
GOEF(IG,N)=GXxX44 


GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF( IH ,N)=GXX 


=N+1 


GCON(I,N) =10*( I-GWIDE)+1 


GCON( IB,N) 
GCON(IC,N) 
GCON(ID,N) 
GCON(IE,N) 
GCON(IF,N) 
GCON(IG,N) 
GCON( IH,N) 


10*( IB-GWIDE)+1 
10*( IC-GWIDE)+1 
10*( ID-GWIDE)+1 
10*( IE-GWIDE)+1 
10*( IF-GWIDE)+1 
10*( IE-GWIDE)+1 
10*( IF-GWIDE)+1 


BACK COEFFICIENT 


GOEF(I,N)=GFXB41 

GOEF( IC ,N)=GFXB42 
GOEF( IE ,N)=GFXB43 
GOEF( IG,N)=GFXB44 


GOEF( 1B,N)=GFXB 
GOEF( ID,N)=GFXB 
GOEF( IF ,N)=GFXB 
GOEF( IH,N)=GFXB 


=N+1 


GCONCI,N) = 7511 


GCON(IB,N) = 7511 
GCON(IC,N) = 7511 
GCON(ID,N) = 7511 
GCON(IE,N) = 7511 
GCON(IF,N) = 7511 
GCON(IG,N) = 7511 
GCON(IH,N) = 7511 
ENDIF 


TOP COEFFICIENT 


GOEF(1I,N)=GZE41 
GOEF( IC ,N)=GZE42 
GOEF( IE,N)=GZE43 
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GOEF(IG,N)=GZE44 
GOEF(1IB,N)=GZC41 
GOEF(ID,N)=GZC42 
GOEF( IF ,N)=GZC43 
GOEF( IH,N)=GZC44 


N=N+1 
GCON(I,N) = 7511 
GCON(IC,N) = 7511 
GCON(IB,N) = 10*I+1 
GCON(ID,N) = 10*IC+1 
GCON(IE,N) = 7511 
GCON(IF,N) = 10*IE+1 
GCON(IG,N) = 7511 
GCON(IH,N) = 10*IG+1 

C 

C BOTTOM COEFFICIENT 

C 
GOEF(I,N) = GZE41 
GOEF(IC,N) = GZE42 
GOEF(IE,N) = GZE43 
GOEF(IG,N) = GZE44 
GOEF(IB,N) = GZC41 
GOEF(ID,N) = GZC42 
GOEF(IF,N) = GZC43 
GOEF(IH,N) = GZB 
N=N+1 
GCON(I,N) =10*(I+GPL)+1 
GCON(IB,N) = 10*(IB+GPL)+1 
GCON(IC,N) = 10*( IC+GPL)+1 
GCON(ID,N) = 10*( ID+GPL)+1 
GCON(IE,N) = 10%*( IE+GPL)+1 
GCON(IF,N) = 10*(IF+GPL)+1 
GCON(IG,N) = 10*( IG+GPL)+1 
GCON(IH,N) = 10*( IH+GPL)+1 

C 


C HEAT INPUT 
PaCGEATCGIH(1),.GJHC1)).NE.0.0) THEN 
GOEF(I,N) = GEAT(GIH(I),GJH(I)) 

=N+1 
GCON(I,N)=9991 
ENDIF 


CXCXCXC 

CXCXCXC 

C 

C 

C LEFT AND RIGHT EDGES EXCLUDING CORNERS 
C 


ELSEIF((GIH(1). EQ. 1. OR. GIH(1). EQ. GDEEP). AND. (GJH(1). EQ. 1. OR. GJH(I) 
+.EQ.GWIDE)) THEN 


an 


DETERMINE COEFFICIENTS FOR TOP LAYER 
IF (GEAT(GIH(1),GJH(1)).NE. 0.0) THEN 
GCON(I,N) = 7 
ELSE 
GCON(I,N) = 6 


Log 


ENDIF 


C CONNECTIONS 
GCON(IB,N) 
GCON(IC,N) 
GCON(ID,N) 
GCON(IE,N) 
GCON( IF,N) 
GCON(IG,N) 
GCON(IH,N) 


OR EPOXY LAYER 


hun dw wow ws 


R 
6 
6 
6 
6 
6 
6 
6 


LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 


LEFT EDGE 

IE (GJ CT). Boe) THEN 
LEFT COEFFICIENT 
GOEF(I,N) = GYLR41 


Q C2 C2) Gao) 


GOEF(IC,N) = GYLR42 
GOEF(IE,N) = GYLR43 
GOEF(IG,N) = GYLR44 
GOEF(IB,N) = GYLR 
GOEF(ID,N) = GYLR 
GOEF(IF,N) = GYLR 
GOEF(IH,N) = GYLR 
=N+1 

GCON(I,N) = 7551 
GCON(IB,N) = 7551 
GCON(IC,N) = 7551 
GCON(ID,N) = 7551 
GCON(IE,N) = 7551 
GCON(IF,N) = 7551 
GCON(IG,N) = 7551 
GCON(IH,N) = 7551 


C | RIGHT COEFFICIENT 
GOEF(I,N) = GYY41 


GOEF(IC,N) = GYY42 
GOEFCIE,N) = GYY43 
GOEFCIG,N) = GYY44 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
GOEFCIF,N) = GYY 
GOEF(IH,N) = GYY 
N=N+1 


GEONCICN) = 10 Clee 1 


GCON(IB,N) = 10* (IB+1) + 1 
GCON(IC,N) = 10 * (IC+1) + 1 
GCON(ID,N) = 10 * (ID+1) + 1 
GCON(IE,N) = 10 * (IE+1) + 1 
GCON(IF,N) = 10 * (IF+1) + 1 
GCON(IG,N) = 10 * (IGt1) + 1 
GCON(IH,N) = 10 * (IH+1) + 1 


C RIGHT EDGE 
ELSEIF (GJH( 1). EQ. GWADE) 


: 


C LEFT COEFFICIENT 
GOEF(I,N) = GYY41 
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C 


C 


C 


GOEF(IC,N) 
GOEF(IE,N) 
GOEF(IG,N) 
GOEF(IB,N) 
GOEF(ID,N) 
GOEF(C IF ,N) 
GOEF(IH,N) 
N=N+1 
GCONCI,N) 
GCONC(CIB,N) 
GCONCIC,N) 
GCONCID,N) 
GCONCIE,N) 
GCONC IF ,N) 
GCONCIG,N) 
GCON(IH,N) 


GOEF(1I,N) 
GOEF(IC,N) 
GOEF(IE,N) 
GOEF(IG,N) 
GOEF(IB,N) 
GOEF(ID,N) 
GOEF(IF,N) 
GOEF( IH,N) 
N=N+1 
GCON(I,N) 
GCON(IB,N) 
GCON(IC,N) 
GCON(ID,N) 
GCON( IE,N) 
GCON( IF ,N) 
GCON(IG,N) 
GCON(IH,N) 
ENDIF 


GYY42 
GYY43 
GYY44 
GYY 
GYY 
GYY 
GYY 


10*(I-1)+1 

10*( IB-1)+1 
10*(IC-1)+1 
10*(ID-1)+1 
10*( IE-1)+1 
10*( IF-1)+1 
10*( IG-1)+1 
10*( IH-1)+1 


RIGHT COEFFICIENT 


GYLR41 
GYLR42 
GYLR43 
GYLR44 
GYLR 
GYLR 
GYLR 
GYLR 


7541 
7541 
7541 
7541 
7541 
7541 
7541 
7541 


FRONT COEFFICIENT 


GOEF(I,N) =GXX41 
GOEF( IC ,N)=GXX42 
GOEF( IE ,N)=GXX43 
GOEF(IG,N)=GXX44 
GOEF( 1B ,N)=GXX 
GOEF(ID,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF( IH,N)=GXX 


N=N+1 


GCONCI,N) =10*( I-GWIDE)+1 


GCON( IB,N) 
GCON(IC,N) 
GCON(ID,N) 
GCON(IE,N) 
GCON( IF ,N) 
GCON(IG,N) 
GCON(IH,N) 


10*( IB-GWIDE)+1 
10*( IC-GWIDE)+1 
HO*CTD-GWIDE )+1 
10*( IE-GWIDE)+1 
10*( IF-GWIDE)+1 
10*( IG-GWIDE)+1 
10**( IH-GWIDE)+1 


BACK COEFFICIENT 
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C 


GOEF(1I,N)=GXX41 
GOEF( IC ,N)=GXX42 
GOEF( IE ,N)=GXX43 
GOEF(1IG,N)=GXX44 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF( IH ,N)=GXX 
N=N+1 


GCON(CI,N) =10*( I+GWIDE)+1 


GCON(IB,N) = 10*( IB+GWIDE)+1 
GCON(IC,N) = 10*( IC+GWIDE)+1 
GCON(ID,N) = 10*( ID+GWIDE)+1 
GCON(IE,N) = 10%*( IE+GWIDE)+1 
GCON(IF,N) = 10%*( IF+GWIDE)+1 
GCON(IG,N) = 10*( IG+GWIDE)+1 
GCON(IH,N) = 10%*( IH+GWIDE)+1 


TOP COEFFICIENT 


GOEF(I,N)=GZE41 

GOEF(IC,N)=GZE42 
GOEF(IE,N)=GZE43 
GOEF(IG,N)=GZE44 
GOEF(IB,N)=GZC41 
GOEF(ID,N)=GZC42 
GOEF( IF ,N)=GZC43 
GOEF(IH,N)=GZC44 
N=N+1 

GCON(I,N) = 7511 


GCON(IC,N) = 7511 
GCON(IB,N) = 10*I+1 
GCON(ID,N) = 10*IC+1 
GCON(IE,N) = 7511 
GCON(IF,N) = 10*IE+1 
GCON(IG,N) = 7511 
GCON(IH,N) = 10*IG+1 


BOTTOM COEFFICIENT 


GOEF(I,N) = GZE41 
GOEF(IC,N) = GZE42 
GOEF(IE,N) = GZE43 
GOEF(IG,N) = GZE44 
GOEF(IB,N) = GZC43 
GOEF(ID,N) = GZC43 
GOEF(IF,N) = GZC43 
GOEF(IH,N) = EZB 

N=N+1 

GCON(I,N) = 10*(I+GPL)+1 
GCON(IB,N) = 10*(IB+GPL)+1 
GCON(IC,N) = 10*(IC+GPL)+1 
GCON(ID,N) = 10*(ID+GPL)+1 
GCON(IE,N) = 10*(IE+GPL)+1 
GCON(IF,N) = 10*( IF+GPL)+1 
GCON(IG,N) = 10*( IG+GPL)+1 
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GCONCIH,N) = 10*( IH+GPL)+1 


oO 


C HEAT INPUT 
IF(GEAT(GIH(I) ,GJH(I)).NE.0.0) THEN 
GOEF(I,N) = GEAT(GIH( 1) ,GJH(I)) 
N=N+1 
GCON(CI,N)=9991 
ENDIF 
C 
Cx***DETERMINE COEFFICIENTS FOR ALL NODES NOT TOUCHING AN EDGE***eerx 
C 
Pee lhCcGluiCtje Ne. 1.OR. GIH( 1). NE. GDEEP). AND. (GJH(1).NE. 1. OR. GJH(I) 
+. NE. GWIDE)) THEN 


C DETERMINE CONNECTIONS FOR TOP LAYER 
IF (GEAT(GIH(1),GJH(1I)).NE.0.0) THEN 
GCON(I,N) = 7 


SE 
GCON(I,N) = 6 


ENDIF 
C 
C CONNECTIONS FOR EPOXY LAYER 
GCON(IB,N) = 6 
GCON(IC,N) = 6 
GCON(ID,N) = 6 
GCON(IE,N) = 6 
GCON(IF,N) = 6 
GCON(IG,N) = 6 
GCON(IH,N) = 6 
e 
C 
C LEFT COEFFICIENT 
GOEF(1,N) = GYY41 
GOEF(IC,N) = GYY42 
GOEF(IE,N) = GYY43 
GOEF(IG,N) = GYY44 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 
GOEF(IH,N) = GYY 
N=N+1 
GCON(I,N) = 10 * (I-1) + 1 
GCON(IB,N) = 10* (IB-1) + 1 
GCON(IC,N) = 10 * (iC-1) + 1 
GCON(ID,N) = 10 * (ID-1) + 1 
GCON(IE,N) = 10 * (IE-1) + 1 
GCON(IF,N) = 10 * (IF-1) + 1 
GCON(IG,N) = 10 * (IG-1) + 1 
GCON(IH,N) = 10 * (IH-1) + 1 
C RIGHT COEFFICIENT 
GOEF(I,N) = GYY41 
GOEF(IC,N) = GYY42 
GOEF(IE,N) = GYY43 
GOEF(IG,N) = GYY44 
GOEF(IB,N) = GYY 
GOEF(ID,N) = GYY 
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C 


C 


C 


GOEF(CIF,N) = GYY 

GOEF(IH,N) = GYY 

N=N+1 

GCONCI,N) = 10 * (I+1) + 1 
GCONCIB,N) = 10* CIB+1) + 1 
GCON(CIC,N) = 10 * (IC+1) + 1 
GCONCID,N) = 10 * (ID+1) + 1 
GCONCIE,N) = 10 * (CIE+1) + 1 
GCONCIF,N) = 10 * (IF+1) + 1 
GCONCIG,N) = 10 * (IG+1) + 1 
GCONCIH,N) = 10 * (IH+1) + 1 


FRONT COEFFICIENT 


GOEF( I ,N)=GXX41 
GOEF( IC ,N)=GXX42 
GOEF( IE ,N)=GXX43 
GOEFCIG,N)=GXX44 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
GOEF(C IF ,N)=GxXxX 
GOEF(IH,N)=GXX 
N=N+1 


GCON(I,N) = 


GCON(IB,N) 
GCON( IC,N) 
GCON(ID,N) 
GCON( IE,N) 
GCON( IF ,N) 
GCON(IG,N) 
GCON( 1H,N) 


O*( I -GWIDE)+1 


1 


10*( IB-GWIDE)+1 
10*( IC-GWIDE)+1 
10*( ID-GWIDE)+1 
fe CIE -GW IE 
10*( IF-GWIDE)+1 
10*( IG-GWIDE)+1 
10*( IH-GWIDE)+1 


BACK COEFFICIENT 


GOEF(I,N)=GXX41 
GOEF( IC ,N)=GXX42 
GOEF( IE ,N)=GXX43 
GOEF( 1H,N)=GXX44 
GOEF( 1B ,N)=GXX 
GOEF( ID ,N)=GXX 


GOEFC IF ,N)=GxXxX 
GOEF( IH,N)=GXX 


N=N+1 

GCONCI,N) =10*( I+GWIDE)+1 
GCON( IB, N) = 10*( IB+GWIDE)+1 
GCONCIC,N) = 10*( IC+GWIDE)+1 
GCON( ID,N) = 10*( ID+GWIDE)+1 
GCONCIE,N) = 10*( IE+GWIDE)+1 
GCONCIF,N) = 10%*( IF+GWIDE)+1 
GCONCIG,N) = 10*( IG+GWIDE)+1 
GCON(IH,N) = 10*( IH+GWIDE)+1 


TCRSCOERPICTENT 


GOEF(I,N)=GZE41 

GOEFC IC ,N)=GZE42 
GOEFC( IE ,N)=GZE43 
GOEF(IG,N)=GZE44 
GOEF(IB,N)=GZC41 
GOEFCID,N)=GZ2C42 
GOEFCIF,N)=GZC43 
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GOEF( IH,N)=GZC44 


N=N+1 
GCON(I,N) = 7511 
GCON(IB,N) = 10*I+1 
GCON(IC,N) = 7511 
GCON(ID,N) = 10*IC+1 
GCON(IE,N) = 7511 
GCON(IF,N) = 10*IE+1 
GCON(IG,N) = 7511 
GCON(IH,N) = 10*IG+1 

G 

¢ BOTTOM COEFFICIENT 

C 
GOEF(I,N) = GZE41 
GOEF(IC,N) = GZE42 
GOEF(IE,N) = GZE43 
GOEF(IG,N) = GZE44 
GOEF(IB,N) = GZC41 
GOEFCID,N) = GZC42 
GOEF(IF,N) = GZC43 
GOEF(IH,N) = EZB 
=N+1 
GCON(I,N) =10*( I+GPL)+1 
GCON(IB,N) = 10*(IB+GPL)+1 
GCON(IC,N) = 10*( IC+GPL)+1 
GCON(ID,N) = 10%*( ID+GPL)+1 
GCON(IE,N) = 10*( IE+GPL)+1 
GCON(IF,N) = 10%*( IF+GPL)+1 
GCON(IG,N) = 10%*(IG+GPL)+1 
GCON(IH,N) = 10*( IH+GPL)+1 


C HEAT INPUT 
IF(GEAT(GIH(1) ,GJH(1)).NE.0.0) THEN 
GOEF(I,N) = GEAT(GIH(1) ,GJH(I)) 
N=N+1 
GCON(I,N)=9991 
ENDIF 
ENDIF 

0 CONTINUE 


° 
C 
C GENERATE DATA FILE 
C 


OPEN (3, FILE=NAME ,FORM=' FORMATTED’ , ACCESS=' DIRECT’ , RECL=108 , STATUS 
+='NEW' ) 
WRITE(3,909) DATAF 
909 FORMAT(1X,A79) 
WRITE(3,908) COUNT,CONTEMP, ZER, ZER, ZER, ZER, ZER, ZER,GSEL 
908 FORMAT(2X,9(13,5X)) 
WRITE(3,907) ZER,ZER,ZER 
907. FORMAT(2X,3(13,5X)) 
WRITE(3,9081) NMAX, TMAX,HTRS,D1,D2,D3,D4,D5 ,D6,D7 
9081 FORMAT(2X,9(13,5X)) 
WRITE(3,905) ACC,DAMP ,MAXIT,CONFAC ,GIBT 
905  FORMAT(1X,2(F9.7,1X),14X,12,1X,F9. 7,1X,F9.5) 
WRITE(3,906) GPRT,GFT,GBT,GRT,GLT,GWRT 
906 FORMAT(1X,6(F12. 3,1X)) 
DO 112 I=1,LOCVAR 
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WRITE(3,9100) GCON(1,1),GCON(I,2),GCON(1,3),GCON(I,4),GCON(I,5),GC 
+ON(I,6),GCON(1,7),GCON(I,8) 

9100 FORMAT(I4,3X,7(15,7X)) 
WRITE(3,9110) GOEF(1I,1),GOEF(1,2),GOEF(1,3),GOEF(1,4),GOEF(I,5),GO 
+EF(1,6),GOEF(1I,7) 

9110 FORMAT(7(F9.3,3X)) 

112 CONTINUE 


CLOSE? (3) 
CALL CLS 
WRITE(*,999) NAME 
999 FORMAT(////,' THE OUTPUT DATA HAS BEEN PLACED IN A FILE 
+NAMED | ,A6,/////,° <PRESS ENTER TO CONTINUE>' ) 


READ(*,5912) ANS 
5912 FORMAT(A1) 


END 
CC RERKERKEKKEKKEKKERKRERERREERERERERERRERERKRRERERERERERERRERERERREERRERERERREREREE 
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